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I.  INTRODUCTIW 


CIMAG2  is  a  user  interactive  FORTRAN  program  which  can  be  used  to 
produce  the  color  images  of  the  targets  using  either  measured  or 
simulated  scattered  field  data.  The  program  is  also  capable  of 
processing  either  frequency  or  time  domain  data  prior  to  the  image 
processing. 

The  original  version  of  the  program  is  written  by  Dr.  J.  Young 
and  the  new  version  is  written  and  modified  by  S.  Smithberger  and  G. 
Dural.  An  access  for  the  'file  read'  routines  of  the  program  'FTRAN' 

[11  is  provided  for  the  users  who  already  processed  their  data  by  FTRAN. 

Color  images  are  displayed  on  a  Tektronix  429  color  CRT  display 
and  hard  copy  unit  via  the  computer  program  CLRPL  (Appendix  D).  CIMAG2 
accesses  the  program  CLRPL  via  the  'DCL  commands'  facility  of  the 
program. 

Theory  related  to  the  imaging  procedure  is  described  in  detail  in 
another  report  entitled  "Polarimetric  ISAR  Imaging  Using  Either  Measured 
or  Calculated  Transient  Signatures"  [2]. 

Chapter  II  contains  the  User's  Manual  of  the  program.  The 
Programmer's  Manual  is  included  in  Chapter  III.  Chapter  IV  covers  the 
conclusions.  Program  listings  and  a  list  for  the  'HELP'  library  are 
contained  in  the  appendices. 
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II.  CIMAG2  USER'S  MANUAL 

This  manual  shows  the  user  how  to  use  the  CIMAG2  program.  It  can 
be  read  in  bits  and  pieces  as  the  user  finds  a  need  for  information  but 
it  is  recommended  that  the  user  take  the  time  to  sit  down  and  read  all 
of  the  information  given  here.  CIMAG2  is  a  very  sophisticated  program 
into  which  a  lot  of  automation  has  been  installed.  The  basic  routines 
can  be  used  to  get  perfectly  usable  image,  but  the  user  will  be  able  to 
handle  more  data  more  efficiently  if  he  takes  the  time  to  learn  the 
automatic  features  of  the  program.  Not  only  will  he  save  time  in 
generating  the  images  but  he  will  also  find  it  easier  to  keep  track  of 
what  he  has  done.  It  is  also  recommended  that  the  user  read  this  manual 
while  sitting  in  front  of  a  terminal  in  order  to  get  a  better  grasp  and 
feel  for  what  is  going  on.  There  are  many  examples  of  command  sequences 
in  this  manual  which  can  be  used  for  illustrative  purposes. 

This  manual  stresses  processes  in  which  a  series  of  commands  have 
to  be  used  rather  than  the  individual  commands  themselves.  If  the  user 
needs  more  information  about  the  commands  available,  he  can  do  the 
following: 

-  get  into  the  program  by  typing 

$RUN  user2: (DURAL. CIMAG2)CIMAG2 

-  type  the  command 

OHELP 

At  this  point  the  computer  will  list  a  menu  of  commands  used  by  the 
CIMAG2  program.  The  user  simply  types  the  command  that  he  is  interested 
in  and  the  computer  will  give  him  the  information  that  he  seeks.  When 
he  is  finished  he  simply  keeps  hitting  the  <RETURN>  until  he  gets  back 


to  the  'O'  prompt.  This  command  is  most  useful  when  one  is  in  the 
middle  of  doing  something  and  can' t  remember  what  a  command  is  called  or 
can't  remember  what  a  certain  command  does. 

2.1  PUTTING  THE  DATA  TOGETHER 
2.1.1  Appending  Data  Files 

Some  of  the  data  used  by  the  CIMAG2  program  will  come  in  chunks. 
In  other  words,  sometimes  the  data  for  a  given  look  angle  couldn't  be 
taken  all  at  once  so  it  was  necessary  to  make  more  than  one  data  file  in 
order  to  get  all  the  available  data.  If  we  want  to  use  all  of  the  data 
that  is  available  or  at  least  more  than  one  data  file  for  a  given  look 
angle  then  we  will  have  to  append  the  data  files  together.  This  seems 
simple  enough  since  a  given  element  of  the  array  is  assigned  a  certain 
frequency.  However,  many  times  the  files  may  overlap.  One  file  may 
contain  data  from  1.5-6. 5  GHz  and  the  other  may  contain  data  from  6-12 
GHz.  When  we  put  the  two  files  together  we  don't  want  the  overlapping 
part  of  the  files  to  add  up  so  we  will  have  to  use  what  is  called  the 
FGT  command  which  puts  a  trapezoidal  gate  around  the  information  that  we 
are  trying  to  get  to  and  in  effect  chops  off  the  information  that  we 
don't  want.  After  the  files  have  been  properly  gated  then  we  will  use 
the  COM  (combine)  command  to  put  the  files  together.  The  following  is  a 
general  procedure  for  putting  two  files  together.  The  text  following  the 
exclamation  points  are  comments  to  inform  the  reader  what  is  going  on. 
The  exclamation  marks  and  comments  are  not  to  be  typed  into  the  program. 


ORDFL 


OSBF 

BUFFER  #  1 
ORDFL 


I  A  frequency  domain  read  command 
!  This  command  will  be  covered  in 
!  the  following  sections. 

!  The  store  in  buffer  command 


OSBF 

BUFFER  #  2 


!  A  good  way  of  deciding  which 
!  elements  of  the  array  should 
!  be  modified  is  through  the  use 
!  use  of  the  TYP  command.  This 
I  command  will  list  any  range 
!  of  the  array  elements  and 
!  their  respective  values. 


ORBl 


OFGT 


! 


HARMONIC  FOR  START  OF  HIGH  FREQUENCY  CUTOFF 
598 

HIGH  FREQUENCY  ROLL-OFF  IN  DB  PER  HARMONIC 
50.0 


HARMONIC  FOR  START  OF  LOW  FREQUENCY  CUTOFF 
1 


LOW  FREQ  ROLLOFF  IN  DB 

0.01 


Put  the  contents  of  buffer  1  in 
the  main  buffer. 

This  is  the  gating  command 

The  array  element 
where  the  cutoff 
begins. 

This  is  quite  a 
bit  since  we 
want 

a  sharp  cutoff. 
This  element  is 
in  the  region 
that  is  never 
used  close 
to  DC. 

PER  HARMONIC  !  A  very  slow  roll-off 


OSB3 


I  It  is  stored  in  buffer  3  in  so 
!  that  if  an  error  was  made  we 
!  can  try  again. 


! 

We  will  use  the  combine  command 

I 

to  put  the  two  files  together 

! 

but  the  files  have  to  be 

f 

converted  to  the  time  domain 

ORB2 

1 

before  we  can  do  that. 

OIFF 

! 

Inverse  Fast  Fourier  transform 

1 

does  this  conversion. 

OCOM 

STORAGE  BUFFER  NUMBER (USE  0  TO  FINISH):  4  !  Combine  the  two 

MULTIPLIER  =1.0  I  files. 

STORAGE  BUFFER  NUMBER(USE  0  TO  FINISH);  5 
MULTIPLIER  =  1.0 

STORAGE  BUFFER  NUMBER(USE  0  TO  FINISH):  0  !  A  "0"  means  we're 

!  finished. 

!  At  this  point  it  is  recommended 
!  to  plot  the  waveform  and  see 
!  whether  it  is  continuous.  We 
!  don't  want  a  spike  or  null  to 
I  occur  at  the  point  we  combine 
!  two  files. 

OSB6 


5 


I 


$ 


ORB6 

OFFT 

OSB7 


!  At  this  point  we  should  check 
!  the  frequency  domain  again  to 
!  make  sure  we  did  the  gating 
!  right. 

!  If  data  overlapped  then  there 
!  will  be  a  spike  where  the  two 
!  files  were  joined.  If  there 
!  was  a  hole  left  in  the  data 
!  then  a  null  will  appear  where 
I  the  two  files  were  joined. 


I  Fast  Fourier  Transform  converts 
!  the  time  domain  back  to  the 
I  frequency  domain. 


ORB6 

OWRI 


I  This  command  is  pretty  well 
!  self  explanatory. 


2.1.2  Processing  Routines  Nov  In  Use 

These  are  the  processes  that  are  currently  being  used  on 
the  data  before  an  image  Is  generated. 

2. 1.2.1  Read  Comaands 

There  are  four  read  commands  of  the  two  types  of  reads  to 
be  performed  by  this  program.  The  data  that  is  taken  from  the 
radar  range  is  stored  in  a  frequency  domain  format.  This  is  the 
raw  data.  For  imaging  purposes  we  will  be  using  the  (filtered) 
time  domain  waveform.  CIMAG2  can  only  store  time  domain 
waveforms.  Most  of  the  time  the  user  will  process  a  frequency 
domain  waveform,  multiply  it  with  frequency  transform  it  to  the 
time  domain,  and  then  store  it.  Later  we  will  come  back  and  read 
in  several  of  these  time  domain  waveforms  and  use  them  to  make  an 
image.  With  this  procedure  in  mind  we  will  now  explain  the  four 
read  commands  —  two  of  which  (FTREA,FFREA)  are  used  to  read  the 
files  processed  by  FTRAN  |1J. 

RDFL  Frequency  Domain  Read 

The  RDFL  command  was  designed  to  read  in  the  calibrated 
frequency  domain  files  from  the  ESL  database.  The  following  are 
questions  the  that  the  program  will  ask  and  what  they  mean  to  the 
user: 

FREQUENCY  SAMPLING(l)  OR  . IKL  SAMPLING(O)? 

FREQUENCY  INCREMENT  IN  MHZ 


1 


If  the  user  chooses  to  use  frequency  sampling  then  the 
program  will  ask  for  the  frequency  increment  to  be  used  for  each 
sample. 

INPUT  MAJOR  AXIS  DIMENSION  IN  INCHES 


If  the  user  chooses  to  use  .IKL  sampling  then  CIMAG2  will 
ask  this  question. 

SELECT  THE  TYPE  OF  INTERPOLATION 

INPUT  0  — >  TWO-POINT  INTERPOLATION  ;  NO  SMOOTHING 

INPUT  1  — >  INTERPOLATION  AND  SMOOTHING  USING  A  COSINE  WINDOW 

Most  of  the  time  there  will  be  no  need  for  smoothing 
since  the  user  will  more  than  likely  be  using  data  that  has 
already  been  calibrated.  Usually  the  data  is  smoothed  in  the 
calibration  phase  but  if  the  data  has  not  been  sufficiently 
smoothed,  then  the  facilities  are  here  to  do  more  smoothing. 

ASSUMED  INPUT  AMPLITUDE  IS  IN  DB/SQUARE  CM 
NORMALIZE  TO:  SQ  CM(1),SQ  M(2) ,PI*L*L/A(0)? 

TYPE  DATA  FILE  NAME 

This  is  obviously  the  file  name  of  the  data.  If  the  data 
file  is  not  in  the  current  directory  be  sure  to  fully  specify  the 
file  name.  The  program  will  then  type  out  the  header  of  the  file 
so  that  the  user  can  be  sure  he  has  the  right  file. 
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REA  Tiae  Doaain  Read 


The  REA  command  is  a  much  simpler  command.  It  simply 
asks  for  the  file  name  and  then  types  out  the  header.  This 
conunand  is  only  intended  to  read  files  that  were  output  from  the 
CIMAG2  program  therefore  it  has  been  taken  for  granted  that  the 
file  was  stored  in  whatever  form  it  was  needed.  However,  the 
file  can  be  processed  more  after  it  has  been  read  in. 

FTREA  or  FFREA 

These  are  the  time  and  Frequency  domain  r'cr'd  commands  for 
the  files  processed  by  the  program  FTRAN  [1].  Since  FTRAN  is 
capable  of  reading  either  750  or  11/23  type  data,  the  user  must 
be  careful  about  the  format  of  the  data  file.  Default  is  11/23. 
The  FTY  command  can  be  used  to  switch  from  one  format  to  another. 

FTY 

To  control  the  format  of  the  data  file  the  user  should 
type  the  command  FTY.  The  program  then  will  ask  about  the  file 
type.  Enter  T  for  750,  and  F  for  11/23  type.  Default  is  11/23 
type  when  the  program  is  started  unless  no  FTY  command  is  used 
then. 


2. 1.2.2  IFF  Inverse  Fourier  Transform 
The  INVERSE  FOURIER  TRANSFORM  command  converts  a  data 


file  from  the  frequency  domain  to  the  time  domain.  This  command 


<v*v  r-s^r  w: 


must  be  performed  before  the  file  can  be  used  for  an  image  and 


before  the  file  can  be  output  using  the  WRI  command. 


2. 1.2.3  ROT  Rotate 


The  ROTATE  command  takes  all  the  elements  of  a  time 


domain  file  and  moves  them  either  in  the  positive  direction  or 


the  negative  direction.  When  an  Inverse  Fourier  Transform  is 


used  to  convert  a  waveform  from  the  frequency  domain  to  the  time 


domain  a  waveform  is  created  which  repeats  throughout  time.  The 


4096  element  array  is  in  effect  a  time  window  which  shows  one 


complete  Iteration  of  this  waveform.  The  main  pulse  of  a  time 


domain  waveform  may  occur  around  the  zero  time  position.  When 


this  happens  part  of  the  pulse  will  be  plotted  at  the  beginning 


of  the  plot  and  part  of  the  pulse  will  be  plotted  at  the  end  of 


the  plot.  Using  this  command  we  can  move  the  pulse  to  some  other 


part  of  the  plot  (usually  to  element  1024)  so  that  the  entire 


pulse  can  be  plotted  in  one  area  of  the  plot. 


2. 1.2.4  VND  Windowing 


The  frequency  domain  files  that  we  are  using  for  data  are 


band  limited.  On  either  side  of  the  valid  data  that  lies  within 


this  bandwidth  is  a  null  value.  This  creates  a  very  sharp  change 


in  the  amplitude  of  the  waveform  where  the  data  starts  and  ends. 


The  technique  of  Fourier  Transform  assumes  that  the  frequency 


domain  waveform  is  continuous.  The  sharp  changes  in  the  data 


file  are  very  discontinuous.  Thus  when  these  discontinuities  are 


I 


transformed  into  the  time  domain  they  cause  an  oscillation  of  the 
time  domain  waveform.  The  WINDOW  command  reduces  this  problem  by 
convolving  the  frequency  domain  waveform  with  a  cosine.  This 
greatly  reduces  the  change  in  amplitude  that  occurs  at  the 
beginning  and  the  end  of  the  valid  data  thereby  reducing  the 
oscillation  that  results  in  the  time  domain. 

2. 1.2.5  DCV  Dotmconversion 

Since  the  bandwidth  of  the  data  that  we  are  using  lies  in 
a  region  well  above  DC,  many  of  the  resulting  time  domain 
waveforms  take  the  appearance  of  a  modulated  waveform  much  like 
that  received  from  a  radio  station.  To  get  rid  of  this  effect  we 
are  currently  doing  just  what  the  radio  receivers  do:  move  the 
center  of  the  bandwidth  to  DC.  That  is  what  the  DOWNCONVERT 
command  does.  It  respond  with  the  number  that  represents  the 
middle  of  the  bandwidth. 

eg.  : 

If  we  are  using  a  bandwidth  of  1-12  GHz,  then  the 
center  of  the  bandwidth  is  6.5  GHz.  ((12-l)/2  +1)  This 
corresponds  to  element  650. 

2. 1.2. 6  MJW  Multiply  by  jw 

Multiplying  by  jw  in  the  frequency  domain  is  the  same  as 
differentiation  in  the  time  domain.  Another  way  of  thinking  of 
it  is  to  say  that  it  will  shift  the  phase  of  the  waveform  ninety 


degrees  and  act  as  a  high  pass  filter.  One  can  see  that  the 
higher  the  frequency  the  higher  the  amplitude. 

2. 1.2. 7  NVR  Multiply  by  v 

This  command  is  used  as  a  high  pass  filter.  The  higher 
the  frequency  is  the  the  bigger  the  multiplying  factor  will  be. 
It  is  a  required  step  in  data  processing  with  the  algorithm 
described  in  [2]. 

2. 1.2.8  VRI  Vrite 

The  VRITE  command  is  the  only  way  to  output  a  file.  Many  times 
after  the  user  has  processed  a  file  he  will  want  to  save  it.  Then  when 
he  wants  to  use  it  again  he  simply  reads  it  in  and  the  file  can  be  used 
in  a  pre-processed  form.  This  command  can  only  write  out  time  domain 
files  so  somewhere  along  the  line  the  file  will  have  had  to  been 
transformed  using  the  IFF  command.  If  the  user  wishes  to  save  the  file 
and  then  later  when  he  reads  it  in  perform  some  more  processing  in  the 
frequency  domain,  he  simply  reads  it  in  in  the  time  domain  and  then 
performs  a  FFT  command  which  will  convert  it  back  to  the  frequency 
domain. 

2.1.3  Summary  Of  Data  Processing 

Here  is  the  sequence  of  commands  that  we  are  currently 
using  to  process  the  data: 


1  RDFL 


I 

% 

s 

*•. 

f 

I 

I 


2  VND 


3  MWR 


4  DCV 


5  IFF 


-  frequency  sampling  every  ten  degrees 

(Different  sampling  intervals 
can  be  used  to  scale  the 
image. ) 

-  no  smoothing 

-  normalized  to  square  meters 


-  use  a  Hanning  Window 


6  ROT  (Optional) 

-  usually  by  1024 

7  WRI 


2.2  FROCBOURE  OBFINITKmS 

2.2.1  What  Are  Procedure  Definitions? 

Most  of  the  time  when  an  image  is  made  many  data  files 
will  be  needed.  This  means  that  many  data  files  will  have  to  be 
processed  using  the  same  command  sequences  with  the  same 
parameters.  This  is  time  consuming,  redundant,  costly,  and 
bothersome.  CIMAG2  has  a  way  to  get  around  this  problem.  A 
procedure  definition  is  a  series  of  commands  that  are  to  be 
performed  on  a  number  of  data  files.  It  will  take  a  list  of 
input  files  and  a  list  of  output  files,  process  each  file  and 
then  put  the  results  in  the  corresponding  output  file.  This 
allows  the  user  to  go  through  the  procedure  once  and  then  let  the 
computer  do  all  the  work. 


2.2.2  Hov  to  Use  the  PROC  Coaaand 

The  PROC  command  is  a  very  powerful  command  but  the  user 
must  take  great  care  in  using  it.  It  is  a  good  idea  to  go 
through  the  first  data  file  the  regular  way  to  make  sure  the 
method  of  data  processing  gives  the  results  that  you  are  looking 
for.  Then  you  can  go  ahead  and  define  a  procedure  with  an  input 
list  and  an  output  list. 

COMPARE  THE  RESULTS  WITH  THE  FILE  MADE  THE  REGULAR  WAY 

TO  MAKE  SURE  THE  PROCEDURE  IS  DOING  WHAT  YOU  THINK  IT  IS  DOING!! 

It  is  very  possible  to  think  you  have  defined  one  thing 
when  in  fact  you  have  defined  something  else.  This  may  take  some 
time  but  not  as  much  time  as  processing  all  the  files  manually  or 
redoing  an  image  because  you  have  used  data  that  is  trash. 

Let's  say  that  you  have  already  done  the  test  case  and 
you  know  the  exact  sequence  of  commands  that  you  want  to  perform 
on  each  data  file.  This  is  the  sequence  of  events  that  will 
occur  when  you  define  the  procedure: 


OPROC 


>5ia 


f 


CM 


Do  you  have  a  procedure  definition  file 
for  this  process  already?(y  or  N) 


!  If  you  answer  this  question  with  a  Y  then  you 
!  will  be  asked  for  a  filename.  If  you  answer 
!  N  then  you  will  be  given  the  following  prompt. 

Enter  the  process  using  regular  commands  and 
NAHE.DAT  for  a  filename.  For  the  new 
filename  use  NEVNAME.DAT.  When  finished 
defining  the  process  use  the  command  DONE. 

(Warning:  a  filename  must  be  listed  for 
each  time  it  is  used.) 


!  At  this  point  the  user  will  type  in  something  that 
!  looks  like  the  following.  Notice  how  the  filenames 
1  NAHE.DAT  and  NEVNAME.DAT  have  been  used.  The  lines 
!  that  the  user  input  have  been  marked  with  a 


ORDFL 

FREQUENCY  SAMPLING(l)  OR  .IKL  SAMPLING (0)? 

1 

FREQUENCY  INCREMENT  IN  MHZ 

10. 

SELECT  THE  TYPE  OF  INTERPOLATION 
INPUT  0  — >  TWO-POINT  INTERPOLATION  ;  NO  SMOOTHING 
INPUT  1  — >  INTERPOLATION  AND  SMOOTHING  USING  A  COSINE  WINDOW 
0 

ASSUMED  INPUT  AMPLITUDE  IS  IN  DB/SQUARE  CM 
NORMALIZE  TO:  SQ  CM(1),SQ  M(2) ,PI*L*L/A(0)? 

2 

TYPE  DATA  FILE  NAME 
NAHE.DAT 

Dummy  file  for  procedure  definition 
a  6  in.  sphere 

NL1200  FF=  1000  IN=  9  frequency  domain 
OWND 

INPUT  HARMONICS; ST ART, END, TYPE  OF  WINDOW 
TYPE:0=HANNING,1=HAMMING,2=GAUSSIAN,  0,N,M,=TEST 
100,1200,0 
ODCV 

INPUT  THE  HARMONIC  NUMBER  TO  BE  MOVED  TO  DC 


fk  * 

650 

u  * 

OMWR 

”  * 

OIFF 

★ 

OROT 

ROTATE  BY  INCREMENTS  OF: 

* 

1024 

•k 

m 

OWRI 

FILE  NAME  ? 

f 

NEWNAME.DAT 

«  a 

DO  YOU  WANT  TO  CHANGE  THE  FILE  HEADER  Y=l,  N-0 

* 

0 

t**  * 

ODONE 

V!. 

Do  you  wish  to  save  this  procedure  de£inition?(Y 

k 

Y 

Filename: 

^  * 

PROC.DEF 

Do  you  have  a  data  list  file?{Y  or  N) 

k 

S 

N 

Enter  the  list  of  data  files,  following 
each  with  <RETURN>.  When  finished  type  the 

word  DONE. 

* 

DATA1.DAT 

★ 

DATA2.DAT 

i  * 

DATA3.DAT 

DONE 

Do  you  wish  to  save  this  data  list?(Y  or  N) 

★ 

Y 

>" 

y 

Filename: 

y  * 

INPUT.DAT 

Is  there  an  output  filename  list?(Y  or  N) 

■ 

N 

i 

Enter  a  list  of  the  output  file  names  in 
the  order  they  are  to  be  used.  When 

finished  type  DONE. 

★ 

0UT1.DAT 

b  * 

0UT2.DAT 

0UT3.DAT 

&  * 

DONE 

Do  you  wish  to  save  this  list?(Y  or  N) 

X-  * 

Y 

★ 

Filename: 

0UTPUT.DAT 

Your  data  is  being  processed. 

I  At  this  point  the  computer  processes  the  data 
1  according  to  what  you  have  told  it  to  do. 

!  CHECK  THE  OUTPUT  TO  MAKE  SURE  IT  IS  WHAT  YOU 
^  !  WANTED! ! 


is 


SI 


2.3  CREATING  AN  IMAGE 

The  main  objective  of  this  program  is  to  produce  an  image 
of  the  target  on  a  computer  monitor.  The  idea  is  that  we  can  get 
this  image  to  the  point  where  we  can  identify  the  arget  with  the 
image.  This  is  the  sequence  of  events  used  to  for.n  an  image  on 
the  screen: 


-  read  in  and  store  all  the  data  files  needed 

for  the  image  in  the  buffers 

-  CMI  command 

-  IMG  command 

-  $RUN  USER2: [DURAL. IMAGE ICLRPL 


2.3.1  Reading  In  the  Data 

Assuming  that  the  data  has  been  processed  using  some 
method  based  on  methods  presented  earlier  in  this  manual  filtered 
and  stored  in  the  time  domain,  we  can  use  a  REA  (or  PTREA  command 
if  data  are  not  processed  by  CIHAG2)  to  read  in  each  data  file 
and  then  we  can  use  the  SBF  command  to  store  all  the  data  files 
in  separate  buffers.  The  SHO_BUF  command  can  be  used  to  give  a 
listing  of  all  the  buffers  and  their  contents.  However  the  user 
will  have  to  remember  the  polarization,  buffer#,  look  angle,  and 
center  element  for  each  data  file  so  he  will  probably  want  to 
keep  track  of  these  things  on  a  sheet  of  paper  while  going 
through  this  process. 


m 


2.3.2  The  CHI  Couand 


The  CMI  command  sets  up  the  data  structures  for  the  rest 
of  the  imaging  processes.  First  the  program  will  ask  the  user 
for  the  number  of  files  in  a  given  polarization.  It  will  look 
like  this: 

NUMBER  OF  W  TIME  DOMAIN  WAVEFORMS  TO  BE  USED? 

W  stands  for  vertical  polarization.  The  user  should 
answer  with  an  integer.  Then  it  will  ask  the  user  these  questions 
for  each  of  the  files  for  that  given  polarization: 

BUFFER  NUMBER  FOR  W  FILE  # 

LOOK  ANGLE  IN  DEGREES  FOR  FILE  # 

CENTER  ELEMENT  NUMBER  FOR  FILE  # 

The  center  element  of  the  file  will  be  zero  unless  the 
user  has  used  the  rotate  command  on  the  data.  After  all  the 
information  for  a  given  polarization  has  been  accumulatedi  then 
CIMAG2  will  repeat  the  sequence  for  HH  (horizontal)  and  HV  (cross 
polarization) . 

2.3.3  The  IMG  Cowand 

The  IMG  combines  all  the  data  down  into  a  hundred  by 
hundred  matrix.  These  are  the  questions  that  will  be  asked: 


SIZE  OF  THE  IMAGE, (1  TO  4096)=? 

Generally  a  good  value  for  this  is  around  300.  The  user 
shouldn't  really  use  any  bigger  value  than  this  since  the  current 
resolution  is  only  hundred  by  hundred.  When  a  higher  resolution 
device  is  connected  up  to  this  software,  it  might  be  better  to 
get  a  larger  window. 

LOOK  ANGLE  OF  THE  IMAGE, DEGREES? 

The  computer  is  able  to  spin  the  image  that  is  on  the 
screen  so  that  it  may  be  easier  to  see  certain  things  but  this 
doesn't  really  have  any  real  affect  on  the  image. 

POLARIZATION  OF  THE  IMAGE,  (1=W,2=HH,3=HV,4=ALL) 

This  is  self  explanatory. 

After  this  command  has  finished  executing  then  the  final  array  is 
ready  to  be  imaged.  After  the  image  is  generated  the  computer 
asks, 

DO  YOU  WANT  TO  STORE  THE  IMAGE?  Y=1,N=0 

If  the  answer  is  "1"  then  the  computer  asks  for  the  name  of  the 
data  file  to  store  the  image  and  the  frequency  increment  for  the 
frequency  domain  signal  (usually  10  MHz)  which  is  used  for 
calculating  the  time  axes  in  the  plot. 


2.3.4  OUTPUT  THE  IMAGE 

In  order  to  output  the  image  the  user  should  enter 
$RIIN  USER2:IDURAL.IMAGE]CLRPL 

2.4  THE  LOG  &  FILE  COMMANDS 

2.4.1  Creating:  the  History  of  An  laage 

The  sequence  of  events  necessary  to  create  a  given  image 
can  be  long  and  complex.  Many  times  it  is  advantageous  to  keep 
track  of  exactly  what  has  been  done  to  a  file  or  an  image.  This 
can  be  especially  helpful  if  an  error  occurs  or  you  want  to 
reevaluate  a  couple  of  different  procedures.  The  history  of  the 
procedure  can  also  be  used  to  regenerate  results.  This  is 
especially  useful  in  the  case  of  images.  The  image  can  be 
reproduced  much  faster  using  automated  techniques  rather  than 
manually  entering  all  the  commands  that  are  necessary.  The  final 
reason  for  generation  of  a  history  is  that  it  is  easy  to  do  and 
if  the  user  can  take  advantage  of  a  history  with  little  effort 
then  why  not  do  it? 

To  make  the  history  is  a  simple  matter.  The  user  simply 
types  the  command  LOG  before  he  types  in  the  procedure.  CIMAG2 
will  prompt  for  a  filename  in  which  to  dump  the  history.  Then 
the  user  uses  the  program  the  same  way  he  would  if  there  were  no 
history  being  kept.  When  he  is  finished  with  the  procedure  he 
types  the  command  STOLOG. 

This  closes  the  history  file  given  by  the  LOG  command. 


2.4.2  Regenerating  Results  (the  PILE  coaaand) 

As  was  mentioned  in  section  5.1  the  history  of  a 
procedure  can  be  be  used  to  regenerate  the  results  of  that 
procedure.  To  do  this  the  user  will  use  the  FILE  command.  The 
user  then  types  in  the  name  of  the  history  and  the  computer  will 
take  over  and  perform  the  procedure  stored  in  the  history  file. 
When  the  computer  is  finished  executing  this  file  it  will  give 
the  user  this  message: 

CONTROL  HAS  RETURNED  TO  THE  TERMINAL 

At  this  point  the  user  is  free  again  to  do  as  he  wishes. 

If  the  user  would  like  to  see  the  normal  prompts  that  the 
computer  outputs  to  the  screen  for  each  command  so  that  he  can 
follow  the  execution  of  the  history  file  he  can  use  the  ECHO 
command  before  he  uses  the  FILE  command  and  the  computer  will 
output  this  information. 


2.5  QUICK  REFERENCE  FOR  SOME  COMMANDS  OF  GENERAL  USE 

CLR  -  clears  the  main  buffer 

CLRBUF  -  clears  all  the  buffers 

CMI  -  initialize  buffers  for  imaging 

COM  -  combine  command  adds  time  domain  waveforms  together 
DCV  -  downconversion 

DEF  -  sets  the  default  directory  for  input  files 

DEL  -  deletes  a  file  from  the  directory 
ECHO  -  sets  the  echo  for  the  FILE  command 
FFREA  -  FTRAN  Read  Command  (Frequency) 

FFT  -  fast  fourier  transform 

FGT  -  trapezoidal  gate 

FILE  -  executes  a  history  file 

FTREA  -  FTRAN  Read  Command  (Time) 

FTY  -  File  type  for  the  input  format 

IFF  -  inverse  fourier  transform 

IMG  -  sets  up  the  overall  parameters  for  image 

LOG  -  creates  a  log  or  history  file 

MJU  -  multiply  by  jw 

N0_ECH0  -  stops  the  echo  for  the  FILE  command 
PLO  -  plot  command  plots  a  data  file  on  one  of  the  plotting 
devices 

PROC  -  procedure  definition 

RBF  -  read  buffer 

RDFL  -  read  frequency  domain  file 

REA  -  read  time  domain  file 

RLB  -  relabel  the  file  header 

ROT  -  rotate  a  time  domain  file 
SBF  -  store  in  buffer 

SH0_BUF  -  sho  the  contents  of  all  the  buffers 
STO_LOG  -  stop  the  log  or  history 

TYP  -  type  the  value  of  all  the  elements  of  a  data  file 
WND  -  window 

VfRI  -  write  time  domain  data  file 


CIHAG2  PROGRAMMER'S  MANUAL 


CIMAG2  is  a  very  large,  non-trivial  program.  It  also  has 
some  distinct  methods  of  control.  It  is  the  aim  of  this  manual 
to  familiarize  the  user  with  the  methods  and  data  structures  used 
in  CIMAG2.  If  more  information  is  needed  the  Vax  Fortran  Manual 
will  have  more  detailed  information  on  file  handling.  The 
purpose  of  this  manual  is  to  inform  the  reader  of  names  and 
structures  unique  to  this  program.  It  is  suggested  that  the 
programmer  read  the  CIMAG2  User's  Manual  first.  This  will 
familiarize  the  programmer  with  the  commands  that  have  been 
implemented  and  give  him  background  for  the  discussions  that 
follow. 

3.1  Linking  CIMAG2 

CIMAG2  is  a  very  complicated  program.  It  is  divided  up 
into  many  smaller  programs  and  data  files.  The  data  files  not 
only  include  data  files  created  by  users  but  also  data  files  for 
such  structures  as  error  messages,  help  libraries  etc.  For  some 
of  these  files  it  is  necessary  to  link  them  with  the  program  but 
others  must  be  linked.  As  a  result  the  linking  for  CIMAG2  has 
become  rather  cumbersome.  To  get  around  this  the  linking  for  the 
program  has  been  placed  in  the  command  file 

USER2: (DURAL. CIMAG2 ]CIMAG2. COM.  This  means  that  there  are  only 
two  commands  necessary  to  convert  a  new  version  of  CIMAG2.F0R 
into  an  execution  file: 


$FOR/LIST  CIMAG2 


$@CIMAG2 

The  first  command  compiles  the  main  program  and  creates  a 
listing  file.  The  second  command  executes  the  linking  command 
file.  A  listing  of  CIMAG2.COM  is  given  in  Section  (3-4).  It  is 
recommended  that  programmers  first  make  changes  to  the  NEVMAG.FOR 
file  first.  Then  after  the  new  routine  has  been  debugged  and 
tested  NEVMAG.FOR  can  be  copied  over  into  CIMAG2.F0R.  This  means 
that  a  working  version  of  the  program  will  always  be  in  CIMAG2. 
NEVMAG  can  be  linked  the  same  way  that  CIMAG2  is  through  the  use 
of  user2: [DURAL. CIMAG2]NEVMAG. COM.  This  command  file  links 
NEVMAG  in  the  same  a  way  shown  below. 

$F0R/LIST  NEVMAG 
$@NEVMAG 

A  listing  of  NEVNAG.COM  can  also  be  found  in  Section  (3-4). 

3.2  Process  Control  Structure 

The  Vax  treats  files  and  devices  the  same  way.  In  other 
words  the  input  from  the  terminal  looks  like  a  file  to  the  Vax. 
This  fact  been  used  in  a  few  of  the  most  powerful  commands  in 
CIMAG2;  namely  the  FILE,  PROC,  LOG,  and  ST0_L0G  commands.  Ve  can 
transfer  control  from  the  terminal  to  an  input  file  and  back 
again.  This  gives  the  user  the  ability  to  use  processes  that 


have  already  been  defined  and  define  new  processes  himself.  Here 
are  the  basics  of  how  this  process  control  is  implemented  in 
CIMAG2. 

All  devices  and  files  that  are  to  be  used  are  assigned 
vhat  is  called  Logical  Unit  Numbers  (LU#).  Each  file  or  device 
will  be  referenced  by  this  number  in  read  and  write  statements. 

(  When  a  programmer  writes  new  software  for  CIMAG2  he  should  not 
use  any  ACCEPT  or  TYPE  statements  in  his  code.  These  two 
statements  will  defeat  the  purpose  of  this  whole  method  of 
control.)  Here  are  the  three  basic  control  structures: 

C0M_UNIT  ,  !  LU#  t<  r  the  command  input 

lOUT  ,  I  LU#  for  the  program  output 

LOGUNIT  ,  !  LU#  for  the  log  file 

COM_UNIT  and  lOUT  can  be  set  to  either  the  terminal  or  a 
file.  The  L0G_UNIT  can  either  be  set  to  a  file  or  to  the  null 
device.  Using  these  three  structures  the  I/O  for  a  given  routine 
would  be  written  like  this: 

READ  (COMUNIT,*)  X 

WRITE  (LOG_UNIT,*)  X 

WRITE  (lOUT,*)  X 

Notice  how  the  input  was  immediately  written  to  the 
LOGUNIT.  All  input  should  be  done  this  way.  This  allows  the 
routine  to  be  used  in  procedure  definitions  and  logging  files. 

The  format  was  used  in  this  example  but  any  numbered  format 
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may  be  used  just  as  in  any  Fortran  programming.  For  these 
structures  to  be  used  in  this  vay  they  must  be  assigned  values 
that  correspond  to  the  various  devices  and  files.  Here  are  some 
variable  names  that  are  assigned: 

TERM_UNIT  ,  !  LU#  for  the  terminal 

FILE_IINIT  ,  !  LU#  for  the  command  file 

NULL_UNIT,  !  LU#  of  null  device 

STO_UNIT  !  LU#  for  buffer  storage 

Using  these  variables  this  is  the  way  that  the  structures 
are  initialised: 

C0M_UNIT  =  TERM_UNIT 
lOUT  =  TERM  UNIT 
LOGUNIT  =  NULL_UNIT 

Initially  the  input  is  coming  from  the  terminal,  the 
output  is  going  to  the  terminal,  and  we  are  logging  to  the  null 
unit  which  means  we  really  aren't  logging  anywhere. 

3.3  Installing  a  New  Routine 

In  order  for  a  routine  to  be  used  by  this  program  it  must 
be  written  in  a  particular  format.  This  format  is  outlined  in 
section  3.2.  PROCESS  CONTROL. 

Once  the  routine  is  properly  formatted  it  will  probably 
need  to  change  values  of  some  of  the  variables  included  in  common 
blocks.  The  following  is  a  list  of  the  files  in  which  the  common 
blocks  used  by  this  program  are  stored: 

MAGCMN.FOR  Contains  all  the  working  variables 

and  arrays  of  the  program. 


MAGCMN2 . FOR 


The  program  control  variables. 


HEADER. CMN  All  the  variables  that  define  the 

data  header  fields. 

FTRN.FOR  Common  variables  used  in  FTRAN  read 

routine. 


You  may  also  want  to  include  something  in  the  helps.  All 
the  helps  are  in  a  file  called: 


CIMAG2.HLP 


To  convert  this  into  a  library  file,  enter: 


LIBRARY/CREATE/HELP  CIMAG2.HLB  CIMAG2.HLP 
CIMAG2.HLB 


(Content  of  the  existing  HELP  file  is  listed  in  Appendix  A). 
There  is  another  utility  this  program  uses  called  the 
message  utility.  If  you  wish  to  use  this  utility  to  generate 
error  messages,  the  existing  message  file  is: 


CIMAGMSG.MSG 


When  the  utility  is  ran  the  output  will  be  put  into  a  file  named 


CIMAGMSG.OBJ 
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More  details  on  this  utility  are  included  in  the  Vax 


manuals. 


3.4  Listing  of  the  Linking  CoBaand  Piles 


$  ! 

$  !  This  command 
$  ! 

$  LINKXNOTRACE 
Cl MAG 2 
INTER, 

RDFLE, 

FORT, 

DDMPB, 
CIMAGMSG, 
'GRPllLIB' , 

' PL0T0LD2 ' 

$  ! 

$  EXIT 
$  ! 

$  !  This  command 
$  ! 

$  LINKXNOTRACE 
NEWMAG, 
INTER, 

RDFLE, 

FORT, 

DDMPB, 
CIMAGMSG, 
'GRPllLIB' , 
'PL0T0LD2' 


procedure  links  the  modules  for  the  CIMAG2  program. 


_!Main  program 

-•Frequency  Domain  Read 

-{Frequency  Domain  Read  (called  by  INTER) 

-•Fourier  Transform 

-*BSC  Read 

-{Program  Error  Messages 
-{Plotting 
{Plotting 


procedure  links  the  modules  for  the  NEWMAG  program. 


-{Frequency  Domain  Read 
-{Frequency  Domain  Read  (called  by  INTER) 
-{Fourier  Transform 
-{BSC  Read 

-{Program  Error  Messages 
-{Plotting 
{Potting 
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APPENDIX  A 


CIHAG2  HELP  LIBRART 


1  DCL  Commands 

The  user  may  use  DCL  commands  while  he  is  still  in  the 
program  by  simply  typing  '$'  before  the  command  he  wants  to 
use.  In  this  way  he  may  use  the  $DIR  command  to  see  the  files 
that  are  available  in  a  given  directory,  run  a  calibration 
program  and  then  use  this  new  data  with  the  data  that  he 
already  had  in  the  program,  etc... 

Warning: 

The  user  will  not  be  able  to  use  the 
$SET  DEFAULT  command  due  to  the  way  in  which 
the  DCL  commands  are  enabled.  The  default  will 
be  set  to  whatever  the  default  was  when  the 
user  entered  the  program.  This  means  he  will  to 
fully  specify  directory  names  if  he  wishes  to 
use  other  directories.  When  in  the  program 
however  he  may  use  the  DEF  command  to  set  the 
default  for  the  program's  read  statements. 

2  SSPAWN 

If  the  user  uses  the  $SPAWN  command  he  can  in  effect 
suspend  the  program  and  open  up  a  new  terminal.  From  here  he 
can  do  anything  he  wants.  When  he  is  through  he  simply  logs  out 
and  then  he  will  find  himself  back  in  the  program  in  the  same 
spot  he  left. 


1  BSC 

This  routine  reads  data  from  a  basic  scattering  code 
calculation. 

1  CHANGE 

This  command  allows  the  user  to  change  the  value  of  any 
harmonic  in  a  file.  This  was  originally  intended  for  use  in 
creating  test  files  but  the  user  may  find  other  uses  for  it. 

2  Parameters 

The  computer  will  ask: 

Which  harmonic  do  you  wish  to  change? 

(Answer  with  an  integer  between  1  and  4096  representing 
the  number  of  the  harmonic  you  wish  to  change.) 

Then  it  will  give  you  some  information  and  then  ask  you 
for  a  new  value: 


Current  Value:  [The  harmonic  number]  [The  current  value] 

New  Value: 

(Answer  with  a  DECIMAL  NUMBER.  This  is  most  important. 
If  the  user  forgets  to  put  in  the  decimal  point  then 
the  program  will  give  erroneous  results.) 

It  is  a  good  idea  to  use  the  TYP  command  to  check  the 
file  to  make  sure  that  the  results  turn  out  to  be  what  they 
were  intended  to  be. 

1  CREATE 

Creates  a  blank  time  domain  file  in  the  main  buffer. 

The  values  of  all  the  harmonics  will  be  zero. 

2  PARAMETERS 

The  program  will  ask  the  user  for  a  three  line  header. 
The  user  simply  types  in  each  line  of  the  header  finishing  each 
with  a  carriage  return.  When  all  three  lines  have  been  typed  in 
the  program  will  give  the  user  the  command  prompt. 


1  CLR 
1  CLR_BUF 
1  CMI 

image. 


Clears  the  main  buffer  by  setting  all  values  to  zero. 
Clears  all  buffers  including  the  main  buffer. 

Routine  to  set  up  the  parameters  for  a  two  dimensional 


1  COM 

The  combine  command  calculates  any  linear  combination 
of  the  data  in  any  number  of  the  storage  buffers  and  puts  the 
result  into  the  main  buffer. 

IMPORTANT  NOTE:  TIME  domain  waveforms  are  LINEAR,  so  combine 

is  ADDITION  or  SUBTRACTION.  FREQUENCY  spectra 
are  LOGARITHMIC,  so  combine  is  MULTIPLICATION 
or  DIVISION  of  spectra.  (Same  as  convolution 
or  deconvolution  in  the  time  domain.) 


2  Parameters 

It  starts  by  clearing  the  main  buffer.  Then  it  asks: 

STORAGE  BUFFER  NUMBER(USE  0  TO  FINISH): 

(reply  with  a  buffer  number  n  or  0  to  finish 
the  combination  process  and  CARRIAGE  RETURN) 
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MULTIPLIER  = 

(reply  with  a  floating  point  number  a  and 
CARRIAGE  RETURN). 

Then  the  main  buffer  will  be: 

BUF(0)=al*BUF(nl)+a2*BUF(n2)+a3*BUF(n3)  + . 

1  DATA 

This  command  puts  the  user  in  the  Database  program. 
This  is  a  separate  program  designed  to  find  data  files  using 
the  header.  It  has  it's  own  HELP  command  once  the  user  gets 
into  the  program.  When  the  user  wishes  to  return  to  the  CIMAG2 
program  he  simply  types  EXIT  and  he  will  find  himself  wherever 
he  left  off  before  he  entered  the  Database. 

1  DCV 

Routine  to  downconvert  a  spectrum  to  a  dual-polarity 
envelope.  The  user  chooses  the  harmonic  number  to  call  his 
"center  frequency"  which  gets  converted  to  the  DC  term. 

The  routine  works  on  a  frequency  spectrum  in  the  main 
buffer  and  returns  the  result  to  the  main  buffer. 

1  DEF 

Sets  the  default  directory  to  a  new  default  for  the 
program's  input  read  statements.  This  gives  the  user  the 
ability  to  input  data  files  from  anywhere  in  the  computer 
without  having  to  fully  specify  devices  and  directories  every 
time, 

2  Parameters 

It  will  ask: 

DEFAULT? 

The  user  should  answer  the  same  way  that  he  would 
for  the  DCL  command: 

Device: [directory] 
eg: 

USER1:[ STEVIE. 747] 

1  DEL 

Deletes  a  data  file  in  VAX  memory  device. 

2  Parameters 

ENTER  THE  FILE  NAME 

(reply  with  up  to  50  characters,  CARRIAGE 
RETURN) 


1  DIF 


The  differentiate  command  calculates  the  derivative  of 
the  TIME  DOMAIN  waveform  in  main  buffer  and  places  the  result 
in  main  buffer. 


1  DJW 

Divides  the  frequency  spectrum  in  main  buffer  by 
j*2*pi*f  (equivalent  to  integration  in  the  time  domain)  and 
places  the  result  back  in  main  buffer. 


1  ECHO 

For  use  with  a  command  file.  If  the  user  wishes  to  have 
the  commands  displayed  as  the  command  file  executes  them  he  may 
do  so  by  using  this  command.  To  turn  the  echo  off  again  he 
simply  uses  the  NO  ECHO  command. 


1  EXI 

Exit  Command 

1  FFT 

This  command  performs  the  "fast  Fourier  Transform"  on  a 
4096  point  time  domain  waveform  that  resides  in  the  main  buffer 
(buffer  0).  It  then  places  the  resulting  2048-harmonic  spectrum 
back  in  the  main  buffer. 

1  FGT 

This  frequency  spectrum  gating  routine  performs  a 
trapezoidal  gate  on  the  log  amplitude  components  (0$iC048) 
while  leaving  the  phase  spectral  components  (2049^<4096) 
unchanged.  It  modifies  the  spectrum  in  main  buffer. 

2  Parameters 

HARMONIC  FOR  START  OF  HIGH  FREQUENCY  CUTOFF 

(reply  with  integer  "a",  CARRIAGE  RETURN) 

HIGH  FREQUENCY  ROLL-OFF  IN  DB  PER  HARMONIC 

(reply  with  a  floating  point  value,  CARRIAGE 
RETURN) 

HARMONIC  FOR  START  OF  LOW  FREQUENCY  CUTOFF 

(reply  with  integer  "b",  CARRIAGE  RETURN) 

LOW  FREQ  ROLLOFF  IN  DB  PER  HARMONIC 

(reply  with  a  floating  point  value,  CARRIAGE 
RETURN) 

The  amplitude  spectrum  in  main  buffer  is  unchanged 
between  harmonics  a  and  b.  It  is  attenuated  at  the  rates 
specified  for  harmonics  above  and  below  those  values. 


33 


1  FILE 


This  conunand  gives  control  of  the  program  to  a 
specified  command  file.  When  the  new  file  takes  control  it  puts 
the  old  command  device  on  a  stack.  When  the  commard  file  is 
finished  it  may  then  return  control  to  the  old  device  or  it  may 
give  control  to  another  command  file  by  using  the  FILE  command 
itself.  However  it  must  eventually  pass  control  back  to  the 
device  that  gave  it  control.  Vhen  the  user  can  again  enter 
commands  he  will  be  given  the  prompt: 

CONTROL  HAS  RETURNED  TO  THE  TERMINAL 

2  Commandf iles 

The  command  files  are  just  what  the  name  implies:  a 
list  of  commands.  The  file  however  must  also  include  everything 
the  user  would  normally  type  in.  So  all  of  the  prompts  must  be 
answered  such  as: 

BUFFER#? 

The  easiest  way  to  do  this  is  through  the  use  of  the 
LOG  command. 

2  Parameters 

FILE?:  This  is  asking  for  the  name  of  the  new  command  file. 

1  FFREA 

Reads  a  data  file  written  by  [DOMI.DAT]FTRN.  (FTRN  REA 
Command) -Frequency  Domain. 


1  FRD 

Reads  a  frequency  domain  data  file  set  in  the  old 
standard  bands  of  1-2,2-4,4-8,8-12  GHz  into  the  main  buffer  in 
place  of  its  present  contents.  See  Bill  Leeper  for  further 
info. 


1  FTREA 


Reads  a  data  file  written 
Command) -Time  Domain 


by  [DOMI.DAT] FTRN. 


(FTRN  REA 


1  FTY 

File  type  control  for  the  data  processed  by  FTRAN. 

T=750  F=ll/23 

(Program  initially  sets  to  11/23) 

1  GAT 

This  command  performs  a  trapezoidal  gate  of  the  data  in 
the  main  buffer. 


2  Parameters 

START  GATE  OPEN (SAMPLE  NO.) 

(reply  with  integer  "a",  1^<4096) 

START  GATE  CLOSE  AT  SAMPLE 

(reply  with  integer  "b",  b>a, 1 <b<4096) 

RAMP  LENGTH  IN  SAMPLES  = 

(reply  with  integer  "c") 

Then  if  the  main  buffer  is  X(n),  and  the  result  to  be 
placed  in  the  main  buffer  is  X'(n): 

X' (n)=0  for  n<a 

X' (n)=X(n)*(n-a)/c  for  a<n<(a+c) 

X'(n)=X(n)  for  (a+c)^<b 
X' (n)=X(n)*(n-b)/c  for  b<n<(b+c) 

X'(n)=0  for  (b+c)^<4096 

1  GRID 

Creates  a  time  domain  file  that  has  a  pulse  at  every 
given  number  of  nanoseconds.  When  this  file  is  imaged  it 
effectively  gives  the  user  a  time  scale. 

WARNING:  Due  to  the  resolution  of  the  system  a  pulse 
will  usually  have  to  be  several  harmonics  wide,  depending  on 
the  widow,  in  order  for  it  to  be  picked  up.  It  is  a  good  idea 
to  use  the  TYP  command  to  see  how  many  pulses  you  should  see  in 
a  given  window  to  make  sure  that  they  were  all  picked  up  in  the 
scan. 

2  Parameters 

How  many  nanoseconds  per  division? 

(Answer  with  a  DECIMAL  number.  Make  sure  the  decimal  is 
Included. ) 

Type  in  header : (Three  lines) 

(Type  in  the  three  header  lines  ending  each  with  a 
CARRIAGE  RETURN.) 

How  many  harmonics  wide  should  the  line  be? 

(Answer  with  an  integer.  Remember  this  system  is  set  up 
for  continuous  wave  forms  and  not  pulses  so  the  pulse 
may  have  to  be  several  harmonics  wide;  i.e.,  I  have 
found  7  harmonics  to  work  best  for  a  window  size  of 
300.  As  the  window  gets  smaller  fewer  harmonics  will  be 
needed . ) 


This  command  performs  the  inverse  of  the  "fast  Fourier 
Transform"  on  a  2048-harmonic  spectrum  in  the  main  buffer  and 
places  the  resulting  4096-point  time  domain  waveform  back  in 
the  main  buffer. 

1  IMG 

Routine  to  form  the  image  for  one  polarization  from 
time  domain  waveforms  as  set  up  by  the  CMI  command. 

1  INT 

The  integrate  command  calculates  the  integral  of  the 
TIME  DOMAIN  waveform  in  main  buffer  and  places  the  result  in 
main  buffer. 

1  LOG 

Logs  the  user  input  into  a  command  file.  This  is  a  good 
way  to  build  command  files  for  the  FILE  command.  You  can  stop 
logging  with  the  ST0_L0G  command. 

2  Parameters 

LOG  FILE?:  Give  a  name  for  the  file  to  which  the  log  is 

to  be  sent. 

2  Command_files 

The  log  command  is  the  easiest  way  to  make  command 
files  for  the  FILE  command.  It  automatically  takes  all  the 
input  from  the  input  device  and  makes  a  file  out  of  it. 

If  the  user  wishes  to  put  comments  in  the  log  file  he 
may  do  so  by  using  ' ! '  as  the  first  character  of  a  command.  The 
program  will  ignore  it  but  it  will  be  logged  into  the  log  file. 
This  is  a  very  good  way  of  identifying  what  a  given  command 
file  does  or  even  what  a  section  of  one  does. 

The  only  other  thing  to  remember  when  making  command 
files  is  that  it  is  a  good  idea  to  set  the  default  in  the  very 
first  line  through  the  use  of  the  DBF  command.  In  this  way  a 
given  command  file  can  be  run  regardless  of  what  the  default 
was  set  at  previously. 


1  MANUALS 


There  are  two  manuals  for  the  CIMAG2  program: 


1  MJW 


1  MWR 


-  CIMAG2  USER'S  MANUAL  USER.TXT 

-  CIMAG2  PROGRAMMER'S  MANUAL  PROG.TXT 

The  user's  manual  gives  information  on  the  procedures 
used  to  create  images  and  manipulate  data.  The  programmer's 
manual  gives  information  that  might  be  helpful  to  someone 
wanting  to  make  changes  or  additions  to  the  program.  To  get  a 
copy  of  either  text  the  user  simply  needs  to  use  either  the  NEC 
command  to  get  a  copy  from  the  spinwriter  or  the  PRINT  command 
to  get  a  copy  from  the  printer.  This  can  be  done  from  inside 
the  program  or  from  the  monitor.  (See  DCL  Commands) 


$LASER  USER.TXT 


Multiplies  the  frequency  spectrum  in  main  buffer  by 
j*2*pi*f  (equivalent  to  differentiation  in  the  time  domain)  and 
places  the  result  back  in  main  buffer. 


Multiplies  the  frequency  domain  spectrum  in  main  buffer 
by  2*pl*f  and  places  the  result  back  in  the  main  buffer. 


1  NO  ECHO 


Turns  off  the  echo  that  was  enabled  by  the  ECHO 


1  NOR 


1  PIM 


1  PLO 


command . 


This  normalize  command  calculates  the  mean  of  the  4096 
point  waveform  in  the  main  buffer  and  shifts  the  main  buffer 
waveform  so  its  mean  is  zero. 


This  routine  plots  an  isometric  view  of  a  single 
polarization  two-dimensional  target  image  on  the  plot  device  In 
Isometric  form  with  no  shadowing. 


The  PLOT  command  allows  the  user  to  plot  any  one  of  the 
buffers  on  one  of  the  plotting  devices  listed. 


2  Parameters 


j 


•r-'i 


WAVEFORM  BUFFER  NUMBER? 

(reply  with  the  number  of  the  buffer  in  which 
the  desired  waveform  is  stored) 

DO  YOU  WANT  A  NEW  WINDOW?  (1=Y,0=N) 

(this  allows  you  to  choose  what  section  of  the 
wave  you  want  to  look  at  or  you  may  look  at  the 
whole  thing.  If  your  answer  is  1  then  it  will 
ask  the  for  the  range  of  data  numbers  you  want 
to  look  at  (1^<4096).  If  your  answer  is  0  then 
it  will  default  to  the  window  that  you  used 
last . ) 

DO  YOU  WANT  NEW  AXES? 

(this  allows  you  to  fix  the  labeling  for  the 
graph) 

INPUT  TITLE  FOR  PLOT? 

(this  is  the  title  for  the  top  of  the  plot. 

Type  in  any  title  desired.) 

Then  it  gives  you  a  list  of  devices  on  which  you  may 
output  the  plot.  Just  type  the  number  of  the  device 
which  you  wish  to  use. 


1  PROC 


This  operation  allows  the  user  to  specify  a  process  to 
be  performed  on  a  group  of  files.  Then  it  will  ask  for  the  list 
of  data  files  and  a  list  of  the  output  filenames  desired.  It 
will  then  run  each  of  the  data  files  through  the  defined 
process  and  make  output  files  as  requested. 

2  Procedure 

The  user  is  given  explicit  instructions  throughout  the 
procedure.  There  are  only  a  few  points  that  need  to  be  stressed 
here: 

-  The  filename  NAME.DAT  will  be  used  for  all 

input  data  files.  This  has  to  be  typed 
in  capital  letters. 

-  The  filename  NEWNAME.DAT  is  used  for  all 

output  file  names.  It  also  has  to  be 
typed  in  capital  letters. 


-  When  making  either  of  the  filename  lists  if  a 

given  file  is  used  mere  than  once  it 
has  to  be  listed  more  than  once. 

-  Filenames  are  used  in  the  order  that  they 

appear  in  the  lists. 

-  Whenever  the  user  finishes  defining  a  process 

or  a  list  he/she  will  type  the  word 
"DONE"  in  capital  letters  and  the 
procedure  will  move  on  to  the  next 
phase. 

WARNING:  This  is  a  very  powerful  command  but  the  user  must  be 
very  careful  when  using  it  so  that  the  output  files  will 
have  the  desired  content  and  not  some  other  content  without 
the  users  knowledge. 


1  PSM 


Point  smooth  command  fits  a  cubic  curve  to  points 
surrounding  a  small  bad  region  of  time  waveform  or  spectrum  for 
the  data  in  main  buffer. 

2  Parameters 

For  TIME  waveforms  it  asks: 

FIRST  BAD  POINT  INTEGER* 

(reply  with  the  integer  for  the  first  point  to 
be  replaced,  CARRIAGE  RETURN) 

LAST  BAD  POINT  INTEGER* 

(reply  with  the  integer  for  the  last  point  to 
be  replaced. This  can  be  the  same  as  the  first 
point) 

For  FREQUENCY  spectrum  in  main  buffer, it  asks: 

FIRST  BAD  HARMONIC 

(reply  with  a  harmonic  #  between  0  and 
2048,  CARRIAGE  RETURN) 

LAST  BAD  HARMONIC 

(reply  with  a  second  harmonic  #,  perhaps  same 
as  first) 

For  both  cases, the  bad  section  of  data  is  replaced  in  main 
buffer  with  values  which  join  surrounding  values  with 
continuous  slope. 


1  RBF 


This  command  reads  the  contents  of  a  storage  buffer 
into  the  main  buffer. 


2  Parameters 

BUFFER#:  There  are  35  buffers.  This  number  is  used  to  specify 
which  one  is  to  be  used. 

1  RDFL 

Reads  a  modern  frequency  domain  file  into  the  main 
buffer  in  place  of  its  present  contents.  It  asks: 

INPUT 

1  REA 

Reads  a  time  domain  file  in  "Jon  Young  format"  from  a 
VAX  storage  device  into  the  main  buffer  (destroys  what  was  in 
main  buffer  before). 

2  Parameters 

INPUT  FILE  ? 

(reply  with  up  to  50  characters  for  a  VAX  file 
name) 

WAVEFORM  NUMBER? 

(reply  with  an  integer, since  "Jon  Young  format" 
files  will  accept  more  than  one  waveform  per 
file;  normally  1  is  used  with  those  files  which 
have  only  one  waveform  in  them) 

If  the  name  is  unrecognized  or  if  any  other  error  happens,  the 
routine  aborts  and  the  "ERROR  IN  COMMAND"  message  is  printed. 

1  RLB 

The  re-label  command.  It  will  print  out  the  first  line 
of  the  old  title  block  and  then  wait  for  the  user  to  type  in 
the  new  line. 

OLD  TITLE  BLOCK: 

Whatever  the  line  is 

NEW  TITLE  BLOCK: 

The  user  then  simply  types  in  whatever  the  new  line  is  to  be. 

1  ROT 

Rotates  the  data  in  main  buffer  by  "a"  increments  to 
the  right(+a)  or  lef t(-a) ,with  values  shifted  beyond  0  or  4096 
appearing  on  the  other  end  of  the  waveform.  It  asks: 

ROTATE  BY  INCREMENTS  OF: 

(reply  with  "a",  CARRIAGE  RETURN) 


This  command  stores  the  contents  of  the  main  buffer 
into  one  of  the  storage  buffers. 


1  SBF 


2  Parameters 

BUFFER#:  There  are  35  buffers.  This  number  is  used  to  specify 
which  one  is  to  be  used. 

1  SCT 

This  smooth  cutoff  tail  routine  attaches  a  shifted 
cosine  amplitude  spectrum  attenuation  function  to  the  spectrum 
in  main  buffer  to  eliminate  Gibbs  phenomenon  in  the  time 
waveform. 

2  Parameters 

INPUT  THE  STARTING  AND  END  POINT  FOR  THE  FREQ  ROLLOFF 
(reply  with  a  pair  of  integers, 

"a",  "b",  0^,b,<4096) 

The  spectrum  is  unchanged  for  n^.  It  is  -100  db  for  n>b.  And 
for  a^<b,the  spectrum  has  a  cosine  rolloff, 

1  SH0_BUF 

This  command  shows  the  contents  of  all  of  the  buffers 
that  have  anything  in  them.  The  buffer  number  is  given  for  each 
entry.  The  domain  of  the  file  is  also  given  (time  or  freq). 

This  is  the  format: 

BF#  DOMAIN 
DESCRIPTION  . . . 

1  ST0_L0G 

The  stop  logging  command  stops  the  input  from  going  to 
the  logging  file.  It  assumes  that  the  user  was  writing  to  a 
logging  device. 

1  TYP 

Type  a  portion  of  data  in  a  buffer.  This  data  is  typed 
out  in  numerical  form  so  that  the  user  can  see  the  actual  data 
values  of  Individual  points. 

2  Parameters 

START  AT  ELEMENT  NUMBER: 

(reply  with  an  integer  from  1  to  4096,  CARRIAGE 
RETURN) 

END  AT  ELEMENT  NUMBER: 

(reply  with  a  second  integer,  1  to 
4096,  CARRIAGE  RETURN) 

BUFFER  NUMBER(0=MAIN): 

(reply  with  a  buffer  integer  #,0  to  8,  CARRIAGE 
RETURN) 

The  specified  portion  of  the  specified  buffer  is  then  typed  out 
on  tne  CRT  screen. 


This  routine  attaches  a  smooth  cutoff  tail  to  the 
spectrum  in  order  to  reduce  Gibb's  Phenomenon. 

2  Parameters 

INPUT  HARMONICS; START, END, TYPE  OP  WINDOW 

TYPE : 0-HANNING , 1 =HAMMING , 2-GAUSSI AN ,  0 , N , M-TEST 
(reply  with  a  pair  of  integers, 

"a",  '’b'',0^,b,<4096  and  then  0,  1,  or  2 
depending  on  the  type  of  rolloff  desired) 

1  WRI 

Writes  a  time  domain  waveform  in  the  main  buffer  to  a 
specified  VAX  memory  device  in  "Jon  Young  format".  The  data  in 
main  buffer  is  not  changed. 

2  Parameters 

FILE  NAME  ? 

(reply  with  a  name  up  to  50  characters, 
CARRIAGE  RETURN) 

1  YSH 

Shifts  the  data  in  buffer  up(+)  or  dovn(-)  by  Y  units. 

2  Parameters 

YSHIFT  IN  UNITS(-2048<y<2048)= 

(reply  with  the  shift  value,  CARRIAGE  RETURN) 


non  nn  nn  n  nn  n  n  onnn 


APPENDIX  B 


THE  COLOR  IMAGING  PROGRAM  'CIMAG2' 


CIMAG2 


CHARACTER»80  COMMAND  !  the  command 

INCLUDE  'USER2 : ( DURAL . Cl MAG2 ) FTRN . FOR ' 

INCLUDE  ’USER2 : I  DURAL . C I MAG2 1 MAGCMN . FOR ' 

INCLUDE  ’USER2.- (DURAL. CIMAG2iMAGCMN2. FOR'  !  common  block;. 

INCLUDE  ' USER2 : [ DURAL . Cl MAG2  j HEADER . CMN ' 


EQUIVALENCE  ( HEADER, HEAD ( 1 , 1 ) ) , ( LINEl ( 1 ) , HEADER (1 ) ) 

EQUIVALENCE  ( FLTP, IDS( 1 ) ) , (HEADER! 61 ) ,LINE3( 1 ) ) , (HEADER! 3 1 ) , L I NE2 ( 1  t  I 


INTEGER*4  FOR_RETCODE 
INTEGER*4  FOR_EOF  /  -1  / 
INTEGER*4  RETCODE 
INTEGER*4  LIBSGET  LUN 


1  fortran  return  code 
!  end  of  file  code 
1  return  code  for  RTL  loutine^ 
1  'get  LU# '  RTL  routine 


CALL  START 


initialize  all  the  vaiiables 


RETCODE  -  LIB$GET_LUN(  TERM_UNIT  ) 

OPEN!  FILE-'lf!'  ,  1  open  terminal  foi  J  (' 

+  UNIT-TERM_UNIT  , 

+  STATUS- 'NEW'  ) 


RETCODE  -  l:B$GET_LUN!  NULL_UN1T  )  !  null  unit  for  logging 

OPEN!  FILE-'NL:'  , 

♦  UNIT-NULL  UNIT  , 

+  STATUS- 'NEW'  ) 


COM_UN I T-TERM_UNI T 
IOUT-TERM_UNIT 
LOG  UNIT-NULL  UNIT 


1  get  commands  from  terminal  initially 
I  send  output  to  the  terminal 

I  log  to  null  unit 


C  from 


t 


DO  WHILE  (COMMAND  .NE.  'EXI') 

IF  (COM  UNIT  .EQ.  TERM  UNIT)  THEN  1 


erminal 

lOUT-TERM  UNIT 

I 

then  all 

WRITE!  UNIT-IOUT, 

FHT-5 

) 

ELSE 

IF  (ECHO)  THEN 

I OUT-TERM 

UNIT 

I 

1 

ELSE 

(or  not 

I  do  until  exit  command 
if  waitng  for  a  command 

prompts  go  to  terminal 

if  using  a  command  filo 
output  can  either  be 
(displayed  on  terminal 
be  displayed  at  all 


lOUT-NULLUNIT 

END  IF 


C 

C  control 


END  IF 

READ  (  UNIT-COMUNIT  ,  FMT-10  , 

IOSTAT-FOR_RETCODE  )  COMMAND 

IF  (FORRETCODE  .EQ.  FOREOF)  THEN 

device 


!  read  the  command 
on  end  of  file  get  old 


CALL  BACK 

IF! PROCFLAG)  THEN  I  if  coming  back  from  a  pioceduip 
CALL  PROC3  1  finish  off  the  PROC  piocrdut. 


END  IF 


CALL  STR$UPCASE(  COMMAND  ,  COMMAND  )l  convert  to  uppercase 


i 


C  not  equl  to 
+ 

C  run  conaand  file 


C  the  input 


ir  ((COMMAND  .NE.  'STO_LOG')  .AND. 

(COMMAND  .NE.  'FILE')  .AND. 


(  if  command  is 
t  stop  log  or 


(COMMAND  .NE.  'DONE'))  THEN 

WRITE(  UNIT-LOG_UNIT  ,  FMT-IO  )  1  log 

COMMAND 


IF  (COMMAND  .EQ.  'BSC')  THEN  I  execute  the  command 
CALL  BSC 

C  ELSE  IF  (COMMAND  .EQ.  'CARD')  THEN 

C  CALL  CARDC 

ELSE  IF  (COMMAND  .EQ.  'CHANGE')  THEN 
CALL  CHANGE 

ELSE  IF  (COMMAND  .EQ.  'CLR')  THEN 
CALL  CLR 

ELSE  IF  (COMMAND  .EQ.  'CLRBUF')  THEN 

CALL  CLR_BUF(HEAD,  LINEl ,  LINE2,  LINE3) 

ELSE  IF  (COMMAND  .EQ.  'CMI')  THEN 
CALL  CMI 

ELSE  IF  (COMMAND  .EQ.  'COM')  THEN 

CALL  COM(HBAD,  LINEl,  LINE2 ,  LINE3 ) 

ELSE  IF  (COMMAND  .EQ.  'CREATE')  THEN 

CALL  CREATE(HEAD,  LINEl,  LINE2,  LINE3) 

ELSE  IF  (COMMAND  .EQ.  'DATA')  THEN 

CALL  LIB$SPAHN(  'RUN  USER2 : ' 

*  //' (DURALIDBMANAGER' )  1  use  the  database 

ELSE  IF  (COMMAND  .EQ.  'DCV')  THEN 
CALL  DCV 

ELSE  IF  (COMMAND  .EQ.  'DEF')  THEN 

CALL  DEFFER  I  set  def  command 

ELSE  IF  (COMMAND  .EQ.  'DEL')  THEN 
CALL  DEL 

ELSE  IF  (COMMAND  .EQ.  'DIF')  THEN 
CALL  DIF 

ELSE  IF  (COMMAND  .EQ.  'DJH')  THEN 
CALL  DJW 

ELSE  IF  (COMMAND  .EQ.  'DONE')  THEN 
CALL  PROC2 

ELSE  IF  (COMMAND  .EQ.  'DVJR')  THEN 
CALL  DWR 

ELSE  IF  (COMMAND  .EQ.  'ECHO')  THEN  1  shows  command 

file  i  n«:  t- met  1  ons  as  they  are 

ECHO-. TRUE.  1  being  executed. 

ELSE  IF  (COMMAND  .EQ.  'FFT')  THEN 
CALL  FFT 

ELSE  IF  (COMMAND  .EQ.  'FGT')  THEN 
CALL  FGT 

ELSE  IF  (COMMAND  .EQ.  'FILE')  THEN 
CALL  FILE 

ELSE  IF  (COMMAND  .EQ.  'FRD')  THEN 
CALL  FRDC 

ELSE  IF  ( COMMAND. EQ. 'FTREA' )  THEN 
CALL  FTREA 

ELSE  IF  ( COMMAND. EQ. 'FFREA' )  THEN 
CALL  FFREA 

ELSE  IF  ( COMMAND. EQ. 'FTY' )  THEN 
CALL  FTY 

ELSE  IF  (COMMAND  .EQ.  'GAT')  THEN 


C  help  files 


C  DCL  command 


CALL  GAT 

ELSE  IP  (COMMAND  . EQ .  'GRID')  THEN 

CALL  GRID(HEAD,  LINEl,  LINE2,  LINE3) 

ELSE  IF  (COMMAND  .EQ.  'HELP')  THEN 

CALL  LIB$SPAWN( 'HELP/PAGE/LIBRARV-USER2; ' 

//' (DURAL. CIMAG2 )CIMAG2 .HLB' )  !  ca) 

ELSE  IP  (COMMAND  .EQ.  'IFF')  THEN 
CALL  IFF 

ELSE  IF  (COMMAND  .EQ.  'IMG')  THEN 
CALL  IMG 

ELSE  IF  (COMMAND  .EQ.  'INT')  THEN 
CALL  INTEG 

ELSE  IF  (COMMAND  .EQ.  'LOG')  THEN 
CALL  LOGGER 

ELSE  IF  (COMMAND  .EQ.  'MJW')  THEN 
CALL  MJW 

ELSE  IF  (COMMAND  .EQ.  'HHR')  THEN 
CALL  HHR 

ELSE  IF  (COMMAND  .EQ.  'NOECHO')  THEN 
ECHO-. FALSE. 

ELSE  IF  (COMMAND  .EQ.  'NOR')  THEN 
CALL  NOR 

ELSE  IF  (COMMAND  .EQ.  'PIM')  THEN 
CALL  PIN 

ELSE  IF  (COMMAND  .EQ.  'PLO')  THEN 
CALL  PLO 

ELSE  IF  (COMMAND  .EQ.  'PROC')  THEN  1  process  data 
CALL  PROC 

ELSE  IF  (COMMAND  .EQ.  'PSH')  THEN 
CALL  PSM 

ELSE  IF  (COMMAND  .EQ.  'RDFL')  THEN 

CALL  BDFLiHEAD,  LINEl,  LINE2,  LINES) 

ELSE  IF  (COMMAND  .EQ.  'REA')  THEN 

CALL  REA(HEAD,  LINEl,  LINE2,  L1NE3) 

ELSE  IF  (COMMAND  .EQ.  'RBF')  THEN 

CALL  RBF(HEAD,  LINEl,  LINE2,  LINES) 

ELSE  IF  (COMMAND(l:2)  .EQ.  'RB')  THEN 

NB  -  ICHAR(  COMMANDO:  3)  )  -  48 
CALL  RB(HEAD,  LINEl,  LINE2,  LINE3) 

ELSE  IF  (COMMAND  .EQ.  'RLB')  THEN 
CALL  RLB 

ELSE  IF  (COMMAND  .EQ.  'ROT')  THEN 
CALL  ROT 

ELSE  IF  (COMMAND  .EQ.  'SBF')  THEN 

CALL  SBF(HEAD,  LINEl,  LINE2 ,  LINES) 

ELSE  IF  (COHNAND(l:2)  .EQ.  'SB')  THEN 

NB  -  ICHAR(  COMMANDO:  3)  )  -  48 
CALL  SB(HCAD,  LINEl,  LINE2,  LINES) 

ELSE  IF  (COMMAND  .EQ.  'SHOBUF')  THEN 
CALL  SHOBUF 

ELSE  IF  (COMMAND  .EQ.  'STO_LOG')  THEN 
CALL  STOLOG 

ELSE  IF  (COMMAND  .EQ.  'TVP')  THEN 

CALL  TyP(HEAD,  LINEl,  LINE2,  LINES) 

ELSE  IF  (COMMAND  .EQ.  'WND')  THEN 
CALL  WND 

ELSE  IF  (COMMAND  .EQ.  'WRI')  THEN 

CALL  HRK HEAD,  LINEl,  LINE2,  LINES) 

ELSE  IF  (COMMAND  .EQ.  'YSH')  THEN 
CALL  YSH 

ELSE  IF  (COMMAND  .EQ.  '  1')  THEN 

CONTINUE  !  clause  on  the  old  REA  command 

ELSE  IF  (COHMAND(l:l)  .EQ.  '$')  THEN  !  If 

CALL  LIB$SPAWN(  COMMAND  )  !  execute  the  DCL 


T  o  n  o  n  n  n  n  n  noo  oi-tnoo 


C  conmand 

ELt,E  IP  (C0HHAND(1:1)  .EQ.  '1')  THEN  I  if 

C  command  file  comment 

IP  (COM  UNIT  .NE.  TERMUNIT)  THEN  !  if 

C  commands  are  not  coming  from  the  terminal 

WRITEdOUT.lO)  COMMAND  I  output  the 

C  comment 

END  IP 

ELSE  IP  (COMMAND  .NE.  'EXI')  THEN  1  If  not  exit 

C  then  it  must  be  undefined 

WRITE( lOUT,*)  '7' 

END  IP 

END  IP 

C 

END  DO 


rORMAT(  lx  ,  'O'  ,  S  ) 
PORMAT(  A  ) 

STOP 

END 


SUBROUTINE  COM  FILE (HEAD,  LlNEl,  LINE2,  LINB3) 


INCLUDE  'USER2: (DURAL. CIHAG2]PTRN. FOR' 
INCLUDE  ' USER2 : ( DURAL . CINAG2 ] NAGCMN . FOR ' 
INCLUDE  'USER2! [ DURAL . CINAG2 ] MAGCNN2 . FOR' 
INCLUDE  'USER2: (DURAL. CIHAG2jHEADER. CNN' 
INCLUDE  'USER2: (DURAL. CIMAG2jMSGBLK. FOR' 


INTEGER*4  SYS$ALLOC  , 
■f  SYS$DALLOC 


!  do  a  fortran  stop 


!  the  er 


1  allocate  routine 

1  deallocate 


INTEGER*4  FOR  RETCODE 
INTEGER*4  RETCODE 

integer*4  ss$  normal  , 

SS$~DEVALLOC 


I  fortran  return  code 
t  return  code  for  system  operations 
I  normal  return 

i  device  already  allocated  to  another  process 


PARAMETER  ( SS$  NORMAL  -  ' 00000001 'X) 
PARAMETER  ( SSS'DEVALLOC  -  ' 00000840 'X) 


THE  INITIALIZATION  ROUTINE 


'J  rxTFnnAf.  CTPL  C  HDLR 

C  CALL  LIB$ESTAgLTSH(CTRL_C_HDLR)  !  Ini'iTi.  .'• 

ENTRY  START 
FTYPE-P 
ECHO-. FALSE. 

PROCFLAG-. FALSE. 

DEFINE  FLAG-. FALSE. 

ISIZE-7096 

IS-ISIZE/2. 

FLTP-0. 

IHDSZ-256 

IDSSZ-1 

PI-3.14159 

IZ-12 

ISM2-( ISIZE/2)-2 
NBUF-30 


IMEG-0 
IWINl-1 
WINl-1 . 

IWIN2-2048 

WIN2-2048. 

WIN3— 1024  . 

WIN4-1024 . 

YHAX-1 
YMIN--1 
YTMI-1 
YLAB-'  ' 

XMAX-2048 
XMIN-1 
XTMI-512 
XLAB- '  ' 

EFLAG  -  .FALSE. 

data  CARRAY/0 ,0,0, 0,0, 0,6, 0,6, 0,7, 0,6, 0,6, 1,6, 1,6, 2, 5, 
2,5,3,4,6,4,6,1,11,1,11/ 

DATA  CLRTAB/6 , 3 , 7 , 2 , 4 , 5/ 


NOW  CONSTRUCT  FFT  SXNE  TABLE 
CALL  FORT ( A, I Z,S, 0,1 TERR) 

TYPE  Instructions  on  CIMAG2  available  through  the  command  HELP' 


RETURN 


1  Initialization  is  finished 


ENTRY  FFT 
DOMAIN- 'FREQ' 

FLTP-1  . 

DO  20  N-1,ISIZE 
A(N)-CMPLX( 1 . ,0. )*P(N,1) 

CALL  F0RT(A,IZ,S,-2,IFERR) 

1F( IFERR.NE.O)  THEN 

CALL  LIBSSIGNAL(MAG_COH) 

RETURN 

END  IF 
RT2-SQRT( 2 . ) 

X-CABS( A( 1 ) ) 

IF(X.LT.1.E-20)  X-1.E-20 

P(  1 , 1  )-10.24»(  10.-*-20.»ALOG10(X)  ) 

IS-ISIZE/2 . 

P( IS+1 ,1 )-512. *ATAN2(AIMAG(A( 1 ) ) ,REAL(A( 1 ) ) )/PI 

PMAX--1.E-20 

no  71  N-2,IS 

nn-is+n 

X-CABS(A(N) ) 

IF(X.LT.1.E-20)  X-l.E-20 

P(N,1 )-10.24*(10.+20.*ALOG10(X*RT2) ) 

1F(P(N,1) .GT.PMAX)  PHAX-P(N,1) 

X-CABS(A(N) ) 

IF(X.EO.0)P(NN,  0 

IF(X.NE.0)P(NN,i)-512.*ATAN2(AIMAG(A(N) ) ,REAL(A(N) ) )/PI 
RETURN 


I FFT  ROUTINE 
ENTRY  IFF 
DOMAIN-'TIME’ 
FLTP-O  . 

RRT2-1 ./S0RT( 2 . ) 


-  -  j.  •  -  •  w , 


m 


P(1,1)-(P{1,1)/10.24)-10. 

PUS+l!n-PI*P(  IS+1,1)/512. 

A(1)-CMPLX(C0S(P(IS+1,1) ) ,SIN(P(1S+1,1) ) )*10.**(P(1,1)*0.05) 

DO  30  N-2,IS 

P(N,1)-(P(N,1)/10.24)-10. 

P{  IS+N,  1  )-PI*P(  IS-t-N,!  )/512. 

A(N)-CMPLX{COS(P( IS+N,1 ) ) , SIN( P( IS+N, 1 ) ) ) *RRT2 *1 0 . • * ( P ( N , 1 ) * . 0 5 ) 
A(IS+1)>CHPLX(0. ,0.  ) 

DO  31  N-2,IS 
NN-ISIZE+2-N 

A(NN)-CMPLX(REAL(A(N) ) ,-AIHAG(A(N) ) ) 

CALL  rORT(A,IZ,S,2,irERR) 

IF( IFERR.NE. 0)  THEN  1  Print  an  error  message 

CALL  LIB$SIGNAL(HAG_COH) 

RETURN 

END  IF 

DO  32  N-1,ISIZE 

NN-ISIZE+l-N 

P(NN, 1 )-REAL(A(NN) ) 

RETURN 

ENTRY  FTY 

WRITE( lOUT, *) 'Enter  file  type  for  FTRAN  read  routines' 

WRITE( lOUT, * ) '  T-750  F-11/23' 

READ(CON_UNIT,* )FTYPE 
WR1TE( LOG_UNIT, * ) FTYPE 
RETURN 

FTRAN  'REA 'COMMAND  (FREQUENCY  DOMAIN) 

ENTRY  FTREA 
FLTP-1 

I r(  .NOT. FTYPE)  CALL  FREA ( INPFI LE , AMPL , PHS ) 

IF( FTYPE)  CALL  REU( INPFILE , AMPL , PHS ) 

RETURN 

FTRAN  'REA'  COMMAND  (TIME  DOMAIN) 

ENTRY  FFREA 
FLTP-O 

IF( .NOT. FTYPE)  CALL  FREA( INPFI LE , AMPL , PHS ) 

IF( FTYPE)  CALL  REU ( INPFILE .AMPL , PHS ) 

RETURN 


DIFFERENTIATE  ROUTINE 
MULT  BY  JW  IN  THE  FREQ  DOMAIN 
ENTRY  MJW 
SN-1  . 

GOTO  40 

INTEGRATE  ROUTINE 

DIVIDE  BY  JW  IN  THE  FREQ  DOMAIN 

ENTRY  DJW 

SN— 1 . 


IS-ISIZE/2. 

IF  (FLTP.NE.l. )  THEN 

CALL  LIB$SIGNAL(MAG_CON) 

END  IF 


41 


C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

45 


46 


c 

c 

c 

c 

c 

c 


DO  41  N-2,IS 
NN-IS+N 

P(NN, 1 )-P(NN, 1 )+SN*256 . 

P(N,1)-P(N,1)+10.24*SN*20.*ALOG10(FLOAT(N-1) ) 
P(l,l)— 100 
P( IS  +  1 , 1  )-0  . 

RETURN 


FILTER  ROUTINE 

MULT  BY  W  IN  THE  FREQ  DOMAIN 

ENTRY  MWR 
SN-1  . 

GOTO  45 

FILTER  ROUTINE 

DIVIDE  BY  W  IN  THE  FREQ  DOMAIN 

ENTRY  DWR 
SN--1 . 


IS-ISIZE/2 . 

IF  ( FLTP .NE. 1 . )  THEN 

CALL  LIB$SIGNAL(MAG_COM) 

END  IF 

DO  46  N-2,IS 

P(N, 1 )-P(N, 1 )+10.24*SN*20. *ALOG10( FLOAT(N-l ) ) 
P(  1  ,  1  )— 100 
P( IS+l,l)-0. 

RETURN 


GATE  COMMAND 

TRAPEZOIDAL  GATE  OF  TIME  DOMAIN  WAVEFORM 
ENTRY  GAT 


48 

FORMAT( FIO . 2 ) 

I 

49 

FORMAT  (15) 

WRITE  (IOUT,50) 

50 

FORMAT!’  START  GATE  OPEN (SAMPLE  NO.)-' 
READ  (COM  UNIT, 49)  lOP 

WRITE! LOG  UNIT, 49)  lOP 

S‘ 

WRITE! lOUT, 51 ) 

51 

FORMAT!'  START  GATE  CLOSE  AT  SAMPLE  ') 

READ  (COM  UNIT, 49)  ICL 

WBTTF  (rnr:  UNIT, 49)  ICL 

WRIIE! lOUl , 52 ( 

52 

FORMAT!'  RAMP  LENGTH  IN  SAMPLES  -') 

'  A 

READ  (COM  UNIT, 49)  IRMP 

WRITE  (LOG  UNIT, 49)  IRMP 

IF( IRMP.EQ.0)IRMP-1 

% 

DO  53  N-1,ISIZE 

F-0  . 

RMP-FLOAT( IRMP) 

NUM-N-IOP 

FNUM-FLOAT(NUM) 

IF(N.GT. IOP)F-FNUM/RMP 

IF(F.GT.l. )F-1. 

NUM-N-ICL 

FNUM-FLOAT(NUM) 

IF(N.GT. ICL)F-1 .-FNUM/RMP 

IF( F.LT.O. )F-0. 

!  log  the  input 


!  log  the  input 


P(N,1 )-P(N,l )*F 
RETURN 


TYPE  COMMAND 

TYPES  A  BUFFER  WAVEFORM  SEGMENT 
ENTRY  TYP(HEAD,  LINEl,  LINE2,  LINE3) 

WRITE( IOUT,81 ) 

FORMAT ( '  START  AT  ELEMENT  NUMBER:',?) 

READ  (COM_UNIT,86)  NP 
WRITE  { LOG_UNIT, 86 )  NP 
WRITE( lOUT, 82 ) 

FORMAT ( '  END  AT  ELEMENT  NUMBER:',?) 

READ  (COM_UNIT,86)  NE 
WRITE  (LOG_UNIT,86)  NE 
IF(NE.GE. ISIZE)NE-ISIZE-1 
WRITE! lOUT, 83 ) 

FORMAT!'  BUFFER  NUMBER 1  0-MAIN ):',? ) 

READ  !COM_UNIT,86)  NNB 
WRITE  !LOG_UNIT,86)  NNB 
NNB-NNB+1 
IF!NNB.LT.l)  THEN 

CALL  LIB?SIGNAL!MAG_COM) 

RETURN 

END  IF 

IF!NNB.GT.NBUF)  THEN 

CALL  LIB?SIGNAL!MAG_COH) 

RETURN 

END  IF 

WRITE! lOUT, 117 )  !HEAD( I ,NNB) , t-1, 30) 

WRITE! lOUT, 117 )  IHEAD! I , NNB ) , 1-31 , 60 ) 

WRITE! IOUT,117)  ! HEAD! I, NNB), 1-61, 90) 

WRITE! lOUT, 85) !N,P! N, NNB) , P(N+1 ,NNB) ,N-NP,NE, 2 ) 
FORMAT!I5,2F10.2) 

FORMAT  !I5) 

RETURN 


BUFFER  STORE  ROUTINES 

STORE  A  WAVE  OR  SPECTRUM  IN  ONE  OF  30  TEMP  LOGS 
ENTRY  SBF!HEAD,  LINEl,  LINE2,  L1NE3) 

WRITE  !IOUT,87) 

FORMAT! IX,'  BUFFER  #',?) 

READ  1C0M_UNIT,93)  NB 
WRITE  !LOG_UNIT,93)  NB 

ENTRY  SBlHEAD,  LINEl,  LINE2,  LINE3) 

IF  !Nn  .GT.  35)  THEN 

CALL  LIB?S1GNAL|MAG_C0H) 

RETURN 

END  IF 

DON! NB) -DOMAIN 
DO  90  I-  1,  30 
BUFFERS !NB, I )  -  LINEl! I) 

CALL  BUFSTRI IHDSZ , IDSSZ , ISIZE ,NB , HEAD, P, IDS ) 
FORMAT  ! I  5 ) 

RETURN 


BUFFER  READ  ROUTINES 

RETRIEVES  WAVE  OR  SPECTRUM  AND  ITS  HEADING  AND  ITS 
DESCRIPTIVE  PARAMETERS  INTO  THE  MAIN  BUFFER 
ENTRY  RBF!HEAD,  LINEl,  LINE2,  LINE3) 

WRITE  !IOUT,88) 


1  log  the  input 


1  log  the  input 


log  the  input 


1  Print  an  error  messag 


1  Print  an  error  messag 


I  log  the  input 
[  SB*  routine 


88 


FORMATdX,'  BUFFER  #’,S) 
READ  (COM_UNIT,94 )  NB 
WRITE  (LOG  UNIT, 94)  NB 


1  log  the  input 


1 

c 

■ 

92 

94 

c 

C 

C 

C 

C 

C 

C 

X* 

C 

101 

? 

102 

w 

103 

•j'" 

It 

104 

105 

I 

107 

108 

110 

111 

> 


■J; 

r. 


112 


113 

114 
C 


115 


ENTRY  RB(HEAD,  LINEl ,  LINE2,  LINE3) 

DOMAIN  -  DOM(NB) 

CALL  BUFRD( IHDSZ , I DSSZ , I S I ZE , NB , HEAD, P , I DS ) 
FORMAT  (15) 

RETURN 


COMBINE  ROUTINE 

COMBINES  FILES  FROM  TEMP  BUFFERS, EACH  MULTIPLIED  BY  A 
FACTOR, INTO  THE  MAIN  BUFFER.  HEADING  AND  DESCRIPTIVE  PAR 
OF  THE  LAST  WAVEFORM  WILL  BE  IN  THE  MAIN  BUFFER 

ENTRY  COM(HEAD,  LINEl,  LINE2,  LINE3) 

DO  101  I-1,ISIZE 
P(I,l)-0 
WRITE( lOUT, 103 ) 

FORMAT ( '  STORAGE  BUFFER  NO. (USE  0  TO  FINISH):') 

READ  ( COM_UNIT, 111 )  IB 

WRITE  ( LOGUNIT, 111 )  IB  !  log  the  input 

IF(  IB.LT.l  .OR.  IB.GT.NBUDGO  TO  110 
DOMAIN  -  DOM( IB) 

IB-IB>1 

WRITE) lOUT, 104 ) 

FORMAT)'  MULTIPLIER-') 

READ  (COM_UNIT,108)  AMULT 

WRITE  ( LOG_UNIT, 108 )  AMULT  !  log  the  input 

DO  105  J-1, IHDSZ 
HEAD) J,1)-HEAD( J,1B) 

DO  107  J-1,IS1ZE 
P(J,1)-P(J,1) +AMULT*  P ( J , I B ) 

IDS) 1 )-IDS( IB) 

GO  TO  102 
FORMAT  (F10.2) 

CONTINUE 
FORMAT  (15) 

RETURN 


ROUTINE  TO  READ  FROM  DISK 

ENTRY  REA(HEAD,  LINEl,  L1NE2,  LINE3) 

DOMAIN-'TIME' 

WRITE) lOUT, 112 ) 

FORMAT)'  INPUT  FILE  7  ',$) 

READ  (COM_UNIT,114)  FILNM 

WPTTP  (lOG  UNIT, 114)  FTLNM  I  t  t,„  i...  ... 

IF) ( PROC_FLAG) .AND. ( FILNM  . EQ .' NAME . DAT' ) )  THEN  !  il  »  pi  ' 

DO  I-1,ISIZE  1  Make  a  dummy  tir'^  . 

P(I,1)-10 

END  DO 
IDS) 1 )-0 

WRITE) lOUT, 113 )  I  tall  user  what  is  going  on 

RETURN 

END  IF 

FORMAT)'  This  is  a  duuy  file  for  PROC  definition.') 

FORMAT(A50 ) 

OPEN(UNIT-l ,ERR-115, FILE-FILNH,TyPE-'OLD' , READONLY, 

+  lOSTAT-  FOR_RETCODE, DEFAULTFILE-DEF, FORM-' UNFORMA TTED ' ) 

IPGF-0 
GOTO  116 

IF  (FORRETCODE  .EQ.  29)  THEN  1  if  file  not  found 
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CALL  LIB$SIGNAL(MAG_riLNOTFOU)  1  print  error 

RETURN 

ELSE 

CALL  LIB$SIGNAL(MAG_COM)  1  Print  an  error  messag 

RETURN 

END  IF 

116  READ( 1 , IOSTAT-FOR_RETCODE. end-162)  ( HEAD ( I . 1 ) , I-l , 2 56 ) 

IF  (FOR_RETCODE  .EQ.  35)  THEN 

CALL  LIB$SIGNAL(MAG_1NPF0R) 

RETURN 

END  IF 

WRITE( IOUT,117)  LINEl 
WRITE(IOUT,117)  LINE2 
WRITE! lOUT, 117 )  LINES 

117  FORMAT! X, 3 0A2 ) 

DO  77812  J-1,ISIZE 
77812  P!J,l)-0. 

READ! 1 , END-162)  !P! J ,1 ) , J-l.ISIZE) 

IDS! 1 )-0 
GO  TO  164 

162  WRITE! TERM_UNIT,*)  'ERROR  END  OF  FILE  ENCOUNTERED  AT  ' , I 

164  CLOSE!UNIT-l ) 

165  FORMAT  !I5) 

RETURN 


TO  DELETE  A  FILE  ON  DISK 


ENTRY  DEL 

WRITE! lOUT, * )  'ENTER  THE  FILE  NAME' 

READ  (COM_UNIT, 114 )  FILNM 

WRITE  !LOG  UNIT, 114)  FILNM  1  log  the  input 

OPEN! UNIT-T, ERR-166 , FILE-FILNM,TYPE-'OLD' , 

+  DEFAULTF1LE-DEF,DISP-'DELETE' ) 

CLOSE(UNIT-l) 

GOTO  167 

66  CALL  LIB$S1GNAL!MAG_C0M)  1  Print  an  error  messag 

67  RETURN 


THIS  ROUTINE  PLOTS  A  WAVEFORM  FROM  A  BUFFER 


FNTRV  pr.o 
WRITE  !I0UT,186) 

READ  !COM_UNIT,187)  ISEGl 
WRITE  !LOG_UNIT,187)  ISEGl 
ISEG-ISEGl+1 

WRITE! lOUT, *) 'DO  VOU  WANT  A  NEW  WINDOW? ( 1-Y, O-NO ) ' 

READ  !COM_UNIT,*)  IANS 

WRITE  |LOG_UNIT,*)  IANS 

IF  !IANS.NE.l)  GO  TO  18500 

WRITE! lOUT, *) 'WINDOW  VALUES  ARE  NOW:' 

WRITE  !IOUT,4247)  WINl ,WIN2 ,WIN3 ,WIN4 

WRITE! lOUT, *) 'START  ELEMENT  FOR  PLOT! O-NO  CHANGE):' 

READ  !COM_UNIT,*)  lANSl 

WRITE  !LOG_UNIT,*)  lANSl 

IF! lANSl .NE.O)  IWINl-IANSl 

WRITE! IOUT,» ) 'END  ELEMENT  FOB  PLOT  ! O-NO  CHANGE):' 
READ  !COM  UNIT,*)  lANSl 


!  log  the  input 


!  log  the  input 


!  log  the  input 
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£ 


A 

A 


18500 


18600 


18601 

4247 


18610 


h 


r 


186 

187 

C 

C 

C 

C 


log  the  input 


log  the  input 


WRITE  (LOG_UNIT,*)  lANSl  !  log  the  input 

IF( lANSl.NE.O)  IWIN2-IANS1 

WRITE( lOUT, * ) 'MAX  PLOT  VALUE  (+1024  NORMAL, O-NO  CHANGE):' 

READ  (COM_UNIT,*)  ANS 

WRITE  (LOG_UNIT,*)  ANS  !  log  the  input 

IF  (ANS.NE.O)  W1N4-ANS 

WRITE( lOUT, * ) 'MIN  PLOT  VALUE  (-1024  NORMAL, O-NO  CHANGE):' 

READ  (COM_UNIT,*)  ANS 

WRITE  (LOG  UNIT,*)  ANS  ! 

IF(ANS.NE.I5)  WIN3-ANS 
WINl-FLOAT( IWINl ) 

WIN2-FLOAT( IWIN2 ) 

WRITE( lOUT, * ) 'DO  YOU  WANT  NEW  AXES? ( 1-Y , O-N )  ' 

READ  (COM  UNIT,*)  ANS 

WRITE  (L06_UNIT,*)  ANS  ! 

IF  (ANS.EQ.O)  GO  TO  18600 

WRITE( lOUT,* )' INPUT  Y-AXIS  MAX, MIN, AND  TIK  MARK  INTERVAL’ 

READ  (COM_UNIT,*)  YMAX , YMIN , YTMI 

WRITE  (LOG_UNIT,*)  YMAX , YMIN , YTMI  !  log  the  input 

WRITE( lOUT, * ) 'Y  AXIS  LEGEND' 

READ  (COM_UNIT, 18601)  YLAB 

WRITE  ( LOG_UNIT, 18601 )  YLAB  !  log  the  input 

WRITEdOUT,  *)' INPUT  X-AXIS  MAX, MIN, AND  TIK  MARK  INTERVAL' 

READ  (COM_UNIT,*)  XMAX , XMIN, XTMI 

WRITE  (LOG_UNIT,*)  XMAX , XMIN , XTMI  !  log  the  input 

WR1TE( lOUT, * ) 'X  AXIS  LEGEND' 

READ  (COM_UNIT, 18601)  XLAB 

WRITE  ( LOG_UNIT, 18601 )  XLAB  !  log  the  input 

WRITE( lOUT, * )  'INPUT  TITLE  FOR  PLOT’ 

READ  (COM  UNIT, 18601)  TITLE 

WRITE  (L06  unit, 18601)  TITLE  !  log  the  input 

FORMAT  (A6U) 

FORMAT('  STRT',r8.2,'  END  ' , P8 . 2 , '  MIN  Y',Fe.2,'MAX  Y',F8.2) 

YINC-WIN4-WIN3 

YZERO-WIN3+YINC/2 . 

XINC-W1N2-W1N1 

IXINC-INT(XINC) 

DRAW  THE  CURVE  USING  A  SQUARE  7  IN  PLOT 
CALL  VPLOTS( 0,0,0) 

PSIZE-7 . 

CALL  PLOT(0.5,5.5,-3) 

DO  18610  I-1,IXINC 
X-I*PSIZE/XINC 

Y-PSIZE*(P( IWINl +I , ISEG)-YZER0)/YINC 
CALL  PLOT(X,Y,2) 

CALL  PLOT( 0,-PSIZE/2. ,-3) 

XDS-(XMAX-XMIN)/PSIZE 

XSP-PSIZE*XTMI/(XMAX-XMIN) 

CALL  FFAXIS( 0. ,0. , »REr ( XLAB ) , -60 , PSIZE , 0 . , XMIN, XDS , XSP , 1 , 1 . , 1 ) 

YDS- (YMAX- YMIN) /PSTZE 
YSP-PSIZE*YTMI/( YMAX-YHIN) 

CALL  FFAXIS(0. ,0. , %REF(YLAB) ,60,PSIZE,90. , YMIN , YDS , YSP , 1 , 1 . , 1 ) 

CALL  TEXT ( 0 . 5 , PSIZE+0 . 5 , 0 . 1 , TITLE, 0 . ) 

CALL  PLOT(0,0,999) 

CALL  PLOTNOW( IMSG) 

IF(IMSG.EQ.O)  THEN 

CALL  LlBSSIGNAL(MAG_COM)  1  Print  an  error  message 

RETURN 

END  IF 

FORMAT  ( IX, 'WAVEFORM  BUFFER  NO. ?( O-HAIN) ' ) 

FORMAT  (15) 

RETURN 
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THIS  ROUTINE  READS  DATA  FROM  A  BASIC  SCATTERING  CODE  CALCULATION 


ENTRY  BSC 

CALL  BSCREA(P,COM_UNIT) 
FLTP-1 . 

RETURN 


THIS  ROUTINE  ATTACHES  A  SMOOTH  CUTOFF  TAIL  TO  THE  SPECTRUM 
IN  ORDER  TO  REDUCE  THE  GIBBS  PHENOMENOM 


FROM  LEEPER'S  COSINE  ROLOFF  18FEB83  IN  [ LEEPER . X ] DUMPZ . FOR 

ENTRY  WND 
WRITE( lOUT, * ) 

'INPUT  HARMONICS ; START, END, TYPE  OF  WINDOW' 

WRITE( lOUT, * ) 

' TYPE : 0-HANNING , 1-HAMMING , 2-GAUSSIAN ,  0 , N , M , -TEST ' 
READ(COM_UNIT,*)  IFS,IFE,ITYP 

WRITE( LOG_UNIT, * )  IFS,IFE,ITYP  1  Lc 

IF( IFS.GE. IFE)  THEN 

CALL  LIB$SIGNAL(MAG  COM)  I  Print  an  error  message 

RETURN 

END  IF 

IF(ITYP.E0.2)  THEN 

WRITE( IOUT,»)  'INPUT  ALPHA' 

END  IF 

IF( ITYP.EQ.2)  THEN 

READ(COM_UNIT,*)  ALPHA 

WRITE(LOG_UNIT,»)  ALPHA  !  Lo 

END  IF 

PI-3.1415927 

NW-IFE-IFS+1  1(  -  t  POINTS,  ODD  OR  EVEN  ) 

XN2-NW/2 . 

IF(XN2.NE.INT(XN2) )GO  TO  178 
IFE-IFE+1 
GO  TO  180 
IN2-XN2 

T r ( I rc . MP . n iro  TO  190 
O’)  tnj  1-1,13 
P( 1+lS, l)-0. 

P(I,1)-102.4  (FILLS  P  ARRAY  WITH  UNITY  FOR  TEST 

GO  TO  179 
DO  198  1-IFS,IFE 
N-I-IFS 

WRITE! lOUT, * )  I 

WFACT-0.5*(1.-COS(2*PI*N/NW) )  (HANNING 

IF(ITYP.EQ.l)  WFACT-.54-.46*C0S(2*PI*N/NW)  (HAMMING 

IF( ITYP.NE.2)  GO  TO  1195 
N-N-IN2 

AVAR-- . 5* ( ALPHA*N/IN2 ) **2 
WFACT-EXP(AVAR)  (GAUSSIAN 

IF(WFACT.EQ.O)  WFACT-3.E-6 

P(I,1)-P(I,1)/10. 24-10.  (-20LOG(ORIGINAL  VOLTAGE) 

P(I,l)-10.**(P{I,l)/20. )  (-LINEAR  VOLTAGE 


(  Log  tlie  input 


Log  tlie  input 


a 


.V 


V 


.s 


I 


P( I , 1 )-P( I , 1 )*WFACT  IWEIGHT  WITH  SHIFTED  COSINE 

198  P(I ,1)-10.24*(10.+20.*ALOG10(P(I,1» ) ) 

IFE-IFE+1 

IFS-IFS-l 
DO  199  I-IFE,IS 

199  P( I , 1 )--1024 . 

DO  188  I-1,IFS 

188  P(I,1)--1024. 

RETURN 

C******»fill  remaining  P(K,1)  with  -llOdB  -  (-1024/10.24 )-10. 

C 

C 

C 

C 

C 

C  GATE  IN  THE  FREQ  DOMAIN  ROUTINE: FGT 

C  STARTS  AT  HARMONIC  N 

C  ATTENUATES  AT  X  DB  PER  HARMONIC 

C 

c 


191 

192 

193 

194 

195 

196 
200 
2"! 


ENTRY  FGT 
WRITE  (IOUT,191) 

FORMAT!'  HARMONIC  FOR  START  HIGH  FREQ  CUTOFF  ',$) 
READ  ( COM_UNIT, 200 )  IHAR 
WRITE  (LOG_UNIT,200)  IHAR 
WRITE( lOUT, 192 ) 

FORMAT ( '  HIGH  FREQ  ROLL-OFF  IN  DB  PER  HARMONIC', $) 
READ  (COM_UNIT,201 )  XDB 
WRITE  { LOG_UNIT, 201 )  XDB 
WRITE  (IOUT,193) 

FORMAT ( '  HARMONIC  FOR  START  OF  LOW  FREQ  CUTOFF  ',$) 
BEAD  (COM  UNIT, 200)  ILHAR 
WRITE  (LOg_UNIT,200)  ILHAR 
WRITE  (IOUT,194) 

FORMAT!'  LOW  FREQ  ROLLOFF  IN  DB  PER  HARMONIC  ’,$) 

READ  (COM  UNIT, 201)  YDB 

WRITE  (LOg_UNIT,201)  YDB 

DO  195  I-IHAR,IS 

FHAB-FLOAT( I-IHAR) 

P( I , 1 )-P( I , 1 )-10 .24*XDB»FHAR 

IF( P( I , 1 ) .LT.-1024 . )  P( I ,1 )— 1024. 

CONTINUE 

DO  196  I-l, ILHAR 
FHAR-FLOAT{ ILHAR-I ) 

P( I , 1  )-P( I , 1 )-10 .24*YDB*FHAR 
IF(P(I,1).LT.-1024.)P(I,1 )— 1024. 

CONTINUE 
FORMAT ( I  5 ) 

FORMAT! FIO . 2 ) 

RFTtIPN 


I 


1 


C 

C 

c 

c 


210 


TIME  DOMAIN  INTEGRATE  ROUTINE 

ENTRY  INTEG 

IF! IDS! 1 ) .NE.O )  THEN 

CALL  LIBSSIGNAL!MAG_COH) 
RETURN 

END  IF 

PIZ-2 . ‘PI/ISIZE 
P! 1 , 1 )-P! 1 , 1 ) *PIZ 
DO  210  N-2,ISIZE 
P!N,1)-P(N-1,1)+P(N,1)*PIZ 
RETURN 


Print  an 


log  the  input 


log  the  input 


log  the  input 


log  the  input 


error  message 


C 


TIME  DOMAIN  DIFFERENTIATE  ROUTINE 


ENTRY  DIF 

IF( IDS(l) .NE.O)  THEN 

CALL  LIB$S1GNAL( MAG_COM)  !  Print  an  error  messaqp 

RETURN 

END  IF 

PE-P( ISIZE, 1  ) 

DO  220  N>2, ISIZE 
NN-ISIZE+2-N 

P(NN,1)-(P(NN,1)-P(NN-1,1) )*ISIZE/(2.*PI) 

P(1,1)-(P(1, 1 )-PE)*ISIZE*0. 5/PI 
RETURN 


STATEMENT  550  32583 


RELABLE  COMMAND 


ENTRY  RLB 
WRITE  {IOUT,225) 

WRITE  (IOUT,231)  ITITLE 
FORMAT(/,'  OLD  TITLE  BLOCK:  ') 
WRITE  (IOOT,230) 

FORMAT(/,'  NEW  TITLE  BLOCK:  ') 
READ  ( COM_UNIT, 231 )  ITITLE 
WRITE  (LOG  UNIT, 231)  ITITLE 
FORMAT( 50AT) 

RETURN 


1  sho  the  present  label 
1  prompt 

1  read  the  new  label 
!  log  the  input 


NORMALIZING  TIME  DOMAIN  DATA  TO  ZERO  MEAN 


ENTRY  NOR 
RMEAN-O. 

DO  3111  I-l, ISIZE 
RMEAN-P( I , 1 )+RMEAN 
RMEAN-RMEAN/ISI2E 
DO  3112  I-l, ISIZE 
P( I , 1 )-P( I , 1 )-RMEAN 
RETURN 


'-Trap  thf  main  Pt'FFFR 


ENTRY  CLR 
DOMAIN  -  '  ' 

LINEl(l)  -  0 
DO  3113  I-l, ISIZE 
P(I,l)-0. 

RETURN 

POINT  SMOOTH  ROUTINE 


ENTRY  PSM 

ir{FLTP.EQ.l)  GO  TO  4051 
WRITE  (IOUT,4001) 


4001 


,  V 

.■j' 


V, 


iN- 


's. 


f 


40C2 

4003 


4010 

4051 

4052 


4053 


FORMAT ( '  FIRST  BAD  POINT  INTEGER-  ',$) 

READ  (COM_UNIT,4003)  IQ 
WRITE  (LOG_UNIT,4003)  IQ 
WRITE  (IOUT,4002) 

FORMAT('  LAST  BAD  POINT  INTEGER-  ',$) 

FORMAT  (15) 

READ  ( COM_UNIT, 4003 )  IR 
WRITE  (LOG  UNIT, 4003)  IR 
1F( IQ. LT. 31  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

RETURN 

END  IF 

IF(ia.LT.3)  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

RETURN 

END  IF 

IF( IQ.GT. ISIZE-2 )  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

RETURN 

END  IF 

IF( IR.GT. ISIZE-2  )  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

RETURN 

END  IF 

If  (IR.LT.IQ)  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

RETURN 

END  IF 

PSQ-P( IQ-1,1) 

PSN-IR-lQ+2 

PSS1-P(10-1,1)-P(IQ-2,1) 

PSS2-P(IB+2,1)-P(IR+1,1) 

psdi-p(ir+i,i)-p(:q-i,i) 

PSC-3 . *PSDl/( PSN*PSN)- ( 2 . *PSS1+PSS2 )/PSN 

PSD-(PSS2+PSS1 |/(PSN*PSN)-2.*PSOl/(PSN»PSN*PSN) 

PSA-PSQ 

PSB-PSSl 

IPSN-IR-IQ+l 

DO  4010  1-1,IPSN 

P( lQ+IPSN-1 , 1 )-PSA+PSB*I+PSC*I*I+PSD*I»l*I 

GO  TO  4075 

WRITE  (IOUT,4052) 

FORMAT ( ’  FIRST  BAD  HARMONIC-  ',$) 

READ  (COM_UNIT, 4003 )  IQ 
WRITE  ( LOG_UNIT, 4003 )  IQ 
WRITE  (IOUT,4053) 

FORMAT('  LAST  BAD  HARMONIC-  ',$) 

READ  (COM  UNIT, 4003)  IR 
WRITE  (LOG  UNIT, 4003)  IR 
IF  (IO.LT.T)  THFN 

CALL  LIBSSIGNAL(MAG  COM) 

RETURN 

END  IF 

1F(IR.LT.3)  THEN 

CALL  LIBSSIGNAL(MAG  COM) 

RETURN 

END  IF 

lr( IQ.GT. ISM2)  THEN 

CALL  L1B$SIGNAL(MAG  COM) 

RETURN 

END  IF 

IF( IR.GT. ISM2)  THEN 

CALL  LIB$SIGNAL(NAG  COM) 

RETURN 

END  IF 

IF(IR.LT.IO)  THEN 


!  log  t)ie  input 


I  log  the  input 


1  log  the  input 

!  log  the  i npn f 
I  Print  an  i  i  ■  i  "  p 

!  Print  an  error  me stage 


I  Print  an  error  message 
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CALL  LIB$SIGNAL(HAG_COM) 

RETURN 

END  IF 
RT2-SQRT( 2  .  ) 

DO  4056  N-IQ-2,IR+2 
NN-IS+N 

AMP-( P(N,1)/10.24)-10. 

PHAS-PI»P(NN,1)/512. 

A(N)-CMPLX(COS(  PHAS)  ,Si;,-(PHAS)  )*RT2*10.**(AMP*.05) 

WRITE! lOUT, • )  N,AMP,PHAS,A(N) 

ASO-A( lQ-1) 

ASN-IR-IQ+2 
ASS1-A( IQ-1 )-A( IQ-2 ) 

ASS2-A(  IR+2  )-A(  IR-t-1  ) 

ASD1-A( IR+1 )-A( IQ-1 ) 

ASC-3*ASDl/(ASN*ASN)-( 2. •ASS1+ASS2)/ASN 

ASD-( ASS2+ASS1 )/(ASN»ASN)-2 .•ASD1/(ASN«ASN*ASN) 

ASA-ASQ 

ASB-ASSl 

WRITE! lOUT, * )  ASA,ASB,ASC,ASD 

IASN-IR-IQ+1 

DO  4060  I-1,IASN 

IN-IQ-l+I 

FI -FLOAT! I ) 

A! IN)-ASA+ASB*FI+ASC*FI*FI+ASD*Fl*ri*FI 
WRITE! lOUT, * )  A!IN) 

X-CABS!A! IN) ) 

IF!X.LT.l.E-20)X-l.E-20 
P!IN,1)-10.24*! 10 .+20. *ALOG(X*RT2) ) 

1F!X.EQ.0)  P|IN+IS,l)-0. 

IF!X.NE.0)P! IN+1S,1 )-512.*ATAN2(AIMAG!A( IN) ) .REAL! A!  IN) ) )/PI 

CONTINUE 

RETURN 


YSHIFT  COMMAND 


ENTRY  YSH 
WRITE  !IOUT,251) 

FORMAT!'  YSHIFT  IN  UNITS ! -2O40<Y<2O40 )-  ',$) 

READ  !COM_UNIT,252 )  lY 

WRITE  !LOG  UNIT, 252)  lY 

FORMAT  !I5T 

y-pr  OAT! lY ) 

DO  255  I-1,ISIZE 

P!I,1)-P!I,1)+Y 

RETURN 


1  log  the  input 


ROTATE  COMMAND 


ENTRY  ROT 
WRITE  !IOUT,260) 

FORMAT  !'  ROTATE  BY  INCREMENTS  OF: 

READ  !COM_UNIT,26] '  IX 

WRITE  !LOG  UNIT, 261'  IX 

FORMAT  !I5T 

DO  265  1-1,ISIZE 

J-I+IX 


1  log  the  input 


265 


I 

.V 

Vi 


267 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IF( J.GT. ISIZE) J-J-ISIZE 
IF( J.LT.1)J-J+ISIZE 
Q(J)-P(I,1) 

DO  267  I-l, ISIZE 
PtI,l)-Q(I) 

RETURN 


ROUTINE  TO  DOWNCONVERT  A  SPECTRUM  TO  A  DUAL-POLARITY  ENVELOPE 
USER  CHOOSES  THE  HARMONIC  NUMBER  TO  CALL  HIS  "CENTER  FREQUENCY" 
WHICH  GETS  CONVERTED  TO  THE  DC  TERM 

WORKS  ON  A  FREQUENCY  SPECTRUM  IN  MAIN  BUFFER , RESULT  TO  MAIN  BUFFER 


ENTRY  DCV 

IF  (FLTP.EQ.O)  THEN 

CALL  LIB$SIGNAL(MAG_COM)  t  Print  an  error  message 

RETURN 


END  IF 


WRITEdOUT,*)  'INPUT  THE  HARMONIC  NUMBER 

TO  BE  MOVED  TO  DC' 

READ  (COM  UNIT,*)  IHRDC 

WRITE  (LOG  UNIT,*)  IHRDC 

!  log  the  input 

IF  (IHRDC.  GTdS)  THEN 

CALL  LIB$SIGNAL(MAG  COM) 

1  Print  an  error  message 

RETURN 

END  IF 

•  ■ 

DO  5110  N-1, ISIZE 

5110 

Q(N)-P(N,1) 

, 

DO  5120  N-1, IS 

0(N)-(0(N)/10.24)-10. 

■ 

5120 

Q(N+IS)-PI*Q(N+IS)/512. 

DO  5190  N-1,IS 


M-N-1 


-C" 

K 


y 

yi 


BCX-O . 

CCX-0. 

P(N+IS, 1 )-0. 

IF( IHRDC+H.GT.IS)  GO  TO  5130 

BCX-CMPLX(COS(Q( IHRDC+M+IS) ) , SIN( Q( IHRDC+M+I S ) ) ) * 1 0 . * * ( Q ( I HRDC+M ) * . 05 ) 
5130  IF( IHRDC-M.LT.l )  GO  TO  5140 
AEXP-10.**(0.05*Q(IHRDC-H) ) 

CCX-CMPLX{COS(Q( IHRDC-M+IS) ) ,-SIN(Q( IHRDC-M+IS) ) ) *AEXP 
5140  DCX-BCX+CCX 

X-CABS(DCX) 

IF(X.LT.l.E-20)  X-l.E-20 

P(N, 1 )-10.24*( 10.+20.*ALOG10(X) ) 

T F( X  r F . 1 . E-70 )  GO  TO  5190 

r  '"1.’  .  ‘ATAN2(A1IIAG(DCX)  ,REAL(DCX)  )/PI 

cion  rriNTI  NUE 

RETURN 
C 
C 
C 
C 
C 
C 

c 

C  ROUTINE  TO  SET  UP  THE  PARAMETERS  FOR  A  TWO  DIMENSIONAL  IMAGE 

C 

C 

c 

ENTRY  CMI 

WRITEdOUT,*)  ’NUMBER  OF  W  TIME  DOMAIN  WAVEFORMS  TO  BE  USED?’ 


59 


oooonooo 


u  "'j  w'\r«  J  •u  •u  •Ml 


1 11 IL  ll 


5503 

5505 

5507 

5510 


5513 


5520 


5523 


5530 


READ  (C0M_UNIT,*)  NUMB ( 1 ) 

WRITE  (LOG  UNIT,*)  NUMB(l) 

DO  5510  I-T,NUMB(1) 

WRITE  (IOUT,5503)  I 

FORMAT!'  BUFFER  NUMBER  FOR  W  FILE  t  ',I5) 

READ  (COM_UNIT,*)  BUFN(1,I) 

WRITE  (LOG  UNIT,*)  BUFN(1,I) 

BUFN( 1 , I )-T+BUFN(  1 , 1  ) 

WRITE  (IOUT,5505)  I 

FORMAT!'  LOOK  ANGLE  IN  DEGREES  FOR  FILE  t  ',I5) 
READ  ( COM_UNIT, * )  ANG(1,I) 

WRITE  (LOG_UNIT,*)  ANG(1,I) 

ANG( 1 , I )-ANG( 1, I )*PI/180. 

WRITE  (IOUT,5507)  I 

FORMAT!'  CENTER  ELEMENT  NUMBER  FOR  FILE  I  ',15) 
READ  !COM_UNIT,*)  CNTR ( 1 , I ) 

WRITE  !LOG_UNIT,*)  CNTR! 1,1) 

WRITE! lOUT, * )  'NUMBER  OF  HH  TIME  DOMAIN  WAVEFORMS 
READ  !COM_UNIT,*)  NUMB! 2) 

WRITE  (LOG  UNIT,*)  NUMB! 2) 

DO  5520  I-T,NUMB!2) 

WRITE  !I0UT,5513)  I 

FORMAT!'  BUFFER  NUMBER  FOR  HH  FILE  I  ',I5) 

READ  !C0M_UNIT,*)  BUFN!2,I) 

WRITE  (LOG  UNIT,*)  BUFN!2,I) 

BUFN!2,I)-T+BUFN!2,I) 

WRITE  !IOUT,5505)  I 
READ  !COM_UNIT,*)  ANG!2,I) 

WRITE  !LOG_UNIT,*)  ANG!2,I) 

ANG! 2, I )-ANG! 2, I )* PI/1 80. 

WRITE  !IOUT,5507)  I 

READ  (COM_UNIT,*)  CNTR! 2, I) 

WRITE  !LOG_UNIT,*)  CNTR! 2, I) 

WRITE! lOUT, * )  '  NUMBER  OF  HV  TIME  DOMAIN  WAVEFORMS 
READ  !COM_UNIT,*)  NUMB ! 3 ) 

WRITE  (LOG  UNIT,*)  NUMB! 3) 

DO  5530  1-T,NUMB! 3) 

WRITE  !IOUT,5523)  I 

FORMAT!'  BUFFER  NUMBER  FOR  HV  FILE  I  ',I5) 

READ  (COM  UNIT,*)  BUFN!3,I) 

WRITE  !L05  unit,*)  BUFN!3,I) 

BUFN! 3 , I )-T+BUFN!  3,1) 

WRITE  !IOUT,5505)  I 
READ  !COM_UNIT,*)  ANG!3,I) 

WRITE  !LOG_UNIT,*)  ANG(3,I) 

ANG! 3, I ) -ANG! 3, I )*P 1/180. 

WRITE  !IOUT,5507)  I 

READ  !COM_UNIT,*)  CNTR!3,I) 

WRITE  !LOG_UNIT,*)  CNTR! 3, I) 

RETURN 


1  log  the  input 


!  log  the  input 


I  log  the  input 


1  log  the 
TO  BE  USED?' 


input 
1  log  the  input 


!  log  the  input 


I  log  the  input 


1  log  the 
TO  BE  USED?' 


input 
!  log  the  input 


1  log  the  input 

!  log  the  input 

!  log  the  input 


ROUTINE  TO  FORM  THE  IMAGE  FOR  ONE  POLARIZATION  FROM  TIME  DOMAIN 
WAVEFORMS  AS  SET  UP  BY  THE  CMI  COMMAND 


ENTRY  IMG 


WRITE! lOUT, * )  'SIZE  OF  THE  IMAGE, 
READ  !COM_UNIT,*)  IMSIZE 
WRITE  (LOGUNIT,*)  IMSIZE 
FIMINC-FLOAT! IMSIZE/100) 
PNUM-FIMINC  I 

IF  !M0D!PNUM,2)  .EQ.  0)  THEN 


(1  TO  4096)-?'  I  get  the  window  size 

I  log  the  input 
I  the  size  of  the  increment 
figure  out  the  number  of  harmonics  to 
I  be  averaged  for  an  array  value 


PNUM-PNUM+1 


I  This  number  must  be  odd 


i 


«y. 


A, 


V 

y 


i 

I 


'  «> 

■U: 


.y 


I  when  PNUM  is  3 


the  weights  of  the  other 
! two  numbers  are  0 


get  the 


END  IF 

IF  (PNUM  .LT.  3)  THEN 
PNUM-3 

END  IF 

WRITE! lOUT, * )  'LOOK  ANGLE  OF  THE  I MAGE , DEGREES? ' 
rotation  angle  of  the  image 
READ  (COM_UNIT,»)  ALANGL 
WRITE  (LOG_UNIT,*)  ALANGL 
BLANGL-ALANGL*P1/180 . 

WRITE! lOUT, * )  'POLARIZATION  OF  THE  I  MAGE , ! 1-W , 2-HH , 3-HV , 4 -ALL ) 
polarizations  will  be  used? 


log  the  input 


whi  ch 


!  log  the  input 

if  not  using  all  polarizations  then 
just  do  the  individual  pol . 

I  figure  out  the  X  coordinate 

1  figure  out  the  Y  coordinate 

!  the 


the  Y 


1  clear  the  image  array 


READ  !COM_UNIT,*)  IIPOL 
WRITE  !LOG_UNIT,*)  IIPOL 
IPOL-I IPOL 

IF! I IPOL.NE. 4 )  GO  TO  5660  I 

DO  5699  IIPOL-1,3  1 

5660  DO  5690  1-1,100 

X-! 1-50 . 5 ) ‘FIMINC 
DO  5680  J-1,100 
Y-  !  J-50 . 5  )  *FIMINC 
XX-X*COS! BLANGL)+Y*SIN!BLANGL) 

C  X  coordinate  taking  the  rotation  angle  into  account 
YY--X*SIN! BLANGL)+Y*COS!BLANGL) 

C  coordinate  taking  the  rotation  angle  into  account 
ARRAY! IIPOL, I , J )-0 . 

DO  5670  K-1 , NUMB! IIPOL) 

XXX-XX*COS!ANG! IIPOL, K) )+YY*SlN!ANG! IIPOL, K) )  !  the  X  coordinate 

YYY--XX*SIN! ANG! I IPOL, K ) )+YY»COS!ANG( IIPOL, K) )  1  the  Y  coordinate 

IPOINT-YYY+CNTR! IIPOL, K)  1  the  point  of  intersection 

IF! IPOINT. LT. 0 )  IPOINT-IPOINT+ISIZE  I  the  time  domain  waveform 
C  repeats  throughout  time 

IF! I  POINT. GT. I  SIZE)  IPOINT-IPOINT-ISIZE 
VALUE-0 

DO  ELEMENT-dPOINT-! !PNUM-l)/2) 1 .(IPOINT+((PNUH-l)/2) )  !  take 

C  the  points  around  the  center  point 
ELEM-ELEMENT 

WEIGHT-l+COS! !2»(ELEH-1P01NT)»P1)/(PNUH-1) )  !use  a  weighted 

C  shifted  cosine  average 

1F!ELEM.LT.0)  ELEM-ELEM> I  SIZE  1  the  time  domain  waveform 

C  repeats  throughout  time 

IF! ELEM.GT. ISIZE)  ELEM-ELEM-ISIZE 
PVALUE-WEIGHT«P!ELEM,BUFN! IIPOL, R) ) 

VALUE-VALUE+PVALUE 


END  DO 

VALUE-VALUE/! PNUM-1 ) 

ARRA Y ! 1 1  POL , I , J ) -ARRAY ! 1 1  POL , I , J ) +VALUE/NUHB (IIPOL) 

5670  CONTINUE 

5680  CONTINUE 

ceon  roMTTNUE 

5GOU  '.''Mtinije 

C  WRITE  THE  IMAGE  ARRAY  TO  AN  OUTPUT  FILE 

WRITE! lOUT, *)' DO  you  want  to  store  the  image?  Y-l,N-0  ' 

READ  (COMUNIT,*)  ANS 
WRITE! LOG_UNIT, * )ANS 
I F! ANS. EQ.O) RETURN 

WRITE! lOUT, * )  'Enter  the  output  file  name' 

READ ! COM_UN I T , 1 0 ) FNAME 
WRITE!LOG  UNIT,10)FNAME 
10  FORMAT!A4C) 

C 

WRITE! lOUT, • )  'Enter  the  freq.  increment  of  the  signal  in  MHz' 
READ! COM_UNlT, •  )  FER 
WRITE! LOG_UNIT, * )FER 

PER-1 ./FER*1000  (Period  of  the  time  signal 

C 
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I 


•< 


'.•.v 


OPEN( UNIT-31 ,FILE-FNAME, TYPE-' NEW' , 
FORM- ' UNFORMATTED ' ) 


WRITE( 31 ) IMSIZE.PER 

DO  99  1-1,100 

DO  99  J-1,100 

WRITEOl)  ARRAY(  IPOL,  I ,  J  ) 

CONTINUE 

RETURN 


THIS  ROUTINE  PLOTS  AN  ISOMETRIC  VIEW  OF  A  SINGLE  POLARIZATION 
TWO-DIMENSIONAL  TARGET  IMAGE  ON  THE  PLOT  DEVICE  IN  ISOMETRIC 
FORM  WITH  NO  SHADOWING 


ENTRY  PIM 

WRITE( lOUT, * )  'TYPE  IN  A  POLARIZATION  INDEX ( l-W , 2-HH , 3-HV ) ' 

READ  (COM_UNIT,*)  IIPOL 

WRITE  (LOG_UNIT,*)  IIPOL  !  log  the  input 

IF( IIPOL.GT. 3 )  THEN 

CALL  LIB$SIGNAL( MAG_COM)  1  Print  an  error  message 

RETURN 

END  IF 

IF(IIPOL.LT.l)  THEN 

CALL  L1B$SIGNAL(MAG_C0M)  I  Print  an  error  message 

RETURN 

END  IF 

CALL  VPLOTS( 0,0,0) 

CALL  PLOT(4. ,5. ,-3) 

DO  5790  1-1,100 
DO  5790  J-1,100 
X-(J-20.5)*0.1 
Y-(I-20.5)»0.1 
XX-0.866*(X+Y) 

YY-0.5»( Y-X)+ABRAY( 1 1  POL , I . J )/256 . 

IF(J.EO.l)  CALL  PLOT{XX, YY,3) 

IF(J.NE.l)  CALL  PLOT(XX, YY,2 ) 

CONTINUE 

CALL  PLOT(0, 0,999) 

CALL  PLOTNOW( IMSG) 

RETURN 


ENTRY  FRDC 

CALL  FRD( P,COM_ONIT) 

FLTP-1 

RETURN 


THIS  COMMAND  READS  NEW  MEASUREMET 


ENTRY  RDFL(HEAD,  LINEl,  LINE2,  LINE3) 


DOMAIN-' FREQ' 

CALL  FRDl ( P, EFLAG) 

IF  (EFLAG)  THEN 

EFLAG  -  .FALSE. 

ELSE 

FLTP-1 

END  IF 
RETURN 


!  if  an  error  occurred 
1  reset  error  flag 


READ  THE  CARD  FILE 


ENTRY  CARDC 
CALL  CARD(P) 
RETURN 


SHOW  BUFFER  ROUTINE: 

This  routine  shews  the  contents  of 
all  of  the  buffers. 


ENTRY  SHO_BUF 
DO  1-1,40 

IF  (BUFFERS! I , 1 )  .NE.  0)  THEN 

WRITE! lOUT, 6105)  (I,  DOM(I)) 

WR1TE( IOUT,6110>  ( BUFFERS) 1 , J ) ,  J-1,30) 

END  IF 

END  DO 

FORMAT! '0' ,  14,  3X,  A4 ) 

FORMAT! 4X,  60A2) 


********<r*<<l 


*•*  *•  It  It 


SUBROUTINE  FRD ! P , COM_UNIT ) 

COMPLEX  BA(201) ,KLA,CA!2049) 

CHARACTER*20  INFILE 
INTEGER* 4  COM_UNIT 
REAL  FA(2049) ,P!2048,11  ) 

DO  3410  1-1,1024 
P!I+1024,l)-0. 

P(I,  1)  — 1024. 
no  i.i ,2049 

CA( I )-0 . 

FA! I )-0. 

WRITE! lOUT, * )  'INPUT  MAJOR  AXIS  DIMENSION  (INCHES)' 

READ  (COM_UNIT,*)  DLE 

WRITE  (LOG_UNIT,*)  DLE  !  log  the  input 

FLOW-100. 

FHIGH-0. 

IBASE-0 
DLE-DLE*2 . 54 
TWOPI-ATAN!  1 .  )*8. 

WRITE! lOUT, * )  '***  TYPE  RETURN  IF  THERE  IS  NO  SUCH  FILE  ***' 

WRITE! lOUT, * )  'NOW  INPUT  THE  1-2G  FILE  NAME  ' 

READ  ! COM_UNIT, 7823 )  INFILE 

WRITE  (LOG  UNIT, 7823)  INFILE  1  log  the  input 

FORMAT!A20T 

IF! INFILE. EQ. '  ’)GO  TO  3350 


1 


non 


IF( FLOW.GT. 1 . ) FLOW-1 . 

IF( FHIGH.lt. 2. )FHIGH-2. 

CALL  TDATA( BA, INFILE) 

DO  3360  1-1,201 
CA( I )-BA( I ) 

3360  FA( 1 )-( 1 .+( I-l )/200 . )*TWOPI*DLE/3. 
lBASE-IBASE+200 

3350  WRITE( lOUT, * )  'INPUT  THE  2-4G  FILE  NAME' 

READ  ( COM_UNIT, 7823 )  INFILE 

WRITE  ( LOG_UNIT, 7823 )  INFILE  !  log  the  input 

IF( INFILE. EQ. '  ')GO  TO  3380 
IF( FLOW.GT. 2. ) FLOW-2. 

IF( FHIGH.LT. 4 . )FHIGH-4. 

CALL  TDATA( BA, INFILE) 

DO  9210  1-1,201 
CA( I+IBASE)-BA( I ) 

9210  FA( I  +  IBASE)-( 2 .  +  ( I-l ) /1 00. ) *TWOPl •DLE/3 . 

IBASE-IBASE+200 

3380  WRITE( lOUT, * )  'INPUT  THE  4-8G  FILE  NAME’ 

READ  (COM_UNIT,7823)  INFILE 

WRITE  ( LOG_UNIT, 7823 )  INFILE  !  log  the  input 

1F( INFILE. EQ. '  ')GO  TO  3390 
IF( FLOW.GT. 4 . ) FLOW-4 . 

I F( FHIGH.LT. 8. )FHIGH-8. 

CALL  TDATA( BA, INFILE) 

DO  9220  1-1,201 
CA( 1+1BASE)-BA( 1 ) 

9220  FA( l+IBASE)-( 4 .+( I-l )/50 . ) *TWOPl*DLE/3 . 

IBASE-IBASE+200 

3390  WRITE( lOUT, • )  'INPUT  THE  8-12G  FILE  NAME' 

READ  ( COM_UNIT, 7823 )  INFILE 

WRITE  ( LOG_UNIT, 7823 )  INFILE  !  log  the  input 

IF( INFILE. EQ. '  ' ) GO  TO  3400 

I F ( FLOW . GT . 8 . ) FLOW-8 . 

IF( FHIGH.LT. 12. )FHIGH-12. 

CALL  TDATA( BA, INFILE) 

C  NOTE  THAT  THIS  DO  LOOP  RUNS  201  TIMES 

DO  9230  1-1,201 
CA( I+IBASE)-BA( I ) 

9230  FA( I+IBASE)-(8.+(I-1)/50. ) *TW0PI*DLE/3 , 

3400  1-1 

IF( FLOW. GE.FHIGH) RETURN 
KLOW-INT(TWOPI*FLOW*DLE/3. )+l 
KHIGH-INT(TW0PI*FH1GH*DLE/3. ) 

DO  7720  K-KL0W,KHIGH 

7740  IF( K .GE . FA( I ) .AND. K. LE. FA{ I+l ) )GO  TO  7730 
I-I  +  l 

GO  TO  7740 

7730  IF( K .GT. 1024 )GO  TO  7777 

KT  CA(  I  +  l  )-CA<  I  )  ) /(  FA  (  T  *  1  ) -FA  (  T  '  )  *  (  K-FA  (  I  )  )+CA<  T  ' 

LA*  2  .,  SQP '■(  TC'jr  ( ,  .’.'‘I'LE) 

P( K, 1 )-10 . 24*( 10.+20. *ALOG10(CABS(KLA) ) ) 

7720  P( K+1024,1 )-1024 . *ATAN2 ( AIMAG( KLA) ,REAL(KLA) )/TWOPI 

7777  RETURN 

END 


SUBROUTINE  BUFSTR( IHDSZ , IDSSZ , ISI ZE , IB , HEAD , P , IDS ) 
INTEGER*2  HEAD(256,31) 

DIMENSION  P(  4096, 31  ) , IDS( 31  ) 

IC-IB+1 

DO  10  I-l, IHDSZ 
10  HEAD( I , IC)-HEAD( I ,1 ) 

DO  20  I-1,IS1ZE 
20  P( I , IC)-P( I , 1 ) 


6A 


w 

s 

10 

20 

% 

C 

C 

.N 

C 

IDSnC)-IDSl  1) 

RETURN 

END 


SUBROUTINE  BUFRD( IHDSZ , IDSSZ , I SI ZE , IB , HEAD , P , IDS ) 
INTEGER»2  HEAD(256,31) 

DIMENSION  P(4096,31) ,IDS(31) 

IC-IB+1 

DO  10  I-l, IHDSZ 
HEAD(I,1)-HEAD(1,IC) 

DO  20  I-1,ISIZE 
P(I,1)-P(I,1C) 

IDS( 1 )-lDS( IC) 

RETURN 

END 


SUBROUTINE  TDATA( CA, INFILE ) 

COMPLEX*8  CA(201) 

CHARACTER*20  INFILE 
INTEGER*2  BUFF(1024) 

REAL*4  AM( 201 ) ,PH(  201 ) 

EQUIVALENCE  (AM( 1 ' .BUFFI  53 ) ) , ( PH( 1 ) , BUFFI  4 55 ) ) 

INCLUDE  'SYS$LIBRARY:F0R10SDEF' 

OPEN I  UNI T- 1 9 , NAME- INFILE, STATUS- ' OLD ' , 

DEFAULTFILE-DEF, lOSTAT-IERR.ERR-aiOO ) 

DO  20  1-1,1024,256 
J-I+255 

READ I  19, 30) I  BUFFI K) ,K-I ,  J) 

FORMATI 256A2 ) 

DO  70  1-1,201 

AMI  I )-10**|AM| I )/20. ) 

CAIl)-CMPLX|AMi:)»COStPHU)),AH|l)*SlN(PH|I)  )  ) 

GO  TO  80 

1 F I  1 ERR . EQ . FORS I OS  FI LNOTFOU ) THEN 
WRITEITERM  UNIT,11T2)  INFILE 

FORMAT!'  FILE  :  ',A20,'  DOSE  NOT  EXIST',//,'  ENTER  FILENAME  AGAIN') 
ELSE  IF  |lERR.EQ.FORSIOS_FILNAHSPE)THEN 

WRITEITERHUNIT,*)  ' FILE INFILE, 'WAS  BAD,  ENTER  NEW  FILENAME' 
ELSE 

WRITEITERMUNIT,*)  'UNRECOVERABLE  ERROR  ,  CODE  -',IERR 

STOP 

ENDIF 

GO  TO  810 

CLOSEIUNIT-19,DISP-'SAVE' ) 

RETURN 

END 


CONTROLLER 

This  subroutine  controls  the  input  and  output 
devices  of  the  program. 

SUBROUTINE  CONTROLLER 

INCLUDE  'USER2 : | DURAL. CIMAG2 ) FTRN. FOR' 

INCLUDE  'USER2: (DURAL. CINAG2jMAGCMN2. FOR' 

Include  'USER2: [DURAL. CIHAG2jNAGCMN. FOR' 

INCLUDE  'USER2: [DURAL. Cl MAG2 I HEADER. CMN' 

INCLUDE  'USER2: (DURAL. CIHAG2jHSGBLK. FOR'  I  error  message  common  block 


-  •  -  n_jl  -  J."*  vV  s,  t 


I 

i 


► 

r 


i 

rx 


C 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


2 


4 


6 

8 


C 

C 

C 

c 

c 

c 


c 

c 

c 

c 


INTEGER* 4  FOR_RETCODE 
INTEGER*4  RETCODE 
INTEGER*4  SS$_NORMAL 
INTEGER*4 

LIB$GET_LUN  , 

LIB$FREE_LUN 

INTEGER*4  COH_STACK( 20) 

INTEGER* 4  TOP_COM  /  0  / 

CHARACTER* 40  DEF_STACK ( 20 ) 

INTEGER*4  TOP_DEF  /  0  / 

CHARACTER*70  NAME 

PARAMETER  ( SS$  NORMAL  -  ' 00000001 'X) 


I  fortran  return  code 
1  return  code 
I  "normal"  return  code 

1  'get  LU#'  RTL  routine 
1  'free  LU#'  RTL  routine 

I  the  command  units  stack 
I  top  of  the  command  stack 
1  default  stack 
I  top  of  default  stack 

t  name  of  the  file 


WRI  : 

ROUTINE  TO  WRITE  TO  DISK 

ENTRY  WRI (HEAD,  LINEl ,  LINE2,  LINES) 

WRITE( IOUT,2) 

FORMAT { '  FILE  NAME  ?  ',$) 

READ  (COM_UNIT,4)  FILNM 

WRITE  (LOG  UNIT, 4)  FILNM  1  log  the  input 

FORMAT(A50T 

OPEN(UNlT-l ,ERR-6, FILE-FILNM, TYPE- 'NEW' , 

+  FORM- ' UNFORMATTED ' ) 

GOTO  8 

CALL  LIB$SIGNAL(MAG_COM)  1  Print  an  error  message 

RETURN 

WRITEdOUT,*)  'DO  YOU  WANT  TO  CHANGE  THE  FILE  HEADER  Y-1 ,  N-0' 

READ  (COM_UNIT,*)  lYN 

WRITE  (LOG  UNIT,*)  lYN  1  log  the  input 

IF(IYN.EQ.T)  call  RLB  1  the  relabel  routine 

WRITE(l)  (HEAD( I ,1 ) , I-l ,256) 

WRITEd)  (P(I,1),I-1,ISIZE) 

CLOSE(UNIT-l) 

RETURN 


FILE: 

This  command  gives  control  of  the  program  to 
a  rnmmand  file.  The  output  however  will  still  qn 
wlie  level  I  OUT  is  set. 

ENTRY  FILE 


RETCODE-LIB$GET_LUN(  FILE_UNIT  )  I  get  the  LU#  for  the  command  file 
IF  (RETCODE  .NE.  SS$_NORMAL)  THEN  I  if  no  errors  occurred 

CALL  LIB$STOP(  %VAL ( RETCODE )  ) 

END  IF 

IF( PROC_FLAG)THEN 

NAME-'USER2 :[ DURAL . Cl MAG2 ) PROC2 . DAT' 

ELSE 

WRITEdOUT, 10)  I  prompt  for  file  name 

READ(COM_UNIT,15)  NAME  1  get  the  file  name 

END  IF 
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oo  s>  o  o  o  n  o  o  oooooonnno  o  ro  »-*»-*  n  o  o 

^  o  ai  o 


OPEN( 


I  open  the  command  file 


UNIT-FILE_UNIT  , 

+  FILE-NAME  , 

+  DEFAULTFILE-DEF  , 

+  IOSTAT-FOR_RETCODE  , 

♦  STATUS- ' OLD ' ) 

C 

IF  (FOR_RETCODE  . EQ .  29)  THEN 

CALL  LIBSSlGNAL{MAG_FILNOTFOU) 

RETURN 

ELSE  IF  (FOR_RETCODE  .NE.  0)  THEN 
CALL  LIB$SIGNAL(MAG_COM) 

RETURN 

END  IF 

C 

TOP_COM  -  TOP_COM  +  1  I  push  the  old  command  unit  onto  the  stack 

COM_STACK(  TOP_COM  )  -  COM_ONIT 
C 

TOP_DEF  -  TOP_DEF  +  1  1  push  Old  def  On  Stack 

DEF_STACK(  TOP_DEF  )  -  DEF 
C 

COM_UNIT  -  FILE_UNIT  t  give  control  to  the  file 


1  if  file  not  found 
1  print  error 

!  i f  not  normal 
1  Print  an  error  messag 


FORMATdX,  'COMMAND  FILE?',  $)  1  the  name  prompt 

FORMAT!  A  ) 

FORMATdX, 'ERROR  IN  RETRIEVING  LU# '  ) 

RETURN 


BACK: 

This  command  gives  control  bac)t  to 
the  unit  that  had  control  before  this 
one  took  over. 


ENTRY  BACK 


IF  (COMUNIT  .EQ.  TERMUNIT)  THEN 
RETURN 

END  IF 

F1LE_UNIT  -  COH_UNIT 
CLOSE!  UNIT  -  FILE_UNIT  ) 

COMUNIT  -  COM_STACK!  TOP_COM  ) 
TOPCOM  -  TOPCOM  -  1 

DEF  -  DEF_STACK!  TOP  DEF  ) 

TOP  DEF  -  TOP  DEF  -  I 


I  if  terminal  is  command  input 
1  then  ignore  this  command 

1  retrieve  file  LU* 

!  close  the  file 

J  pop  the  last  com_unit  off  the  stack 
!  pop  old  def 


RETCODE  -  LIB$FREE_LUN!  FILE_UNIT  )  I  give  LU#  back  to  system 

IF  !RETCODE  .NE.  SS$_NORMAL)  THEN  !  if  error 

CALL  LIB5ST0P(  %VAL(  RETCODE  )  )  I  then  Stop  and  give  reason 

END  IF 


IF  !COM_UNIT  .EQ.  TERMUNIT)  THEN  1  if  the  terminal  has  become 

WRITE! lOUT, 25)  I  the  command  unit 

END  IF  1  then  tell  the  user 

FORMAT!  IX,  'CONTROL  HAS  RETURNED  TO  THE  TERMINAL') 

RETURN 
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This  command  sets  the  default  for  all  the 
file  I/O. 


ENTRY  DEFFER 
WRITE( lOUT, 30  ) 

READ(COM_UNIT,35)  DEF 
WRITE  (L0G_UNIT,35)  DEF 
IF  (PROC_FLAG)  THEN 

DEF-'USER2: [ DURAL . CIHAG2 ) ' 
in  this  directory 
END  IF 

FORMATdX,  'DEFAULT?',  $) 

FORMAT( A40 ) 

RETURN 


1  ask  user  for  the  default 

1  log  the  input 
I  if  defining  a  procedure 
I  then  all  files  accessed  aie 


This  command  copies  all  the  user  input 
into  a  command  file. 


ENTRY  LOGGER 

RETCODE-LIB$GET_LUN(  L0G_UN1T  )  1  get  the  LUi  for  the  log  file 

IF  (PROC  FLAG)  THEN  1  if  defining  a  procedure 

NAME-'USER2: (DURAL. CIMAG2)PROC. DAT'  Ithls  is  the  definition  file 

ELSE 

WRITE( lOUT, 40 )  1  prompt  for  file  name 

READ{COM_UNIT,45)  NAME  !  get  the  file  name 

END  IF 

IF  (RETCODE  .EQ.  SS$_NORMAL)  THEN  I  if  no  errors  occurred 

OPEN(  UNIT-LOG_UNIT  ,  1  open  the  log  file 

FILE-NAME  , 

STATUS- 'UNKNOWN' ) 


END  IF 


WRITE( lOUT,  50 ) 

CALL  LIB$STOP( kVAL(RETCODE) ) 


FORMATdX,  'LOG  FILE?',  S) 

FORMAT!  A  ) 

FORMATdX, 'ERROR  IN  RETRIEVING  LUI') 


I  the  name  prompt 


STO  LOG: 


Stop  logging. 


ENTRY  STO_LOG 

CLOSE!  UNIT-LOGUNIT  )  I  close  the  logging  file 

RETCODE  -  LIB$FREE_LUN!  LOG  UNIT  )  I  give  LU#  back  to  system 

IF  ! RETCODE  .NE.  SS$_N0RMALT  THEN  I  if  error 

CALL  LIBSSTOP!  »VAL!  RETCODE  )  )1  then  Stop  and  give  reason 

END  IF 

L0G_UNIT  -  NULL_UNIT  1  log  to  null  device 

RETURN 


onoocjnoon  onnnno 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


C 


C 

C 


C 


SUBROUTINE  COMPLEX_COM: 

This  subroutine  contains  all  of  the  complex 
commands.  These  complex  commands  use  the  other 
command  routines  in  combinations  to  perform 
more  complex  operations. 


SUBROUTINE  COMPLEX_COM 

INCLUDE  ’USER2: [DURAL. Cl MAG2)FTRN. FOR' 

INCLUDE  'USER2: [ DURAL . C IMAG2 ) MAGCMN . FOR ' 

INCLUDE  'USER2: (DURAL. CIMAG21MAGCMN2. FOR'  I  controls  common  block 


REAL 

INTEGER* 4 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 


SEC 

INC 

GWID 

PROC_UNIT 
LIST_UNIT 
NEW_UNIT 
PROC2  UNIT 


I  number  of  nanoseconds  per  division 
!  increment 

1  width  of  a  grid  line 
!  LU#  for  defintion  file 
!  LUI  for  list  file 
I  LUt  for  list  of  output  file  names 
1  The  command  procedure  file 


INTEGER*4 

INTEGER*4 

INTEGER*4 

INTEGER*4 

INTEGER*4 


FORRETCODE 

PROCSTAT 

LISTSTAT 

NEWSTAT 

FOREOF  /  -1  / 


1  fortran  return  code 
1  status  of  a  read  from  def.  file 
I  status  of  a  read  from  list  file 
status  of  a  read  from  output  list  file 
!  end  of  file  code 


CHARACTER*80  OPER 

CHARACTER*!  REPLY 


input  data  from  the  definition  file 
1  answer  to  Y  or  N  question 


CHARACTER* 70  NAME 

CHARACTER* 70  NEWNAME 

CHARACTER* 7 5  COMMAND 


1  name  of  the  logging  file 
t  new  file  name 

I  command  sent  to  a  spawned  procedure 


CHARACTER* 70  PROCNAME 
CHARACTER* 70  LISTNAME 
CHARACTER*70  OUTNAME 


I  definition  file  name 
t  data  list  file 
1  output  list  name 


CLR_BUF: 

~  This  routine  clears  all  the  buffers 


ENTRY  CLRBUFCHEAD,  LINEl,  LINE2,  LINES) 
CALL  CLR 
DO  NB-1,40 

IF  (BUFFERS(NB,1)  .NE.  0)  THEN 

CALT,  SB  (HEAD,  LINEl,  LINE2, 

END  IF 

END  DO 
RETURN 


! 

1 

I 

t.INFI  ) 


clear  the  main  buffer 
do  for  all  the  buffers 
if  buffer  is  not  empty 
!  SP#  rant'  i  np 


GRID: 

This  makes  a  time  domain  file  with  a  harmonic 
value  of  100  at  every  given  number  of  nanoseconds 
in  the  main  buffer. 


ENTRY  GRID(HEAD,  LINEl,  LINE2,  LINES) 

HRITE( lOUT, 10 )  1  prompt  for  the  number  of  nanoseconds  between  lines 
READ( C0H_UN1T, 20 )  SEC  1  The  number  of  nanoseconds  between  lines 
CALL  CREATE(HEAD,  LINEl,  LINE2,  LINE3 )  1  Create  a  new  time  domain 
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C  tile 


10 
.•■  0 
10 


r 

c 


in  the  main  buffer 


TEMP-=100 . 0 
INC-SEC*''.!  !  The 
N-INC 

WRITE! lOUT, 30 ) 
READ(COM  UNIT,* 

DO  whileTn  .LE. 
DO 


!  The  amplitude  of  the  line  in  the  time  domain  file 
size  of  the  increment  (41  harmonics  -  1  nanosecond) 

1  The  current  harmonic  to  be  changed 
1  Prompt  for  the  width  of  the  line 
!  The  width  value 

For  all  the  harmonics  in  the  time  domain  file: 
1  The  line  will  have  a  width  of  GWID 

I  Change  the  value  of  these  harmonics 
IMake  the  line  wide  enough 


GWID 

4096)  1 

I-l ,GWID 

CALL  GRD 
N-N+1 

END  DO 

N-N-GWID  !  Figure 
N-N+INC 


out  the  increment  to  the  next  set  of  harmonics 


END  DO 

FORMAT!'  How  many  nanoseconds  per  division?') 

FORMAT! FlO . 2 ) 

FORMAT!'  How  many  harmonics  wide  should  the  line  be?') 
RETURN 


C  PROC: 

C  This  command  allows  the  user  to  specify  a  list 

C  of  files  to  be  processed  and  how  these  files  are  to 

C  be  processed.  This  facilitates  the  processing  of 

C  of  large  amounts  of  data. 


ENTRY  PROC 

PROC_FLAG- . TRUE .  !  tell  the  rest  of  the 

progiaro  that  Proc  is  being  rsed 

WRITE! lOUT, 31 )  I  does  the  user  have  a  def.  file  already? 

WRITE! lOUT, 32 ) 

READ!  COM_UNIT,  33  )  RET-LY  I  get  answei 

IF  (REPLY  .EQ.  'Y')  THEN  !  if  yes  then 

WRITE! lOUT, 34 )  I  ask  for  filename 

READ! COM_UNIT, 33 )  PROCNAME  !  read  in  file  name 

GOTO  58 


ELSE 

DEFINEFLAG-.TRUE.  1  set  definition  flag 

WRITE! lOUT, 35 )  1  instruction  prompt 

WRITE! lOUT, 40 ) 

WRITE! lOUT, 45 ) 
write! lOUT, 50  ) 

WRITE! lOUT, 52  ) 

WRITE! lOUT, 53 ) 

CALL  LOGGER  !  make  a  definition  file  of  the  commands 

PROCNAME- 'UEER2 :[ DURAL. CIMAG2 )PROC . DAT'  !  make  the 

definition  file 


}  ; 

f-nn  IF 
r  '  'P  (*  r,  r  ( 

'  Do  you  have 

a  procedure  definition  file') 

32 

FORMAT! 

'  for  this 

process  already?!?  or  N)') 

3  3 

FORMAT! 

A  ) 

34 

FORMAT! 

'  Fi lename  :  '  ) 

35 

FORMAT! 

'  Enter  the  p 

roccss  using  regular  commands  and') 

40 

FORMAT! 

'  NAHE.DAT 

for  a  filename.  For  the  new') 

45 

FORMAT! 

'  filename 

use  NEWNAHE.DAT.  When  finished') 

50 

FORMAT! 

'  defining 

the  process  use  the  command  DONE.’) 

52 

FORMAT! 

'  (Warning: 

a  filename  must  be  listed  for  ’) 

53 

FORMAT! 

each  time  it  is  used.)') 

c 

RETURN 

ENTRY  PROC2 

!  reentry  after  DONE  command 

CALL  STO  LOG 

!  close  up  the  definition  file 

DEFINE 

FLAG- . FALSE. 

!  reset  definition  flag 

WRITE! lOUT, 55 ) 

;  ask  whether  he  wants  to  save  def. 

/o 


O  O  'J 


READ! COM_UNIT, 56 )  REPLY  !  read  in  answer 

IF(REPLY  .EQ.  'Y')  THEN  1  if  he  wants  to  save  the  de£ 

WRITE( lOUT, 57 )  1  prompt  for  a  file  name 

READ(COM  UNIT, 56)  NAME  1  read  the  filename 

COMHANDd : 33 )-' SCOPY  USER2: I DURAL.CIMAG2JPR0C.DAT  '!  the  first 
half  of  the  command 

COMMAND! 34 : 75 )-NAME  1  the  filename  that  the  copy  is  output  to 
CALL  LIB$SPAWN{  COMMAND  )  I  Spawn  a  process  to  copy  the  file 

END  IF 

FORMAT!'  Do  you  wish  to  save  this  procedure  definition?! Y  or  N)') 

FORMAT!  A  ) 

FORMAT!'  Filename:') 


WRITE! lOUT, 59 ) 

READ! COMUNIT, 75 )  REPLY 
IF!REPLY  .EQ.  'Y')  THEN 
WRITE! lOUT, 60 ) 

READ!COM  UNIT, 75)  LISTNAME 
CALL  OPEfJER  I  LI  STUNIT,  LISTNAME) 
CLOSE!LIST  UNIT) 


1  As):  if  there  is  already  a  data  list 
I  get  answer 


1  data  list  filename 
!  get  a  LUtt 

I  close  the  unit  back  up 


ELSE 


open  a  file 
prompt 


I  write  name  into  list  file 

i  read  next  name 


L I STNAME- ' USER2 ; ! DURAL . CIMAG2 ) LI  ST . DAT ' 

CALL  OPENER! LI STUNIT, LISTNAME) 

WRITE! IOUT,61 ) 

WRITE!I0UT,65) 

WRITE! lOUT, 70) 

READICOM  UNIT, 75)  NAME  I  read  first  name  in  list 

DO  WHILeTnamE  .NE.  'DONE')  I  when  finished  use 

C  the  word  DONE 

WRITE!LIST  UNIT, 75)  NAM 
read!com_uRit,75)  name 

END  DO 

CLOSE! LIST_UNIT)  I  close  the  list  file 

WRITE! lOUT, 71 )  I  ask  the  user  if  it  should  be  saved 

READ(C0M_UN1T,75)  REPLY 
IF! REPLY  .EQ.  'Y')  THEN 

WRITE! lOUT, 60)  lask  for  file  name 

COMMAND! 1: 34 )-'$COPY  USER2 :[ DURAL . CIMAG2 ) LIST . DAT  ' 

C  the  first  half  of  the  command 

COMMAND! 35: 75) -NAME  )  the  filename  that  the  copy 

C  is  output  to 

CALL  L1B$SPAWN1  COMMAND  )  I  spawn  a  process  to 

C  copy  the  file 

END  IF 

END  IF 

FORMAT!'  Do  you  have  a  data  list  file?!Y  or  N)') 

FORMAT!'  Filename:') 

FORMAT!'  Enter  the  list  of  data  files,  following') 

FORMAT!'  each  with  <CR>.  When  finished  type  the') 

FORMAT! '  word  DONE. ' ) 

FORMAT!'  Do  you  wish  to  save  this  data  list7!Y  ot  tJ'') 

FORMAT!  A  ) 


WRITE! lOUT, 76 )  !  is  there  an  output  name  file 

READ!COM_UNIT,75)  REPLY  I  read  answer 

IF  ! REPLY  .EQ.  'Y')  THEN 

WRITE! lOUT, 77 )  1  prompt  for  a  filename 

READ!COM_UNIT,75)  OUTNAME 

CALL  OPENER !NEW_UNIT, OUTNAME)  !  get  a  LU# 

CLOSE! NEW_UNIT) 

ELSE 

0UTNAME-'0SEB2:(CIMAG21NEWLIST.DAT'  1  list  of  new  file  names 

CALL  OPENER! NEWUN IT, OUTNAME)  I  open  newlist  file 

WRITE! lOUT, 80 )  1  prompt 

WRITE!IOUT,85) 


on  n  VO  n  n 


I  initial  read 

1  put  the  new  name  in  file 
1  read  next  file  name 


I  save? 


WRITE) IOUT,90) 

READ(COM  UNIT, 95)  NAME 
DO  WHILETnAME  .NE.  'DONE') 

WRITE) NEW  UNIT, 95)  NAME 
READ)COM_UNIT,95)  NAME 

END  DO 

CLOSE )NEW_UNIT) 

WRITE) lOUT, 78 ) 

READ)COM_UNIT,75)  REPLY 
IF  )REPLY  .EQ.  'Y')  THEN 
WRITE) IOUT,77 ) 

READ) COM  UNIT, 75)  NAME 

C0MMAND)T:37)-'$C0PY  USER2 ;[ DURAL . CIMAG2 ) NEWLI ST . DAT  ’ 
C  the  first  half  of  the  command 

COMMAND) 38: 75) -NAME  I  the  filename 

C  that  the  copy  is  output  to 

CALL  LIB$SPAWN)  COMMAND  )  1  spawn  a  process 

C  to  copy  the  file 

END  IF 


END  IF 

76  FORMAT)'  Is  there  an  output  filename  list?)Y  or  N)') 

77  FORMAT)'  Filename:') 

78  FORMAT)'  DO  you  wish  to  save  this  li6t?)Y  or  N)') 

80  FORMAT) '  Enter  a  list  of  the  output  file  names  in' )  1  The  prompt 

85  FORMAT)'  the  order  they  are  to  be  used.  When') 

90  FORMAT)'  finished  type  DONE.') 

95  FORMAT)  A  ) 


WRITE) lOUT, 97) 

7  FORMAT)'  Your  data  is  being  processed.') 

CALL  OPENER)PROC_UNIT,PROCNAME) 

OPEN)  UNIT-LIST_UNIT  , 

♦  FILE-LISTNAME  , 

+  STATUS- ' OLD ' ) 

OPEN)  UNIT-NEW  UNIT  , 

+  FILE-OUTNAME  , 

+  STATUS- ' OLD ' ) 

NAME-'USER2!(DURAL. CIMAG2JPR0C2.DAT'  I  The  command  procedure 

CALL  OPENER) PROC2  UNIT, NAME) 


C 

c 

c 

c 

c 


READ)  LISTUNIT,  100,  lOSTAT-LISTSTAT  )  NAME  I  Initial  read  from 
data  file 

DO  WHILE)LISTSTAT  .NE.  FOREOF )  1  do  while  there  is  still  data 

to  be  processed 

READ)  PROCUNIT  ,  100  ,  lOSTAT-PROCSTAT  )  OPER  1  initial  read 
from  definition 

DO  WHILE) PROCSTAT  .NE.  FOR_EOF)  I  do  entire  definition 

IF)  OPER  .EQ.  'NAME.DAT'  )  THEN  1  if  the  data 

file  name  is  needed  then  get  the  name 

IF)LISTSTAT  .EQ.  FOREOF )  THEN  !  if  there  is 

no  filename  print  error 

CALL  LIB$SIGNAL)HAG  COM) 


Cl  read  next  data  entry 
C  file  is  needed 
Cl  read  output  file  name 
C  print  message 


END  IF 
OPER-NAME 

READ)  LIST  UNIT,  100,  lOSTAT-LISTSTAT  )  NAME 


ELSE  IF)  OPER  .EQ.  'NEWNAHE.DAT'  )  THEN  1  if  output 
READ)  NEWUNIT,  100  ,  lOSTAT-NEWSTAT  )  NEWNAME 
IF)NEWSTAT  .EQ.  FOREOF )  THEN  1  if  error 
CALL  LIB$SIGNAL)MAG  CON) 
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C  new  name 


END  IF 

OPER>NEWNAHE 


1  give  command  procedure  the 


i 


9 


5 


END  IF 

WRITE( PROC2_UNIT, 100 )  OPER 

READ(  PR0C_UNIT,  100,  lOSTAT-PROCSTAT )  OPEN 

C  next  def.  statement 
END  DO 

CLOSE(PROC  UNIT) 


read 


100 

C 


OPEN(  UNIT-PROC_UNIT  , 

+  FILE-PROCNAME  , 

+  STATUS-'OLD' ) 

END  DO 
FORMAT (  A  ) 

CLOSE( PROCUNIT) 

CLOSE( LIST_UNIT) 

CLOSE (NEW_UN1T) 

CLOSE( PROC2_UNIT) 

CALL  LIB$FRE'_LUN(  PROC_UNIT  ) 
CALL  LIB$FREE_LUN(  LIST_UNIT  ) 
CALL  LIBSrREE_LUN(  NEWUNIT  ) 
CALL  LIB$FREE_LUN(  PROC2_UNIT  ) 

CALL  FILE 
RETURN 


!  Start  definition  over 
I  start  definition  over 


1  close  all  the  files 


1  f  ree  up  the  LU(t '  s 


1  execute  the  command  procedure 


II 

I 


W 

5; 


ENTRY  PROC3 

IF  (PROCNAME  .EQ.  ' USER2 :( DURAL . CIMAG2 ] PROC . DAT ' )  THEN  !  delete  all 
C  the  utility  files 

OPEN ( UNI T-1, FILE- 'PROC. DAT' .TYPE- 'OLD' , 

+  DEFAULTFILE-'USER2; IDURAL.CIMAG2) ' ,DISP- 'DELETE' ) 

CLOSE(UNIT-l) 

END  IF 

IF  (LISTNAME  .EQ.  ' USER2 : t DURAL. CIMAG2 ) LIST. DAT' )  THEN 
0PEN(UNIT-1,FILE-'L1ST.DAT' ,TYPE-'0LD' , 

DEFAULTFILE-  'OSER2  :  (  DURAL .  CIMAG2  ]  '  ,  DISP-  '  DELETE  '  ) 
CLOSE(UNIT-l ) 

END  IF 

IF  (OUTNAME  .EQ.  ' USER2 :{ DURAL . CIMAG2 JNEWLIST . DAT ' )  THEN 
OPEN(UNIT-l,FILE-'NEWLIST.DAT' ,TYPE-'OLD' , 

+  DEFAULTFILE-'USER2: [DURAL. C1HAG2 ) ' , DISP- ' DELETE '  ) 

CLOSE(UNIT-l ) 

END  IF 

OPEN ( UNI T-1 , FILE- 'PROC2.DAT' , TYPE- 'OLD' , 

+  OEFAULTFILE-'USER2:(DURAL.CIMAG2J' , DISP- ' DELETE ' ) 

CLOSE( UNIT-1 ) 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


PROC  FLAG-. FALSE. 

RETURN 

END 

SUBROUTINE  OPENER ( LUNIT, NAME ) 

This  routine  opens  up  a  file  and  gets  a  LUI 


turn  off  flaa 


1NTEGER*4  FOR_RETCODE 
INTEGER»4  RETCODE 
INTEGER*4  SS$  NORMAL 


1  fortran  return  code 
!  return  code  for  system  operations 

!  normal  return 


INTEGER*4  L_UNIT 
CHARACTER* 70  NAME 

PARAMETER  ( SS$  NORMAL  -  ' 00000001 'X) 


RETCODE-LIB$GET_LUN{  L  UNIT  ) 

IF  (RETCODE  .EQ.  SS$_N(5rMAL)  THEN 
OPEN(  UNIT-L_UNIT  , 

+  FILE-NAME  , 

+  STATUS- • UNKNOWN ' ) 


WRITE! IOUT,55) 

CALL  LIB$STOP( %VAL( RETCODE) ) 


FORMAT! IX, 'ERROR  IN  RETRIEVING  LU#') 

RETURN 

END 


!  logical  unit  # 
!  file  name 
1  normal  return  code 


1  get  the  LU#  for  the  list  file 
1  if  no  errors  occurred 
1  open  the  list  file 


1  error  conditions 


SUBROUTINE  OUTSIDER: 

These  are  routines  are  more  basic  commands 
such  as  those  found  in  COM  FILE. 


SUBROUTINE  OUTSIDER 

I NCLUDE  ' USER2 : t  DURAL . C 1 MAG2 ) FTRN . FOR ' 
INCLUDE  ’USER2: (DURAL. CIMAG2)NAGCMN. FOR’ 
INCLUDE  'USER2: ( DURAL . CIMAG2 )MAGCMN2 . FOR' 
INCLUDE  'USER2 : ( DURAL . CINAG2 ] HEADER . CNN’ 


CREATE: 

Creates  a  blank  time  domain  file  with  a 
header  and  stores  it  in  the  main  buffer. 


ENTRY  CREATE(HEAD,  LINEl,  LINE2,  LINED 

WRITE! lOUT, • )  'Type  in  the  header :( three  lines)' 
READ! COMUNIT, 10 )  LINEl 
no*.^  h»ad»r 

pr.MM  r-Mi  UMIT,  10)  LINE2 
READ!COM  UNIT, 10)  LINE3 
FORMAT! eCAl ) 

DO  1-1,ISI2E 

P!I,l)-0  !  The  ori 

END  DO 
IDS! 1 )-0 

DOMAIN-'TIME'  !  This  file  is  Ir 

RETURN 


!  prompt 
!  read  in 


!  The  original  array  is  'O’s 
!  This  file  is  in  the  time  domain 


CHANGE : 


This  routine  changes  a  given  harmonic  in  a  data 
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ENTRY  CHANGE 

WRI TE ( lOUT , 30 )  1  Prompt  for  the  harmonic  to  be  changed 

READ( COM_UNIT, * )  N  t  Read  in  the  value  for  the  harmonic 

WRITE( lOUT, 40 )  (N,P(N,1))  I  Give  the  current  value  of  the  harmonic 

WRITE( lOUT, 50 )  I  Prompt  for  the  new  value 

READ(COM_UNIT,60 )  TEMP  I  The  new  value 

ENTRY  GRD 

P(N,1)-TEMP  1  Put  the  new  value  in  the  specified  harmonic 

FORMAT(X, 'Which  harmonic  do  you  wish  to  change?') 

FORMAT ( '  Current  Value:  ',I5,F10.2) 

FORMAT ( '  New  Value:  ') 

FORMAT (FI  0.2) 

RETURN 


(;**********************************«******************************** 

C  FTRAN  READ  COMMANDS  ( BY  A .  DOMINEK) 

(^************************4******************************************** 


THE  CALLING  PROGRAM  MUST  CONTAIN  THE  FOLLOWING 
COMMON  BUFFER, NDIM,ANST,AINC 
DIMENSION  AMPL( 5000 ) ,PHS( 5000) 

BYTE  BUFFER( 35000 ) 

INTEGER*2  INPFILE(15) 

USER  SUPPLIES  THE  FILE  NAME  IN  VARIABLE  'INPFILE' 

IN  FREQUENCY  DOMAIN 

AMPL  CONTAINS  AMPLITUDE  OF  DATA  IN  DB 

PHS  CONTAINS  PHASE  OF  DATA  IN  DEGREES 

NDIM  13  THE  NUMBER  OF  FREQUENCY  SAMPLES 

ANST  IS  THE  FREQUENCY  (MHZ)  OF  THE  FIRST  SAMPLE 

AINC  IS  THE  DELTA  FREQUENCY  (MHZ) 

IN  TIME  DOMAIN 

AMPL  CONTAINS  THE  AMPLITUDE  OF  TIME  WAVEFORM 
PHS  CONTAINS  ZERO 

NDIM  IS  THE  NUMBER  OF  TIME  SAMPLES  -  4096 

ANST  IS  THE  STARTING  TIME  -l/DF/2  ( DF-DELTA  FREQUENCY) 

AINC  IS  THE  DELTA  TIME  *  1.E5  -l/DF/4096*! . E5 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


SUBROUTINE  FREA( INPFILE , AMPL, PHS ) 


PROGRAM  NAME  : USERl : ( DOMI ) REAV . FOR 

THIS  PROGRAM  BEADS  BACKSCATTER  DATA  FILES  STORED 

OM  VAX  DISKS.  WITH  11/23  FORMAT 


INCLUDE 'USER2 : ( DURAL. CIMAG2 IMAGCMN.FOR' 
INCLUDE  'USER2: ( DURAL . CIMAG2 )MAGCMN2 . FOR ' 
INCLUDE  'USER2; (DURAL. CIMAG21FTRN. FOR' 


INTEGER*2  LIN1(30) , LIN2 ( 30 ) , PARAH( 30 ) 
REAL*4  AP(IOOOO) 


DEFINE  BUFFER  STRUCTURE 


EQUIVALENCE( LINl ( 1 ) , BUFFER ( 1 ) ) , ( LIN2( 1 ) ,BUFFER(61 ) ) 
1 , ( PARAM( 1 ) ,BUFFER( 121 ) ) , (AP( 1 ) ,BUFFER( 361 ) ) 

EQUI VALENCE ( LINEK 1 ) ,LINl( 1 ) ) , ( LINE2 ( 1 ) , LIN2 ( 1 ) ) 

1 , ( LINE3( 1 ) ,PARAM( 1 ) ) 


READ  A  FILE 


•m 


CALL  TTR( INPFILE) 
TYPE  105,LIN1 
TYPE  105,LIN2 
TYPE  lOS.PARAM 
FORMAT(X, 30A2 ) 


GET  NUMERICAL  INFORMATION  FROM  THE  THIRD  LINE 
OF  THE  HEADER 


CALL  DDCDE 


DIVIDE  AN  AMP-PHASE  ARRAY  INTO 
AN  AMP  ARRAY  AND  A  PHASE  ARRAY 

DO  199  NN-1,NDIM 
AMPL ( NN ) -AP ( 2  *NN- 1 ) 

IF(AMPL(NN) .GT.35)AMPL(NN)-35. 

PHS(NN)-AP(2*NN) 

CARY  INFORMATION  TO  'CIMAG'  (TIME  DOMAIN  ONLY) 
DO  1  I-1,NDIM 
P( I , 1 ) -AMPL ( I ) 

CHECK  FOR  BAD  DATA  POINTS,  I.E.,  AMPL ( I ) . GT . 99 5 

CALL  ERRF(AMPL,PHS) 

RETURN 

END 


SUBROUTINE  REU ( INPFI LE , AMPL , PHS ) 

THIS  PROGRAM  READS  BACKSCATTER  DATA  FILES  STORED 
ON  VAX  DISKS.  WITH  750  FORMAT 


INCLUDE ' USER2 : ( DURAL . CIMAG2  JMAGCMN . FOR' 

INCLUDE  'USER2; ( DURAL . CIMAG2 )MAGCMN2 . FOR’ 

INCLUDE  ’USER2: [ DURAL . CIMAG2 1 FTRN . FOR ’ 

BYTE  PARAM(60) 

INTEGER*2  Ll ( 30 ) , L2 ( 30 ) , L3 ( 30 ) 

■aL*4  AP(IOOOO) 

EwUIVALENCE  ( L 3 ( 1 ) , PARAM ( 1 ) ) 

KOinVALFNCE( L) ( 1 ) .BUFFER! 1 ) ) , ( L2 ( 1 ) , BUFFER! 61 ) ) 

1, (L3! 1 ) .BUFFER! 121 ) ) , ( AP ( 1 ), BUFFER ( 361 ) ) 

WRITE! lOUT, 2 ) 

FORMAT! ' S ENTER  DATA  FILE  NAME;  ') 

READ(COM  UNIT, 10)  INPFILE 
FORMAT!  1'5a2  ) 

WRITEILOG  UNIT, 10 ) INPFILE 
INPFILEd?  )-0 

OPEN (UNIT-8,  NAME- INPFILE, TYPE- 'OLD' , FORM- ' UNFORMATTED' , 
1READ0NLY,ERR-1 ) 

READ(8)  Ll 
READ(8)  L2 
READ! 8)  L3 
TYPE  100, Ll 
TYPE  100, L2 
TYPE  100, L3 
FORMAT(X, 30A2) 


DECODE ( 4,4, PARAM( 3 ) , ERR-1 001 ) NDIM 
DECODE ( 5, 5, PARAM( 11 ) , ERR-1002 ) lANST 
ANST.FLOAT( lANST) 

DECODE (5,5, PARAM ( 20 ) , ERR-1 003 ) lAINC 
AINC-FLOAT( lAINC ) 


4 

FORMAT! 14 ) 

• 

5 

FORMAT! 15) 

9 

DO  200  I-l, NDIM 
READ!  8)  AMPLd) 
10-2*1-1 

IE-2*I 

AP( IO)-AMPL( 1 ) 

,  PHS!  I  ) 

V. 

200 

AP! IE)-PHS! I  ) 

CLOSE  (UNIT-8, DISP-'SAVE' ) 

V. 

M  > 

CALL  ERRF!AMPL, 
GO  TO  99 

PHS, NDIM) 

•»  1 
*'J' 

1001 

WRITE! lOUT ,* )  ' 

DECODE  ERROR 

NDIM' 

*'1 

1002 

GO  TO  99 

WRITE! lOUT, * )  • 

DECODE  ERROR 

ANST' 

r-^ 

GO  TO  99 

1003 

WRITE! lOUT, • )  ' 

DECODE  ERROR 

AINC' 

99 

RETURN 

SUBROUTINE  TTR(INPFILE) 

BYTE  TBUFF( 512 ) 

INCLUDE  'SYS$LIBRARy:FORIOSDEF’ 

I NCLUDE ' US  ER2 : ( DURAL . Cl NAG2 1 KAGCMN . FOR ' 
INCLUDE  'USER2: [ DURAL . CIMAG2 J MAGCMN2 . FOR ' 
INCLUDE  'USER2: ( DURAL . CIMAG2 ) FTRN . FOR ’ 


WRITE! IOUT,llll ) 

FORMAT! IX,'  ENTER  DATA  FILE  NAME:') 
READ<COM  UNIT, 2222)  INPFILE 
FORMAT!  1'5a2  ) 

WRITE!LOG  UNIT, 2222)  INPFILE 

INPFILE! l5)-0 

IB-1 

ICNT-0 

OPEN!UNIT-8 .NAME- INPFILE, READONLY, TYPE- 


'OLD' , lOSTAT-IERR, ERR-BIOO ) 


SET  BLOCK  LENGTH  IN  BYTES 

IF! IB.EQ. 1 )LEN-512-9*4 
IF! IB.GT.l )LEN-512-26*4 

BEAD  A  BLOCK  OF  512  BYTES 

PFAP! R , PO , END-90 )  TBUFF 
FORMAT! 512A1 ) 

STORE  A  BLOCK  INTO  THE  BUFFER  ACCORDING  TO  ITS  LENGTH 

DO  85  I-1,LEN 

BUFFER! ICNT+I )-TBUFF!  I  ) 

lB-IB+1 

ICNT-ICNT+LEN 

GO  TO  82 

DO  86  I-l , LEN 

BUFFER! ICNT+I )-TBUFF! I ) 

ELIMINATE  BLANK  SPACES  IN  BETWEEN  EACH  CHARACTER 
IN  A  FILE  HEADER 


DO  40  1-1,180 

40  BUFFER( I )-BUFrER( 2*1-1 ) 

GO  TO  331 

8100  IF( IERR.EQ.FOR$IOS_FILNOTFOU)THEN 

WRITE( lOUT, 1112 )  INPFILE 

1112  FORMAT ( '  FILE  ',15A2,'WAS  NOT  FOUND' S' , 

1  'ENTER  FILENAME  AGAIN:  ') 

ELSE  IF  ( IERR.EQ.FORSlOS_FILNAMSPE)THEN 
WRITE( 6 , 1113 )  INPFILE 

1113  FORMAT( 'S' , 'FILE  :  ' , 1 5A2 , '  WAS  BAD,  ENTER  NEW  FILENAME 
ELSE 

TYPE  *, 'UNRECOVERABLE  ERROR,  CODE-',IERR 
STOP 
ENDIF 
GO  TO  6 

331  CLOSE(UNlT-8,DISP-'SAVE' ) 

RETURN 
END 

! 

SUBROUTINE  DDCDE 

'  INTEGER*4  IMIN , I INC , NDIM 

'  INCLUDE'USER2 : ( DURAL . Cl MAG2 ] MAGCMN . FOR ' 

!  INCLUDE  'USER2! (DURAL. CIMAG21MAGCMN2. FOR' 

,  INCLUDE  'USER2: [DURAL. CIMAG2]rTRN. FOR' 

i  C 

C  NO  OF  DATA  POINTS  IS  STORED  IN  FOUR  CHARACTERS,  AND 

'  C  STARTING  ANGLE  AND  ANGLE  INC.  IN  5  CHARACTERS 

C 

CHARACTER*!  CNLl 
CHARACTER* 4  CNL 
CHARACTER*5  CFF,CINC 
CHARACTER*!  ECH,TCAS 
DATA  ECH,ZERO/'-' , '0'/ 

!  EQUIVALENCE  ( BUFFER( 1 2 3 ) , CNL ) , ( BUFFER ( 1 31 ) ,CFF) ,  (  BUFrER(  140 )  ,  CINC ) 

'  EQUIVALENCE  ( BUFFER( 12 3 ) ,TCAS ),( BUFFER! 124 ), CNLl ) 

C 

C 

C  CONVERT  CHARACTERS  INTO  THEIR  NUMERICAL  EQUIVALENTS 

.  C 

1  IF{ECH.EQ.TCAS)  THEN 

decode! 3 , 102 , CNLl , ERR-9 )NDIM 
I  ELSE 

I  DECODE! 4, 101, CNL, ERR-9)NDIH 

'  END  IF 

100  rORMAT!l5) 

1  101  FORMAT! 14) 

102  F0RMAT!I3) 

103  prf’nnF!5,10n,CFF,ERR-'><»)IMIN 

104  DECODE! 5, loo, CINC, ERR-999)IINC 
ANST-FLOAT! IMIN) 

AINC-FLOAT!  UNO 

RETURN 

C 

9  WRITE  !6,200)  'NDIM' 

READ  !5,*,ERR-91)  NDIM 
GO  TO  103 

C 

99  WRITE  (6,200,ERR-991 )  'ANST' 

READ! 5,*)  ANST 
GO  TO  104 
C 

999  WRITE  (6,200,ERR-9991)  'AINC' 

READ! 5,*)  AINC 
RETURN 


78 


c 

200  FORMAT  'HAVING  PROBLEMS  READING  HEADER.  ENTER  ',A4, 

1'  MANUALLY:  ') 

C 

91  WRITE(6,300) 

GO  TO  9 
C 

991  WRITE(6,300) 

GO  TO  99 
C 

9991  WRITE(6,300) 

GO  TO  999 
C 

300  FORMAT( lx, '  “‘INVALID  ENTRY***') 

C 

END 

c  . * . 

SUBROUTINE  ERRF ( AMPL , PHS ) 

COMPLEX  Cl,C2,CD 

INCLUDE' USER2 : ( DURAL . CIMAG2 IKAGCMN . FOR' 

INCLUDE  'USER2! ( DURAL . Cl MAG2 JMAGCMN2 . FOR ' 

INCLUDE  'USER2: [ DURAL . CIMAG2 ) FTRN . FOR ' 

DO  1  I-1,NDIM 

IF(PHS(I) .GT.995. )  WRITE(6,2)  I . AMPL( I ) , PHS ( I ) 

1  CONTINUE 

2  FORMAT! IX, 16HERROB  AT  DATA  PT, 1 1 4 , 4HMAG- , 1 FI 0 . 4 , 4HPHS- , 1  FI 0 . 4 ) 

C  CHECK  LEFT  HAND  END  POINT 

IF( AMPL( 1 ) .GT. 100 . )  THEN 
DO  200  I-2,NDIM 

IF(AMPL( I ) .LE. 100.  .AND.  AHPL( 1*1 ) . LE . 1 00 . )  THEN 
Al-10.**(AMPL(I)/20. ) 

Cl-CMPLX(Al*COSD(PHS( I ) ) ,A1*SIND( PHS( I ) ) ) 

A2-10. **( AMPL ( 1+1 )/20. ) 

C2-CMPLX(A2*COSD(PHS( I+l) ) , A2*SIND( PHS ( 1+1 ) ) ) 

CD-C1-C2 

RD-REAL(CD) 

AD-AIHAG(CO) 

DO  212  11-1,1-1 
RC-REAL(C1 )+RD*II 
AC-AIMAG(C1)+AD*1I 

AHPL( II )-20. *LOG10(SQRT(RC*RC+AC*AC) ) 

PHS( II )-ATAN2D(AC,RC) 

212  CONTINUE 

GO  TO  211 
ELSE 
END  IF 

200  CONTINUE 

ELSE 
END  IF 

C  CHECK  RIGHT  HAND  END  POINT 

211  IF(AMPL(NDIM) .GT.IOO. )  THEN 

DO  220  I-1,NDIH 
J-NDIM-I 

IF(AMPL( J) .LE.IOO.  .AND.  AMPL( J-1 ) . LE . 100 . )  THEN 
Al-10.**(AMPL( J)/20. ) 

Cl-CMPLX(Al*COSD(PHS( J) ) ,A1*SIND( PHS ( J ) ) ) 

A2-10.**(AMPL( J-1 )/20. ) 

C2-CHPLX(A2*COSD(PHS( J-1 ) ) , A2*SIND( PHS ( J-1 ) ) ) 

CD-C1-C2 

RD-REAL(CD) 

AD-AIMAG(CD) 

DO  222  II-J+1,NDIM 
RC-REAL(C1 )+RD* ( II-J) 

AC-AIMAG( Cl ) +AD* (II-J) 

AMPL( II )-20.*LOG10(SQRT(RC*RC+AC*AC) ) 

PHS{ II )-ATAN2D(AC,RC) 
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222  CONTINUE 

GO  TO  221 
ELSE 
END  IF 

220  CONTINUE 
ELSE 
END  IF 

C  CHECK  INTERIOR  POINTS 

221  DO  230  1-2,NDIM-1 
IF(AMPL( I ) .GT.IOO. )  THEN 

DO  240  K>I-t-l,NDIM 
IF(AMPL(K) .LE.IOO. )  THEN 
Al-10.**(AMPL(I-l)/20.) 

Cl-CHPLX(Al*COSD(PHS( 1-1 ) ) ,Al*SIND( PHS( I-l ) ) ) 
A2-10.**(AMPL(K)/20. ) 

C2-CHPLX(A2*COSD(PHS(X) ) ,A2*SIND( PHS ( K ) ) ) 
CD-(C1-C2)/(K-I+1) 

RO-REAL(CD) 

AD-AIHAG(CO) 

DO  241  II-I,K-1 
RC-REAL ( Cl ) +RD* (11-1+1) 

AC-A1MAG(C1 )+AD*( ll-I+l) 
AMPL(11)-20.*LOG10(SQRT(RC*RC+AC*AC) ) 

PHS( 11 )-ATAN2D(AC,RC) 

241  CONTINUE 

GO  TO  230 
ELSE 
END  IF 


240 

CONTINUE 

ELSE 

END  IF 

230 

CONTINUE 

RETURN 

END 


noon 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  MAGCMN.FOR 

cc  PROGRAM  COMMON  BLOCKS 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 


INTEGER*2 

INTEGER*2 

INTEGER*4 

INTEGER*4 

REAL 

LOGICAL 

LOGICAL 

CHARACTER* 4 

CHARACTER* 50 

CHARACTER*60 

CHARACTER* 40 

DIMENSION 

DIMENSION 

DIMENSION 


HEAD(256,31) ,LINE1(30) ,LINE2(30) ,LINE3( 30) 
BUFFERS(400,30) 


NB,  ISIZE 
PNUM , ELEMENT, ELEM 
VALUE , PVALUE , WEIGHT 

EFLAG  terror  flag 

ECHO 

DOMAIN,  DOM(40) 

I NAME , JNAME , F I LNH 
TITLE, XLAB,YLAB 
FNAME 

P(4096,31),IDS(31),S(1024),Q(4096) ,ARRAy(3,100,100) 
NUMB(3) ,ANG( 3,400 ),BUFN( 3,400 ),CNTR( 3, 400) 
CARRAY{4,8) ,CLRTAB(6),ARRAY2( 3,100,100) 

COMPLEX  A(4096) ,CCI,CCX,BCX,DCX 

COMPLEX  BA(201 ) , KLA, CA( 2049 ) ,ASQ,A5S1 ,ASS2 ,ASDl , ASC , ASD , ASA, ASB 
CHARACTER*20  INFILE 
REAL  FA(2049) 

BYTE  MACRO(128) 

INTEGER*2  BFILE(6) 


COMMON/BLKl/INAME,  JNAME,  I  The  main  common  bloctc 

FILNH,  TITLE,  XLAB,  YLAB,  P,  IDS,  S,  Q, 

+  ARRAY,  NUMB,  ANG,  BUFN,  CNTR,  CARRAY,  CLRTAB, 

+  A,  CCI,  CCX,  BCX,  DCX,  BA,  KLA,  CA,  ASQ,  ASSl, 

*■  ASS2,  ASDl,  ASC,  ASD,  ASA,  ASB,  INFILE,  FA, 

*  MACRO,  BFILE,  BUFFERS,  DOMAIN,  DOM,  NB,  ECHO, 

>  ARRAY2,  ISIZE, FNAME 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  MAGCMN2 . FOR 

cc  THIS  IS  THE  PROGRAM  CONTROLS  COMMON  BLOCK 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

INTEGER*^ 

+  COM_UNIT  ,  t  LUt  for  the  command  input 

■f  lOUT  ,  t  LUt  for  the  program  output 

+  TERM_UNIT  ,  i  LUt  for  the  terminal 

+  FILE_UNIT  ,  i  LUt  for  the  command  file 

■f  L0G_UNIT  ,  t  LUt  for  the  log  file 

-f  NULL  UNIT,  I  LUt  of  null  device  for  logging  routine 

■f  STO_Unit  i  LUt  for  buffer  storage 

C 

c 

CHARACTER*40  DEF  I  default  directory 

LOGICAL  PROCFLAG 

LOGICAL  DEFINE_FLAG 

C 

COMMON/BLK2/COM  UNIT,  lOUT,  TERMUNIT,  FILE_UNIT, 

+  DEF,  LOgUNIT,  NULL_ONIT,  PROCFLAG, 

-f  DEFINE  FLAG 


& 


i 

g 


V 

V 


V, 

V 

V 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

CC  FTRN.FOB 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

COMMON  BUFFER, NDIM,ANST,AINC,FTYPE 
DIMENSION  ANPL( 5000) ,PHS(  5000) 

BYTE  BUFFER( 35000) 

INTEGER*2  INPFILE(15) 

LOGICAL  FTYPE 
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o  n  n 


ccccccccaccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

cc 

cc  DEFINE  THE  TAPE  HEADER  FIELDS 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc 

cc 

BYTE  ITITLE(50)  ITITLE 

INTEGER*2  IHED(6)  ! DATE  AND  TIME  AS  MONTH , DAY , YEAR 

C  1  HOURS, MINUTES. SECONDS 

BYTE  ITARG(38)  I  TARGET  LABEL 

INTEGER*2  IANG(3)  I  STARTING  ANGLE,  ANGLE  INCREMENT  FOR  ROTATION 

C  t  AND  NUMBER  OF  ANGLES  IN  WHOLE  FILE 

BYTE  ITYPE(6)  I  ACOUSTIC  DATA  TYPE  PCW.PLFM 

INTEGER*2  IPARAM(6)  I  FREQUENCY  IN  KHZ,  SAMPLE  INTERVAL,  INTERVAL 

C  I  UNITS,  AND  NUMBER  OF  PINGS  AT  A  GIVEN  ANGLE 

INTEGER*2  IDEANG  1  ELEVATION/DECLINATION  ANGLE 

INTEGER*2  IPULTH  1  PULSE  LENGTH  -  DIGITS  ONLY 

INTEGER*2  IPUNIT  lUNlTS  FOR  PULSE  LENGTH 

INTEGER*2  IMODFW  1  MODULATION  BANDWIDTH  FOR  PLFM 

INTEGER*2  ISTRNG  I  SOURCE-TARGET  RANGE  IN  METERS  *  100 

INTEGER*2  ISRRNG  1  SOURCE-RECEIVER  RANGE  IN  METERS  *  100 

INTEGER*2  IXMTVL  1  RMS  TRANSMIT  VOLTAGE  *  100 

INTEGER*2  IRCVGN  I  RECEIVER  GAIN  IN  DB  *  10 

INTEGER*2  IFTRBW  tRECEIVE  FILTER  3-DB  BW  IN  KHZ 

BYTE  IPROJ(20)  IPROJECTOR  DESCRIPTION 

INTEGER*2  ITRV  IXMIT  LEVEL  OF  PROJECTOR  IN  DB/MICRO  PA/V*10 

BYTE  IHYD(20)  IRECEIVER  DESCRIPTION 

INTEGER«2  IRRS  IRECEIVER  SENSITIVITY  IN  DB/V/MICRO  PA*10 

INTEGER*2  IDATR  INUMBER  OF  BIOMATION  SAMPLES  IN  A  SINGLE  PING 


INTEGER*2  HEADER(256) 

COMMON  /HEADER/  ITITLE, IHED, ITARG, lANG, ITYPE, IPARAM , IDEANG , I PULTH , 

2  IPUNIT , IMODFW , ISTRNG, ISRRNG, IXMTVL .IRCVGN , I FTRBW , 

3  IPROJ , ITRV, IHYD, IRRS, IDATR, I OTHER 
EQUIVALENCE  ( HEADER, ITITLE ) 


APPENDIX  C 


CIMAG2  LINKING  SUBROUTINES 


cc 

s 

cc 

cc 

cc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  INTER 

cc  FREQUENCY  DOMAIN  READ 

CC 

CCCCCCCC'CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc 

SUBROUTINE  FRDl ( P , EFLAG ) 

LOGICAL  EFLAG  !  error  flag 

INTEGER*2  HEAD(256) 

REAL  P(  4096,7) ,AM( 2049) , PH (2049) ,FA{  2049 ) 

INCLUDE  'USER2; (DURAL. CIMAG2)MAGCMN2. FOR'  I  program  control  common  block 

COMMON/HEADER/HEAD 

ISIZE-4096 

IS2-ISIZE/2. 

TYPE  FREQUENCY  SAMPLING(l)  OR  .IKL  SAMPLING  (  0 ) 7  ' 

READ  (COM  UNIT,*)  ITS 

WRITE  (LOS_UNIT,*)  IPS  !  log  the  input 

IF  (IFS.EQ.O)TYPE  *, 'INPUT  MAJOR  AXIS  DIMENSION  IN  INCHES' 

IF  (IFS.EQ.l)TYPE  *,' FREQUENCY  INCREMENT  IN  MHZ  ' 

READ  (COMUNIT,*)  DLE 

WRITE  (LOG  UNIT,*)  DLE  !  log  the  input 

PI-4 . *ATANT1 . ) 

IF  ( IFS .EQ. 1 )DLE-590. 551/(PI*DLE) 

TYPE  *,'  SELECT  THE  TYPE  OF  INTERPOLATION' 

TYPE  *, 'INPUT  0  — >  TWO-POINT  INTERPOLATION  ;  NO  SMOOTHING' 

TYPE  *,' INPUT  1  — >  INTERPOLATION  AND  SMOOTHING  USING  A  COSINE  WINDOW' 
READ  (COMUNIT,*)  IS 

WRITE  (LOG  UNIT,*)  IS  !  log  the  input 

TYPE  *,'  ASSUMED  INPUT  AMPLITUDE  IS  IN  DB/SQUARE  CM  ' 

TYPE  *,'  NORMALIZE  TO:  SQ  CM(1),SQ  M( 2 ) , PI *L*L/4 ( 0 ) 7 ' 

READ  (COM_UNIT,*)  INORM 

WRITE  (LOG_UNIT,*)  INORM  !  log  the  input 

IF  (INORM. GT. 2. OR. INORM. LT.O)  GO  TO  5 
DLE-DLE*2. 54/100. 

CALL  RDFLE(AM,PH,NP,FLOW,FHIGH,riNC,COM  UNIT, EFLAG) 

IF  (EFLAG)  THEN  "  I  If  an  error  occurred 

RETURN 

END  IF 

TWOPI-8.*ATAN(l. ) 

RLOW-TWOPI*rLOW*DLE/300. 

RHIGH-TWOPI*FHIGH*DLE/300. 

RINC-. 1 

WD-TWOPI*riNC*DLE/300 . *6 . 

DO  5  I-1,IS2 
P(  I-MS2,1  )-0 
5  P(I,1)— 1024. 

DO  10  I-1,NP 

10  rA( I )-(rLOW+( I-l )*FINC)*'rWOPI*DLE/300. 

CAIL  GP( INORM, IS,AH,PH,FA,NP,WD,DLE,FLOW,RHIGH,RINr ,p) 

RE'IURN 

END 


!  log  the  input 


!  log  the  input 


SUBROUTINE  GP( INORM , IS , AM, PH , FA,NP ,WD,DLE, RLOW, RHIGH , RINC , P ) 
REAL*4  FA( 2049 ) ,P( 4096, 7) ,AM( 2049) ,FH( 2049) 

COMPLEXES  R,CA(2049) 

PI-4.*ATAN(1. ) 

DO  111  I-1,NP 
ATMP-10.*«(AM( I )/20. ) 

PTMP-PH( I )/180.*Pl 

CA( I )-CMPLX(ATMP*COS(PTMP) ,ATHP*SIN(PTHP) ) 

TWOPI-PI*2 , 

I-l 

DO  7720  RF-RLOW, RHIGH, RINC 

IF(RF.GE.FA(I),AND.RF.LE.rA(I+l))GO  TO  7730 


85 


7730 


7720 

4912 

C 

C 


20 


10 

15 


30 


C 

C 


30 


20 


40 

10 


C 

C 


I-I  +  l 

ir( I .GT.NP)GO  TO  4912 

GO  TO  7740 

K-RF*10 

IF(K.GT.2048)GO  TO  4912 

IF( IS.EQ.l )CALL  INTER ( R , CA, FA, NP , NS , I .WD.RF) 

IF(IS.EQ.O)R-(CA(I+l)-CA(I))/(FA(I+l)-FA(I))*(Rr-FA(I) )+CA(I ) 
IF  ( INORM.EQ.0)R-R*0. 02/(SQRT(PI )*DLE) 

IF  (INORM.EQ.2)R-R/100. 

P(K,  1  )-10. 24*(  10.+20.*AI-OG10(CABS(R)  )  ) 

P(K+2048,1)-1024.*ATAN2(AIMAG(R) ,REAL(R) )/TWOPI 

CONTINUE 

RETURN 

END 


SUBROUTINE  I NTER ( R , CA , FA , ICNT , NS , I , WD , RF ) 
COMPLEX*8  R,CA{2049) 

REAL.*4  FA(2049) 

WEI>0 

IS-I 

XTMP-0. 

YTMP-O . 

RFL-RF-WD/2 . 

RFH-RF+WD/2. 

IF(IS.GT.ICNT)GO  TO  10 
ir(FA(IS) .GT.RFH)GO  TO  10 
T-FA( 1S)-RF 

HAMM-.54+.46*COS( 3.1415926*T/WD) 
XTMP-XTMP+HAMM*REAL(CA( IS) ) 
YTMP-YTMP+HAHM*AIMAG(CA( IS) ) 

WEI«WEI>HAMH 
IS-IS+l 
GO  TO  20 
IS-I-1 

ir{ IS.LT.l )GO  TO  30 
ir(FA(IS)  .LT.RFDGO  TO  30 
T-FA( IS)-RF 

HAMM-.54+.46*COS(3.1415926*T/WD) 
XTMP-XTMP+HAMM*REAL(CA( IS) ) 
YTMP-YT«P+HAHM*AIMAG(CA(IS) ) 

WEI-WEI+HAHM 
IS-IS-1 
GO  TO  15 

R-CHPLX( XTMP/WEI , YTMP/WEI ) 

RETURN 

END 


StiBBOUTTflF  PHrRfpPH,  FB,riNC,PH) 
PFAI*'!  Pll(201) 

PI-3.1415926 
IF(DPH.LT.PI )GO  TO  20 
DPH-DPH-Pl 
GO  TO  30 

IF(DPH.GT,-PI )GO  TO  40 

DPH-DPH+PI 

GO  TO  20 

DO  10  1-1,201 

FA-FB+(I-l)*riNC 

PH( I )-PH( I )-(DPH*rA/FB) 

RETURN 

END 


SUBROUTINE  EPH( PHH , PH, NS ,NF, NI ) 


86 


( 


REAL*4  PH(201) 

WEI-0 

PHM-0 

WD-2.*(Nr-NS) 

DO  10  I-NS,NF,NI 

HAMH-. 54+.46*COS( 3. 1415926* (I-NS)/WD) 
WE  I -WE  I -t- HAHN 
PHM-PHM+HAMM*PH(  I ) 

PHH-PHM/WEI 

RETURN 

END 


(Jjl 

ijjl 
iJJ 

J 

m 

•I' 

«*1 

•I' 

•*l 


UH 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  RDFLE 

cc  FREQUENCY  DOMAIN  READ 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc 

cc 

SUBROUTINE  RDFLE ( A , P , NP , FMIN , FMAX , FINC , COM_UNIT , EFLAG ) 

CHARACTER*!  YN 
LOGICAL  EFLAG 
COMMON  BUFF 
BYTE  BUFF(20000) 

INTEGER* 2  LINEK  30) ,LINE2( 30) ,PARAM( 30) ,OFILE(15) 

INTEGER*4  COM_UNIT 

CHARACTER*70  INFILE 

REAL*4  AP(4098) ,A(2049) ,P(2049) 

INTEGER*2  HEADR(2S6) 

COMMON/HEADER/HEADR 

C 

C  DEFINE  BUFFER  STRUCTURE 

C 

EQUIVALENCE (LINE! ( 1 ) ,BUFF(1) ) , (LINE2( 1 ) .BUFF  (61 ) ) 

1,  (PARAMd  )  ,BUFr(  121)  )  ,  (AP(1  )  ,BUFF(  361  )  ) 

CALL  TR( INFILE, IB, EFLAG) 

IF  (  EFLAG  )  THEN  1  if  error  has  occurred 

RETURN  t  then  return 

END  IF 

TYPE  lOB.LINEl 
TYPE  105,LINE2 
TYPE  lOS.PARAM 
105  FORMAT(X,30A2) 


PUT  INFO  IN  HEADR  BLOCK  FOR  TRANSFER  TO  CALLING  PROG 


1 _ 


DO  108  1-1,30 
HEADR(I)«LINE1(I) 
HEADR(  30-t-I  )-LINE2(I  ) 
HEAOR(60>I )-PARAM( I ) 


GET  NUMERICAL  INFORMATION  FROM  THE  THIRD  LINE 
OF  THE  HEADER 

CALL  DCDE(NP, FMIN, FINC, EFLAG) 

IF  (EFLAG)  THEN 
RETURN 

END  IF 

DIVIDE  AN  AMP-PHASE  ARRAY  INTO 
AN  AMP  ARRAY  AND  A  PHASE  ARRAY 

DO  199  NN-1,NP 
A(NN)-AP( 2*NN-1 ) 

P(NN)-AP( 2*NN) 

FORMAT(Al ) 

rORMAT(lX,5(2rl2.3,lH;)  ) 

FMAX-FMIN+(NP-1 )*FINC 

RETURN 

END 


i f  error  return 


‘ 


'f, 


'•5 

t 


SUBROUTINE  TR( INFILE , IB , EELAG) 

INCLUDE  'MAGCMN2 . FOR’ 

LOGICAL  EFLAG  !  error  flag 

INTEGER*2  IBUFF( 10000) 

CHARACTER*70  INFILE 

INTEGER*4  FOR  RETCODE  1  fortran  return  code 

INCLUDE  ' USERli ( DURAL . CIMAG2 JMSGBLK  .  FOR'  I  error  message  declarations 

COMMON  BUFF 

BYTE  BUFF(20000) ,TBUFF(1500) 

EQUIVALENCE  { BUFF ( 1 ) , I BUFF ( 1 ) ) 

INCLUDE  'SYS$LIBRARY:FORIOSDEF' 

WRITE(6,S) 

FORMAT(lX, 'TYPE  DATA  FILE  NAME') 

READ  (COM_UNIT,10)  INFILE 
WRITE  (LOG_UNIT,10)  INFILE 

IF  ((INFILE  .EQ.  'NAME.DAT')  .AND.  DEFINE_FLAG)  THEN  I  if  defining 

a  procedure  then 

INFILE-'USER2;(DURAL.CIMAG2)NAME.DAT'  !  use  the  dummy  file 

END  IF 
FORMAT (  A  ) 

IB-1 

ICNT-0 

OPEN(UNIT-8 , NAME- INFILE, TYPE- 'OLD' , READONLY , IOSTAT-FOR_RETCODE , 
DEFAULTFILE-DEF.ERR-Sl) 

GOTO  82 

IF  (FOR_RETCODE  .EQ.  29)  THEN  I  if  file  is  not  found 

CALL  LIB$SIGNAL(MAG_PILNOTFOU)  1  write  out  file  not  found  error 

EFLAG  -  .TRUE.  1  error  flag 

RETURN 


CALL  L I B$ SIGNAL (HAG  COM) 
EFLAG  -  .TRUE.  ~ 

RETURN 


1  else  command  error 
t  flag  the  error 


END  IF 


S- 

If- 

S' 

E 


^  »*• 


SET  BLOCK  LENGTH  IN  BYTES 

IF( IB.EQ.l )LEN-512-9*4 
1F( IB.GT.l )LEN-512-26*4 

READ  A  BLOCK  OF  512  BYTES 

READ( 8 ,80,END-90) ( TBUFF( I ) , I-l , 512 ) 

FORMAT( 512A1 ) 

STORE  A  BLOCK  INTO  THE  BUFFER  ACCORDING  TO  ITS  LENGTH 

DO  85  I-1,LEN 

BUFF( ICNT+I )-TBUFF( 1 ) 

Tn-TB*1 

ICNT-ICMT+LEN 

GO  TO  82 

DO  86  I-1,LEN 

BUFF( ICNT+I )-TBUFF( I ) 

ELIMINATE  BLANK  SPACES  IN  BETWEEN  EACH  CHARACTER 
IN  A  FILE  HEADER 

DO  40  1-1,180 
BUFF( I )-BUFF( 2*1-1 ) 

CL0SE(UN1T-8,DISP-’SAVE' ) 

RETURN 

END 


SUBROUTINE  DCDE ( NP , FHIN, FINC, EFLAG ) 


89 


COMMON  BUFF 
BYTE  BUFF(20000) 
INTEGER*4  IMIN,IINC,NP 


s 


NO  OF  DATA  POINTS  IS  STORED  IN  THREE  CHARACTERS,  AND 
STARTING  FREQ  AND  FREQ  INC.  IN  5  CHARACTERS 


100 


110 


CHARACTER* 4  CNL 
CHARACTER*5  CFF.CINC 
INTEGER* 4  FOR_RETCODE 
LOGICAL  EFLAG 

INCLUDE  ’USER2: (DURAL. CIHAG21MSGBLK. FOR’ 


1  fortran  return  code 
1  error  messages 
EQUIVALENCE  (BUFF (123) , CNL ) , ( BUFF ( 1 31 ) ,CFF) , ( BUFF ( 140) ,CINC) 


CONVERT  CHARACTERS  INTO  THEIR  NUMERICAL  EQUIVALENTS 

if  an  equal  sign  apperars  as  th  first  character 
change  it  in  to  a  blan)(. 

IF(  CNL(  1:1  )  .EQ.  )  CNL(  1:  1  )  -  '  ' 

READ(UNIT-CNL,  FMT-100,  IOSTAT-FOR_RETCODE ,  ERR-110)  NP 
FORNAT( 15) 

READ(UNIT-CFF,  FMT-100,  IOSTAT-FOR_RETCODE ,  ERR-110)  IMIN 

READ(UNIT-CINC,  FMT-100,  10STAT-F0R_RETC0DE ,  ERR-110)  IINC 

FMIN-IMIN 

FINC-IINC 

RETURN 

IF  (FOR_RETCODE  . EQ .  64)  THEN  I  if  input  format  error 

CALL  LIB$SIGNAL(MAG  INPFOR)  1  PRINT  ERROR 


ELSE 

EFLAG  -  .TRUE. 

1  flag 

an  error 

CALL  LIB$SIGNAL(HAG  COM) 

1  else 

some  other 

END  IF 
RETURN 
END 

EFLAG  -  .TRUE. 

I  flag 

error 

90 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  FORT 

cc  FOURIER  TRANSFORM 

cc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 


SUBROUTINE  FORT { A , H , S , I FS , I FERR ) 

DIMENSION  A( 1 ) ,S( 1 ) ,K( 15) 

IF(M)2,2,3 
3  IF(M-14)  5,5,2 
2  IFERR-1 
1  RETURN 
5  IFERR-0 
N-2**M 

IF(  lABS(lFS)  -  1  )  200,200,10 

WE  ARE  DOING  TRANSFORM  ONLY.  SEE  IF  PRE-COMPUTED 
S  TABLE  IS  SUFFICIENTLY  LARGE 
10  IF(  N-NP  )20,20,12 
12  IFERR-1 
GO  TO  200 

SCRAMBLE  A,  BY  SANDE'S  METHOD 
20  K(1)-2*N 
DO  22  L-2,M 
22  K(L)-K(L-l)/2 
DO  24  L-M,13 
24  K(L+l)-2 

BINARY  SORT- 
Kl-K( 14  ) 

K2-K( 13) 

K3-K(12) 

K4-K( 11 ) 

K5-K(10) 

K6-K(9) 

K7-K(8) 

K8-K( 7 ) 

K9-K( 6 ) 

K10-K( 5) 

Kll-K( 4 ) 

K12-K( 3) 

K13-K( 2 ) 

K14-K( 1 ) 

N2-K( 1 ) 

IJ-2 

DO  30  J1-2,K1,2 
DO  30  J2-J1,K2,K1 
DO  30  J3-J2,K3,K2 
DO  30  J4-J3,K4,K3 
pn  30  J5-J4,k5,K4 
I’U  3'J  J6-J5,K6,K5 
DO  30  J7-J6,K7,K6 
DO  30  J8-J7,K8,R7 
DO  30  J9-J8,K9,K8 
DO  30  J10-J9,K10,K9 
DO  30  J11-J10,K11,K10 
DO  30  J12-J11,K12,K11 
DO  30  J13-J12,K13,K12 
DO  30  JI-J13,K14,K13 
IF( IJ-JI )28, 30, 30 
28  T-A(IJ-1  ) 

A( IJ-1 )-A( JI-1 ) 

A( JI-1 )-T 
T-A( IJ ) 

A( IJ )-A( JI  ) 

A( JI  )-T 


FORT  001 

FORT  056 

FORT  058 
FORT  059 
FORT  060 
FORT  061 
FORT  062 
FORT  063 
FORT  064 
FORT  065 
FORT  066 
FORT  06' 
FORT  068 
FORT  069 
FORT  030 
FORT  071 


FORT  075 


FORT  076 
FORT  077 
FORT  078 
FORT  079 
FORT  080 
FORT  OR' 
FORT  0  6.; 
FORT  083 
FORT  084 
FORT  085 
FO''T  086 
FORT  087 
FORT  088 


FORT  090 
FORT  091 
FORT  092 
FORT  093 
FORT  094 
FORT  (709 
FORT  096 


wv 
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CIHAG2  THE  COMPUTER  PROGRAM  TO  GENERATE  COLOR  IMAGES 
(U>  OHIO  STATE  UNIV  COLUMBUS  ELECTROSCIENCE  LAB 
G  DURAL  ET  AL  NOV  87  ESL-718e48-7  Ne8ai4-BG-K-a2e2 

F/G  12/5  NL 


30  lJ-IJ+2 

IF( 1FS)32,2,36 

C  DOING  FOURIER  ANALYSIS, SO  DIV.  BY 

32  FN  -  N 

DO  34  I-1,N 

A(2*I-1)  -  A(2*I-l)/rN 
34  A(  2*1  )— A(  2*1  )/FN 
C  SPECIAL  CASE-  L-1 

36  DO  40  I-1,N,2 
T  -  A(2*l-1) 

A(2*I-1)  -T  +  A(2*I+1) 
A(2*I+1)-T-A(2*I+1) 

T-A( 2*1 ) 

A(2*I)  -  T  +  A(2*I+2) 

40  A(2*I+2)-  T  -  A(2*I+2) 

IF(M-l)  2,1  ,50 

C  SET  FOR  L-2 

50  LEXPl-2 

C  LEXP1-2**(L-1) 

LEXP-8 

C  LEXP-2**( L+1 ) 

NPL-  2**MT 

C  NPL  -  NP*  2**-L 

60  DO  130  L-2,H 
C  SPECIAL  CASE-  J-0 

DO  80  I-2,N2,LEXP 

11- I  +  LEXPl 

12- 11*  LEXPl 
13  -I2*LEXP1 
T-A( I-l ) 

A(I-l)  -  T  ♦A(12-l) 

A(I2-1)  -  T-A(I2-1) 

T  -A( I  ) 

A(I)  -  T*A(12) 

A(I2)  -  T-A(I2) 

T-  -A(I3) 

TI  -  A(I3-1) 

A(I3-1)  -  A(Il-l)  -  T 
A(I3  )  -  A(I1  )  -  TI 

A(ll-l)  -  A(ll-l)  +T 
80  A(I1)  -  A(I1  )  +T1 

IF(L-2)  120,120,90 
90  KLAST-N2-LEXP 
JJ-NPL 

DO  110  J-4,LEXF1,2 

NPJJ-NT-JJ 

UR-S(NPJJ) 

UI-S( JJ) 

1LAST-J4KLAST 

DO  rno  r-  j,ilast,lexp 

11- l+LEXPl 

12- I1*LEXP1 

13- I2*LEXF1 

T-A( 12-1 )*UR-A( I2)*UI 
TI-A( I2-1)*UI+A( I2)*UR 
A(I2-1)-A(I-1)-T 
A(I2  )-A(I  )  -  TI 

A(I-l)  -A(I-1)+T 
A(I)  -A(I)+TI 
T— A(  13-1  )*UI-A(  I3)*UR 
TI-A(I3-1)*UR-A( 13)*UI 
A(13-1)-A(I1-1)-T 
A(I3)  -A(ll  )-TI 
A(I1-1)-A(I1-1)+T 
100  A(I1)  -A(ll)  +T1 

C  END  or  I  LOOP 


FORT  097 
FORT  098 

M  AND  CONJUGATE.  FORT  099 

FORT  100 
FORT  101 
FORT  102 
FORT  103 
FORT  104 
FORT  105 
FORT  106 
FORT  107 
FORT  108 
FORT  109 
FORT  110 
FORT  111 
FORT  112 
FORT  113 
FORT  114 
FORT  115 
FORT  116 
FORT  117 
FORT  118 
FORT  119 
FORT  120 
FORT  121 
FORT  122 
FORT  123 
FORT  124 
FORT  125 
FORT  125 
FORT  12~ 
FORT  128 
FORT  129 
FORT  130 
FORT  131 
FORT  132 
FORT  133 
FORT  134 
FORT  135 
FORT  136 
FORT  137 
FORT  138 
FORT  139 
FORT  140 
FORT  141 
FORT  142 
FORT  143 
FORT  ’44 
FORT  145 
FORT  ) 

F’  . 

FORT  148 
FOPT  140 
FORT  150 
FORT  151 
FORT  152 
FORT  153 
FORT  154 
FORT  155 
FORT  156 
FORT  157 
FORT  158 
FORT  159 
FORT  160 
FORT  161 
FORT  162 


n  o 


110  JJ-JJ+NPL 
C  END  OF  J  LOOP 

120  LEXP1-2*LEXP1 
LEXP  -  2*LEXP 
130  NPL-NPL/2 
C  END  or  L  LOOP 

140  IF( irS)145,2,l 

C  DOING  FOURIER  ANALYSIS.  REPLACE  A  BY  CONJUGATE. 

145  DO  150  1-1, N 
150  A(2*I)  — A(2*I) 

160  GO  TO  1 
RETURN 

MAKE  TABLE  OF  S ( J )-SlN( 2*PI* J/NP ) , J-1 , 2 - MT-l,NT-NP/4 

200  NP-N 
MP-M 
NT-N/4 
MT-M-2 

IF(MT)  260,260,205 
205  THETA-. 7853981634 
c  theta-pi/2*»(l+1)  for  L-1 
210  JSTEP  -  NT 

C  JSTEP  -  2**(  MT-L+1  )  FOR  L-1 

JDIF  -  NT/2 

C  JDIF  -  2**(MT-L)  FOR  L-1 

S(JDIF)  -  SIN(THETA) 

IF  (MT-2)260,220,220 
220  DO  250  L-2,MT 

THETA  -  THETA/2. 

JSTEP2  -  JSTEP 
JSTEP  -  JDIF 
JDIF  -  JDIF/2 
S( JDir)-SlN(THETA) 

JCl-NT-JDir 
S(JCl) -COS (THETA) 

JLAST-NT-JSTEP2 
IF (JLAST-JSTEP >250,230, 230 
230  DO  240  J-JSTEP,JLAST, JSTEP 
JC-NT-J 
JD-J+JDir 

240  S( JD)-S( J)»S( JC1)+S( JDIF)*S(JC) 

250  CONTINUE 
260  ir(lFS)20,l,20 
END 
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FORT  171 
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FORT  173 
FORT  174 
FORT  175 
FORT  176 
FORT  177 
FORT  178 
FORT  179 
FORT  180 
FORT  181 
FORT  182 
FORT  183 
FORT  184 
FORT  185 
FORT  186 
FORT  187 
FORT  188 
FORT  189 
FORT  190 
FORT  191 
FORT  192 
FORT  193 
FORT  194 
FORT  195 
FORT  196 
FORT  197 
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FORT  199 
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FORT  201 
FORT  202 
FORT  203 
FORT  204 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  DDHPB 

cc  BASIC  SCATTERING  CODE  READ  ROUTINE 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc 

SUBROUTINE  BSCREA(P,COH  UNIT) 

COMPLEX  ET( 1024 ) ,EP( 1021) 

REAL  P( 4096,9) ,DT(1024 ) .DP(1024) 

CHARACTER*20  INFILE 

INTEGER *4  COM_UNIT 

INCLUDE  'SyS$LIBRARY:FORIOSDEP' 

PI-3.14159265 

810  TYPE  INPUT  FILE  NAME  -  ' 

READ  (COM  UNIT, 7823)  INFILE 
7823  FORMAT(A2iy) 

OPEN (UN IT-19,  NAME- 1 NF I LE , STATUS- ' OLD ' , FORM- ' UNFORMATTED ' 

2, IOSTAT-1ERR,ERR-8100) 

DO  3410  1-1,2048 
P( 1+2048, l)-0. 

3410  P(I,1  )— 1024. 

READ(19)  ET 
READ(19)  EP 

REAO(19)  1B,IE,IS,PLS,PLI 
FAC-10 . +BLOG10 ( 4 . *PI ) 

DO  40  I-IB,IE,IS 

DT( I )-20. +BLOG10(BABS(ET( I ) ) )+FAC 
40  DP( I )-20. *BLOG10(BABS(EP(I ) ) )+FAC 

ETM--200. 

EPM— 200. 

ETA-0. 

EPA-0. 

IT-0 

DO  11  I-IB, IE, IS 

1F(DT( I) .GT.ETM)  ETM-DT(l) 

IF(OP(I) .GT.EPH)  EPM-DP(I) 

IT-IT+1 
ETA-ETA+DT( I ) 

EFA-EPA+DP( I ) 

11  CONTINUE 

ETA- ETA/IT 
EPA-EPA/IT 
PLE-PLS+(IE-1)*PLI 
TYPE  1,PLS,PLE,PLI 

1  FORMAT ( '  START-  ' , FI 0 . 5 , 5X , ' END-  ' , FI 0 . 5 , 5X , ' STEP-  ’,F10.5) 

TYPE  3,ETM,EPM 

3  FORMAT ( '  ETM-  ' , FIO . 5 , 5X, ' EPM-  ',F10.5) 

TYPE  2,ETA,EPA 

2  FORMAT('  ETA-  ' , FI 0 . 5 , 5X , ' EPA-  ',rl0.5) 

TYPE  4 

4  FORMATC  ET-1,  EP-2,  IPOL-  ',$) 

READ  (COM  UNIT,*)  .IPOL 
DLE-0.03/T2.*PI*PLI ) 

TYPE  *, '  L-  ' ,DLE 

TYPE  *, 'NORMALIZE  TO:  SQ  CM(1),  SQ  H(2),  PI *L*L/4 ( 0 ) 7 ' 

READ  (COMUNIT,*)  INORM 
DFAC-0. 

IF( INORM. EQ.O)  DFAC— 10.*BLOG10(0.25*PI*DLE*DLE) 

IP( INORM. EQ.l )  DFAC-20. 

KLOW-1 . 5+PLS/PLI 

DO  7720  I-IB, IE, IS 

K-KLOW+I-1 

IF( IPOL.EQ. 1 )  THEN 

P( K, 1 )-10.24*( 10.+DFAC+OT( I ) ) 

P(R+2048,1)-512.*BTAN2(AIHAG(ET(I) ),REAL(ET(I) ) )/PI 


ELSE 

P(K,1)-10.24*(10.+DFAC+DP(I)) 

P(K-f  2048,1) -512.  *BTAN2(AIMAG(EP(  I)  )  ,REAL(EP(I)  )  )/PI 
ENOIF 

7720  CONTINUE 
GO  TO  80 

8100  IF( IERR.EQ.POR$IOS_FILNOTPOU)THEN 

TYPE  1112, INFILE 

1112  FORMAT('  FILE  :  •,A20,'  DOES  NOT  EXIST',// 

2,'  ENTER  FILENAME  AGAIN') 

ELSE  IF  ( IERR.EQ.FOR$IOS_FlLNANSPE)THEN 

TYPE  *, 'FILE:' , INFILE, 'HAS  BAD,  ENTER  NEW  FILENAME' 

ELSE 

TYPE  *, 'UNRECOVERABLE  ERROR  ,  CODE  -',IERR 

STOP 

ENDIF 

GO  TO  810 

80  CLOSE(UNIT-19,DISP-'SAVE' ) 

RETURN 

END 

I NCLUDE  ' ESL_ESLROOT : { GRPl ILIBIBFILES. FOR ' 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cc 

CC  CINAGHSG 

cc  PROGRAM  ERROR  MESSAGES 

CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc 

cc 

•TITLE  CIMAG_MESSAGES  Program  error  meesagea  1  module  name  and 

C  listing  title 

•FACILITY  CIMAG,2a3/PREFIX-MAG_ 

1 

1 

I  Error  Messages 

•SEVERITY  ERROR 

COM  "error  in  command* 

INPFOR  "error  in  input  format" 

FILNOTFOU  "file  not  found  in  this  directory* 

•  END 


^  n  o  o  o  o  n  oooooooo 


APPENDIX  D 


COLOR  IMAGING  PROGRAM  CLRPL 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
cc  cc 

CC  CLRPL  cc 

cc  cc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


THIS  PROGRAM  DEFINES  A  PLOTTING  SOFTWARE  CODE  TO  PLOT 
IMAGE  OF  AN  OBJECT  WITH  POLARIZATION  INFORMATION 

IMAGE  DATA  SHOULD  BE  CODED  BY  A  100*100  MATRIX  AND 
IT  SHOULD  BE  CONTAINED  IN. THE  MATRIX  'ARRAYC 100 , 100 ) ' 

DIMENSION  ICOL(100,100) ,ICOLA(100) 

DIMENSION  GX( 5) ,GY(5) ,CBX(5) ,CBY(5) ,XP( 5) ,YP( 5) ,XT(10) ,YT(10) 
DIMENSION  ARRAYdOO.lOO)  ,COORD(S) 

INTEGER  COLI.WRSTID 
REAL* 4  L 
LOGICAL  FINISHED 
CHARACTER*40  FILNM 
CHARACTER*?  NUMB 

THE  GX  MATRIX  IS  THE  X  COORDINATES  OF  THE  COLOR  GRID 
DATA  GX/. 30,. 85,. 85,. 30, .30/ 

THE  GY  MATRIX  IS  THE  Y  COORDINATES  OP  THE  COLOR  GRID 
DATA  GY/. 12, .12, .67, .67, .12/ 

THE  CBX  MATRIX  IS  THE  X  COORDINATES  OF  THE  COLOR  BOX 
DATA  CBX/. 30, .85, .85, .30, .30/ 

THE  CBY  MATRIX  IS  THE  Y  COORDINATES  OP  THE  COLOR  BOX 
DATA  CBY/. 70, .70, .72, .72, .70/ 

THE  XT  MATRIX  IS  THE  X  COORDINATES  OF  THE  TIK  HARKS 
DATA  XT/. 3, .44, .57, .71, .85, .26, . 26 , . 26 , . 26 , . 26/ 

THE  YT  MATRIX  IS  THE  Y  COORDINATES  OF  THE  TIK  HARKS  ON  THE  Y  AXIS 
DATA  YT/. 08, .08, .08, .08, .08, . 12 , . 26 , . 40, . 54 , . 67/ 


P1-4.*ATAN(1. ) 

KHK-412505 

WKSTID-1 


VMTrr — V  r.f  MnQnjt-iidr  Top  r^nt-^r  line 

YIITCL-CBVl  3  ;+ .02 

ENTER  THE  ARRAY  TO  BE  PLOTTED 
WRITE(6,*) 'ENTER  FILE  NAME' 

READ( 5,100)  FILNM 
)0  FORHAT(A40) 


POLARIZATION  DEPENDENT  PLOT? 

HRITE(6,*)  'DO  YOU  NEED  TO  USE  POLARIZATION  AS  A  PARAMETER?' 
WRITE(6,*)'(Y-1  N-O)' 

REAO(5,*)  POL 
IF(POL.EQ.l)  THEN 
C  ENTER  THE  POLARIZATION 

WRITE(6,*) 'ENTER  POLARIZATION  <VP-l,HP-2)' 

READ(5,*)  PL 


OPEN ( UN I T«1 , r I LE-r I LNH , STATUS- ' OLD ' , FORH- ' UNFORMATTED ' ) 

ENTER  THE  IMAGE  SIZE  AND  PERIOD  OF  THE  TIME  SIGNAL 
READ(l)  IMGSZ.PER 

ENTER  THE  IMAGE  ARRAY 

DO  1000  1-1,100 
DO  1000  J-1,100 
READ( 1 )ARRAY( I ,J) 

ARRAY(1,J)-ABS(ARRAY(I,J)) 

CONTINUE 


NORMALIZE  THE  ARRAY  VALUES 

CALL  SEARCH(ARRAY,AMAX,AMIN) 

WRITE(6,* ) 'MAX-' ,AMAX, 'MIN-' ,AMIN 
WRITE(6,*) 'ENTER  THE  DESIRED  MAX. .AND, MIN. ' 

READ(5,*)  ANHAX.ANMIN 
DO  20  1-1,100 
DO  20  J-1,100 

IF  (ARRAY( I , J) .GT.ANMAX)  ARRAY ( I , J ) -ANMAX 
IF  ( ARRAY ( I, J).LT. ANNIN)  ARRAY( I , J )-ANMIN 
ARRAY ( I , J ) - ( ARRAY ( I , J ) -ANNIN )/( ANMAX- ANM IN ) 

CONTINUE 

OPEN  GKS  ERROR  FILE 
CALL  GOPKS(6,5000) 

FIND  CONNECT  ID 

CALL  GKHGCU'ESL  4129' ,JER»OR,RCONID> 

IF( JERROR.NE.O)  THEN 

WRITE{6,*)  'Can  not  be  a  connection  ID' 

WRITE(6,*)  'Would  you  lilce  to  wait  ?  (Y-1)' 
READ(5,*)ANS 
IF(ANS.NE.l)  THEN 
STOP 
ELSE 

WR1TE(6,*)  'Enter  1  when  ready' 

REAO(5,*)  ANS 
GO  TC  999 
END  IF 
END  IF 

OPEN  WORKSTATON  11 
CALL  GOPWR( 1 ,RCONID,KHK) 

A'JIIVATE  WORKSTATION  II 
CALL  GACWR(l) 

SET  THE  WORKSTATION  WINDOW/VIEWPORT-FULLSCREEN 

GET  MAX  X  AND  Y 

KUNITS-O 

CALL  GODSP ( KWK , KERROR , KUNITS , XSI ZE , YSI ZE , RRASX , KRASY ) 
CALL  GSWKWN(1,0. ,1. ,0. . YSIZE/X5IZE ) 

CALL  GSWKVPd.O.,. 343,0. ,.274) 

GENERATE  THE  COLOR  INDICES 


IF(POL.EO.O)  THEN 
DO  3  COLI-l.lOO 


IF(COLI .LE.5)L-0 

H-95+2.6*COLI 

CALL  HLSRGB(H,L,S,R,G,B) 

CALL  GSCR(1,1+C0L1,R,G,B) 

ICOLA(COLI )-l+COLl 

CONTINUE 

ELSE 

DO  4  COLI-1,100 
L-.5 
S**l 

rC-25*(4.8)**(COLI/100. ) 

IF(COLI .LE.5)  L-0 
IF(PL.EQ.l)  H-225-FC 
IF(PL.EQ.2)H-110-FC 

I F ( PL . EQ . 2 . AND . COLI . LT . 2  5 . AND . COLI . GT . 5 ) L- . 7  5 

IF(H.LT.O)  H-H+360 

CALL  HLSRGB(H,L,S,R,G,B) 

CALL  GSCR(l,l-)-COLI,R,G.B) 

ICOLA( COLI )-l+COLl 
CONTINUE 
END  IF 

GENERATE  COLOR  CODE  FOR  X  AND  Y  COORDINATED  RCS  LEVELS 
DO  30  1-1,100 
DO  30  J-1,100 

ICOL ( I . J ) -ARRAY ( I , J ) *99+2 
CONTINUE 

PLOT  COLOR  LABEL  USING  CELL  ARRAY 
CALL  GCRSG(SEG) 

CALL  GCA(CBX(1),CBY(1),CBX(3),CBY(3) ,100,1,1,1,100,1,ICOLA) 
CALL  GSPLCI(l) 

CALL  GPL(5,CBX,CBY) 


LABEL  COLOR  BAR 
CALL  GSTXAL(2,3) 

CALL  GSTXP(O) 

CALL  GSCHXP(1.25) 

CALL  GSCHSPd.  ) 

CALL  GTX(GX(1),YMTCL,'0.0' ) 

CALL  GTX( .44,YMTCL, '0.25' ) 

CALL  GTX( .57,YMTCL, '0.50' ) 

CALL  GTX( .71,YMTCL, '0.75' ) 

CALL  GTX(GX(2),YMTCL, '1.0' ) 

CALL  GSCHXPd.) 

CALL  GSCHSPd.) 

PLOT  COLOR  NATIX  USING  CELL  ARRAY 
TFST-1 0 

IF  (TEST.FQ.O)GO  TO  21 

CALL  GCA(GXd),GYd),GX(3),GY(3),100,100,l,l,100,100,ICOL) 
CALL  GSPLCId) 

CALL  GPL(5,GX,GY) 

DRAW  GRID  LINE  ADJACENT  TO  EACH  AXES 

XPd)-.30 
XP(2)-.85 
YPd)-.08 
YP(2)-.08 
CALL  GPL(2,XP,YP) 

XPd)-.26 

XP(2)-.26 


CALL  GPL(2,XP,YP) 


C  PRINT  THE  TIK  HARKS 

CALL  GSMK(2) 

CALL  GSHRSC(l) 

CALL  GSPMCI(l) 

CALL  GPM( lO.XT.YT) 

C 

C  SHOW  THE  ASPECT  ANGLES 

XP(l)-.57 
XP(2)-.57 
YP(1)-.16 
YP(2)-.19 
CALL  GPL(2,XP,YP) 

XP(l)-.565 

XP(2)-.57 

XP(3)-.575 

YP(1 )-.185 

YP(2)-.19 

YP( 3)-.185 

CALL  GPL( 3,XP,YP) 

XP( 1 )-.80 
XP(2)-.77 
YP(l)-.40 
YP( 2)-.40 
CALL  GPL(2,XP,YP) 

XP(l)-.775 

XP(2)-.77 

XP( 3)-.775 

YP(l)-.405 

YP(2)-.4 

YP(3)-.395 

CALL  GPL(3,XP,YP) 

CALL  GTX( .57, .14, '0  DEG.') 

CALL  GTX( .57, .02, 'TIME  IN  NANOSECONDS’) 

CALL  GSTXAL(2,1) 

CALL  GSTXPIO) 

CALL  GSCHUP(1.,0.) 

CALL  GTX( .82, .4, '90  DEG.') 

CALL  GTX( .20, .40, 'TIME  IN  NANOSECONDS') 

C 

C  FIGURE  OUT  THE  COORDINATES  ON  EACH  AXES 

CALL  NARRSdNGSZ, PER, COORD) 

C 

C  PRINT  THE  COORDINATES 

X-.30 

CALL  GSCHUP(0. ,1. ) 

DO  40  1-1,5 

WRITE(NUHB,rMT-' (Ff .2) ' )COORD( I ) 

CALL  GTX(X, .06,NUHB) 

40  X-X-)-.14 

CALL  GSCHUPd.  ,0.  ) 

Y-.12 

DO  50  1-1,5 

WRI TE ( NUMB , FHT- ’ ( F5 . 2 ) ' ) COORD ( I ) 

CALL  GTX( .23,Y,NUHB) 

50  Y-Y+.14 

C 

C 

C  COLOR  1  -  WHITE,  0  -  BLACK  (ON  PLOTTER,  REVERSED  ON  SCREEN) 

C 

21  CALL  GCLSG(SEG) 

C 

WRITE( 6, * ) ' Enter  return  to  finish . ' 

READ( 5,1) FINISHED 


2 


FORMAT (Al ) 

CALL  GCLRWK(WKSTID, 1 ) 

DEACTIVATE  THE  WORK  STATION 
CALL  GDAWK(l) 


IClears  the  screen  on  TEKTRONIX 


CLOSE  THE  WORK  STATION 
CALL  GCLWK(l) 

CLOSE  THE  SYSTEM 
CALL  GCLKS 
STOP 
END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CC  THIS  ROUTINE  MAKES  THE  TRANSFORMATION  BETWEEN  COLOR  SYSTEMS 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  RLSRGB ( H , L , S , R,G,B ) 

REAL  H,L,S,R,G,B,M1,M2 
IF  (L  ,LE.  .5)  THEN 
M2-L*(l+S) 

ELSE 

M2-L+S-L*S 
8  END  IF 

M1-2*L-M2 

B«r9b_value  ( Ml  ,H2  ,H-)-120  ) 

R-rgb_value ( Ml ,M2 ,H ) 

G-rgb_value(Ml ,M2,H-120) 

RETURN 

END 


FUNCTION  rgb  value(Nl ,N2,HUB> 

REAL  rgb_value,Nl,N2,HUE 
IF  (HUE  .GT.  360)  THEN 

hue-hue-360 

END  IF 

IF  (HUE  .LT.  0)  THEN 
HUE-HUE+360 
END  IF 

IF  (HUE  .LT.  60)  THEN 

rgb_value-Nl+ ( N2-N1 ) *HUE/60 
ELSE  IF  (HUE  .LT.  180)  THEN 
rgb_value-N2 

ELSE  IF  (HUE  .LT.  240)  THEN 

rgb_value-Nl+(N2-Nl ) • ( 240-HUE )/60 
ELSE 

rgb  value-Nl 

end' I f 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  THIS  ROUTINE  FINDS  THE  MAX.,  AND  MIN.  OF  A  100*100  ARRAY 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 


SUBROUTINE  SEARCH( ARRAY, AHAX, AMIN) 
DIMENSION  ARRAY(100,100) 

AMAX— 1000 
AMIN-1000 
DO  1  1-1,100 

DO  1  J-1,100 


o  n 


IF(ARRAy(I,J) .GT.AMAX)AHAX-ARRAY(I,J) 

IF ( ARRAY ( I , J) .LT.AniN)AMIN-ARRAY( I , J) 

1  CONTINUE 

RETURN 
END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  THIS  ROUTINE  CALCULATES  THE  DIVISIONS  ON  THE  AXES 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


SUBROUTINE  NARKS ( INGSZ , PER, COORD ) 
DIMENSION  COORD(5) 
SIZE-PER*IHGSZ/4096. 
COORD(l)-SIZE/-2. 
COORD(2)-SIZE/-4. 

COORD( 3)-0 
COORD(4)-SIZE/4. 

COORD(5)-SIZE/2. 

RETURN 

END 


( 

1 


1 


